AD-A129  931  VIBRATION  CONTROL  OF  SPACE  STRUCTURES  VCOSS  A:  HIGH  AND  1/ 
LOH-AUTHORITV  HAR.  .  (U)  LOCKHEED  HISSILES  AND  SPACE  CO 
INC  SUNNVVALE  CA  J  N  AUBRUN  ET  AL.  JUL  82  LNSC-D883019 
UNCLASSIFIED  AFHAL-TR-83-2874  F22615-81-C-2228  F/G  28/11  NL 


V 


0^ 

AFWAL-TR-83-3074 

Vibration  Control  of  Space  Structures 
VCOSS  A:  High  and  Low-Authority 
Hardware  Implementations 


Lockheed  Missies  &  Space  Company,  Inc. 
1111  Lockheed  Way 
Sunnyvale,  Cafifbmia  94086 

July  1983 


Final  Report  Period 


Approved  for  pubic  release,  distribution  unfen ited 


IA  ' 


VS-  APR  1  l  1384 

u 

^  A 


u 


Flight  Dynamics  Laboratory 
Air  Force  Wright  Aeronautical  Laboratories  4 
Air  Force  Systems  Command  8 

Wright-Patterson  Air  Force  Base,  Ohio  45433 


04 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other 
than  in  connection  with  a  definitely  related  Government  procurement  operation,  the  United 
States  Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever;  and  the  fact 
that  the  government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


EROME  PEARSON,  Technical  Manager 
Engineering  Technology  Group 


HOWARD  F.  WOLFE,  Acting  Chief 
Structural  Vibration  Branch 


FOR  THE  COMMANDER: 


JAMES//.  OLSEN 

Assiifcint  for  Research  &  Technology 
Structures  &  Dynamics  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or  if  the 
addressee  is  no  longer  employed  by  your  organization  please  notify  AFWAL/FIBG  w_p  AFB, 
OH  45433  to  help  maintain  a  current  mailing  list". 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security  considerations, 
contractual  obligations,  or  notice  on  a  specific  document. 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  of  This  PAGE  (Whmn  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

UKFORE  COMPLETING  FORM 

1 

II 

3  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Submit*) 

Vibration  Control  of  Space  Structures  (VCOSS) 
VCOSS  A:  High-  and  Low-Authority  Hardware 
Implementations 

5  TYPE  OF  REPORT  6  PERIOD  COVERED 

final  Report 

15  Sep  1981  -  15  April  1983 

6  PERFORMING  ORG.  REPORT  NUMBER 

LMSC-D883019 

1  AUTHORS) 

J.  N.  Aubrun,  C.  Z.  Gregory  (ISI) ,  M.  G.  Lyons 
(ISI) ,  R.  L.  Kosut  (ISI),  A.  A.  Woods,  Jr. 

8  CONTRACT  OR  GRANT  NUMBERfs) 

F33615-81-C-3220 

9  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Integrated  Systems,  Inc. 

151  Universitv  Avenue.  Palo  Alto.  CA  94301 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

P.  E.  62201F 

Project  2401  Task  240104 

Lockheed  Missiles  &  Space  Co.,  Box  504,  Sunnyvale 

CA  W.U.  24010415 

1  1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Flight  Dynamics  Laboratory  (AFWAL/FIBG) 

Air  Force  Wright  Aeronautical  Laboratories 
Wright-Patterson  AFB  OH  45433 

12.  REPORT  DATE 

July  1983 

13.  NUMBER  OF  PAGES 

228 

14.  MONITORING  AGENCY  NAME  ft  ADDRESS^//  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (ot  this  report) 

UNCLASSIFIED 

1S».  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

N/A  _ 

16.  DISTRIBUTION  STATEMENT  (ol  thle  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  If  different  from  Report) 

osp L(  '  -/*.*  n o .  SL 

10.  SUPPLEMENTARY  NOTES 

AFWAL  Project  Manager;  Jerome  Pearson 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

— ■  -1  C'°deS  09  02  COMPUTERS 

22  02  SPACECRAFT  15  03  SURVEILLANCE  SYSTEMS 

20  11  SOLID  MECHANICS,  VIBRATION  13  13  STRUCTURAL  ENGINEERING 

20  06  OPTICS  09  05  FEEDBACK  AMPLIFIERS,  sfevO 

17  09  RADAR  ANTENNA  MECHANISMS 

20.  ABSTRACT  (Continue  on  reverse  aide  If  necessary  and  Identify  by  block  number) 

^The  report  considers  some  hardware  aspects  of  large  space  structure  control 
implementations.  Hardware  and  sy stern  performance  requirements  are  overviewed, 
principally  in  the  context  of  vibration  suppression.  Analytical  models  based 
on  soft  and  stiff  versions  of^CSDL  #2  structure  are  used  to  evaluate  vibration 
suppression  performance  for  actively  controlled  and  passive  structures.  The 
effects  of  hardware  bandwidth  limitations  and  system  modeling  errors  are 
evaluated  for  the  active  controls  case.  It  is  shewn  that  hardware  bandwidth 

DD  I  jAN*73  1473'  EDITION  OF  1  NOV  65  IS  OBSOLETE 

-UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dale  Entered) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  PAGEfH?>««i  Data  Entaratl) 


limitations  do  not  pose  significant  control  sign  problems.  The  basic 
elements  for  a  hardware  demonstration  program  to  evaluate  tliese  results 


is  given 


■ 


SECURITY  CLASSIFICATION  of  This  PAGEr»rh»n  Data  Entered) 


FOREWORD 


This  report  was  prepared  by  the  Lockheed  Missiles  and  Space  Company 
Integrated  Systems,  Inc.  The  work  was  done  under  Contract  F33615-81-C-3220 
for  the  Air  Force  Wright  Aeronautical  Laboratories,  Flight  Dynamics 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio.  The  research  was 
conducted  under  Project  2401  and  Task  240104,  "Vibration  Prediction  and 
Control,  Measurement  and  Analysis";  and  was  partially  funded  by  the  Defense 
Advanced  Research  Projects  Agency.  The  Air  Force  technical  monitor  was 
Jerome  Pearson  (AFWAL/FIB6) . 


M 


CONTENTS 


CONTROL  HARDWARE  ASSESSMENT 

1.  Disturbance  Source  and  Performance 
Metric  Characterization 

2.  Hardware  Requirements  Assessment 

3.  Existing /Proposed  Hardware  Devices 

4.  Required  Hardware  Development 
SYSTEM  DESIGN 

1.  Overall  Control  System  Design 

2.  Actuator  Selection  and  Modeling 

3.  Sensor  Selection  and  Modeling 

4.  Signal  Processing 
MODEL  INTEGRATION 

1.  VCOSS  Model  Description 

2.  Actuator /Sensor  Locations  and  Types 

3.  Open-Loop  System  Evaluation  for 
Control  Design 

4.  Control  Designs 
"Cautious"  Controller 

5.  Impact  of  Disturbance  Model  Changes 

6.  Low  Frequency  Actuator  Limitations 

7.  Effects  of  Actuator  Dynamics 

8.  Sensor  Modeling  Effects 

9.  Effects  of  General  Actuator /Sensor 
Modeling  Errors 

10.  Results  Comparison  for  Stiffness 
Designs 

11.  Conclusions 


Page 

1 

1 

5 

10 

12 

13 

13 

18 

34 

39 

47 

47 

61 

64 

82 

113 

107 

113 

125 

162 

165 

177 

179 


1 

1 


CONTENTS  (CONT.) 


Appendix 

A 


REFERENCES 


REDUCTION  OF  LARGE  FLEXIBLE  SPACECRAFT 
MODELS  USING  INTERNAL  BALANCING  THEORY 


B 


CHARLES  START  DRAPER  LABORATORY  - 
REFERENCE  DATA 


LIST  OF  ILLUSTRATIONS 


Ad 


Figure 


Disturbance  Factors 


Disturbance  Power  Spectral  Density 

General  Overview  of  Damping  Devices  Domain 
of  Operation 

System  Design  Methodology 

Classification  Matrix  of  Actuators  for  Vibration 
Control  and  Generic  Examples 

Inertial  Actuators  Range  of  Operation 

Typical  Region  of  Operation  of  CMG's 


Typical  Region  of  Operation  of  a  Proof-Mass 
Actuator 


CMG  Sizing 

PPM  Actuator  Dynamics  Model 

PPM  Actuator  Design -Performance  Region 

PPM  Actuator  Design  Curves 

Typical  Inertial-Grade  Accelerometer  Noise  PSD 

Wordlength  Effect  on  Coefficient  Perturbation 
and  Stability 

Minimum  Number  of  Bits  for  Stable  Control 
Computations 


Deleted 


An  Architecture  for  High-Speed  Control 
VCOSS  Model  Draper  Labs  Model  #2,  Revision  #3 


VCOSS  Finite -Element  Model 


Modified  Mirror  Models 


VCOSS  Optical  System 


Page 


LIST  OF  ILLUSTRATIONS  (CONT) 


§ 

%  Figure 

&  22 


?! 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


34 

35 

36 

37 

38 

39 


Page 


VCOSS  Evaluation  Model  Frequency  vs.  Mode 

Number  67 

PSD  of  Disturbance  to  LOSx  68 

PSD  of  Disturbance  to  LOSy  69 

Integrated  PSD  of  Disturbance  to  Total  LOS  70 

Integrated  PSD  of  Disturbance  to  LOSx  71 

Integrated  PSD  of  Disturbance  to  LOSy  72 

Open-Loop  Modal  Analysis  76 

Controllability /LOS  Levels  for  Actuators  1,  2,  and  3  79 

Controllability /LOS  Levels  for  Actuators  4,  5,  and  6  80 

Controllability /LOS  Levels  for  Actuators  7,  8,  and  9  81 

10-Mode  Control  Design  Selection  83 

//  VCOSS  (+)  Open-Loop  Spectrum,  (X.) 

Closed-Loop  Spectrum,  26th  Order  87 

//  VCOSS  Closed-Loop  Evaluation  88 

II  VCOSS  (+)  Open -Loop  Spectrum,  (X) 

Closed-Loop  Spectrum,  26th  Order  89 

//  VCOSS  (+)  Open-Loop  Spectrum,  (X) 

Closed-Loop  Spectrum,  26th  Order  90 

Closed-Loop  Modal  Costs  for  10-Mode  Controller  92 

Mode- Selection  for  14-Mode  Reduced  Model  93 

(+)  Open-Loop  Spectrum,  (X)  Closed-Loop 

Spectrum,  34th  Order  97 


(+)  Open-Loop  Spectrum,  (X)  Closed-Loop 
Spectrum ,  34th  Order 

(+)  Open-Loop  Spectrum,  (+)  Closed-Loop 
Spectrum,  34th  Order 

VCOSS  Open-Loop  (Open)  and  Closed-Loop 
(Shaded)  Modal  Costs  for  34th  Order  Controller 

(+)  Open-Loop  Spectrum,  (X)  Closed-Loop 
Spectrum,  HAC/LAC 


A 


LIST  OK  ILLUSTRATIONS  (CONT) 


igurc 

Page 

44 

(+)  Open-Loop  Spectrum,  (X)  Closed-Loop 

Spectrum,  HAC/LAC 

103 

45 

(+)  Open-Loop  Spectrum,  (X)  Closed-Loop 

Spectrum,  HAC/LAC 

104 

46 

VCOSS/PPM  Actuator  Power  Requirements 

115 

47 

VCOSS/PPM  Actuator  Power  Requirements 

116 

48 

VCOSS/PPM  Actuator  Power  Requirements 

117 

49 

VCOSS/PPM  Actuator  Power  Requirements 

118 

50 

Actuator  Sizing 

122 

51 

Low- Frequency  PPM  Design 

124 

52 

Frequency  Response  of  Active  and  Passive 

Dampers 

126 

53 

Open-Loop  Modal  Analysis  of  34-Mode  Model 

128 

54 

Open-Loop  Frequencies  of  34-Mode  Evaluation 

Model  and  Design  Model 

129 

55 

Full  Closed-Loop  Spectrum  (X=HAC,  0=HAC/LAC) 

134 

56 

Full  Closed-Loop  System:  X=HAC,  0=HAC/LAC 
(Expanded  Scale) 

135 

57 

Open-Loop  Modal  Coordinate  "Component  Costs" 
with  Corrected  Disturbance 

137 

58 

Open-Loop  Modal  Coordinate  "Component  Costs" 
with  Corrected  Disturbance  -  Expanded  Scale 

138 

59 

Full  Closed-Loop  Spectrum:  X=HAC,  0=HAC/LAC, 
*=Locus  as  First-Order  Rolloff  Decreases 

140 

60 

Effect  of  First-Order  Actuator  Rolloff  Upon  Mode  31 

141 

61 

Effect  of  First-Order  Actuator  Rolloff  Upon  Mode  32 

142 

63 

Effect  of  First-Order  Rolloff  Upon  Mode  34: 
0=HAC/LAC,  +  =  Open  Loop 

144 

64 

Performance  Degradation  Due  to  First-Order 

Actuator  Rolloff 

146 

65 

Effect  of  2nd-Order  Actuator  Dynamics  Upon  LAC 

150 

MST  OK  ILLUSTRATIONS  (CONT) 


Figure 

Page 

66 

Effect  of  2nd-Order  Actuator  Dynamics  Upon 

Mode  34 

151 

67 

Effect  of  2nd-Order  Actuator  Dynamics  Upon 

Mode  33 

152 

68 

Effect  of  2nd-Order  Actuator  Dynamics  Upon 

Mode  32 

153 

69 

Effect  of  2nd -Order  Rolloff  Upon  HAC/LAC: 

X=HAC ,  0=HAC/LAC ,  *  =  Locus  as  Rolloff 
Frequency  Decreases 

155 

70 

Effect  of  2nd-Order  Rolloff  Upon  HAC/LAC: 

Mode  31 

156 

71 

Effect  of  2nd-Order  Rolloff  Upon  HAC/LAC: 

Mode  32 

157 

72 

Effect  of  2nd-Order  Rolloff  Upon  HAC/LAC: 

Mode  33 

158 

73 

Effect  of  2nd-Order  Rolloff  Upon  HAC/LAC: 

Mode  34 

159 

74 

Performance  Degradation  Due  to  2nd-Order 

Actuator  Rolloff 

161 

75 

Accelerometer  Noise  Models 

163 

76 

Block  Diagram  of  Controlled  Spacecraft 

167 

77 

Stability  Margin  (Actuator  Bandwidth  3.7  Hz) 

175 

78 

Performance  Margin  10%  Degradation  (Actuator 
Bandwidth  160  Hz) 

176 

LIST  OF  ABBREVIATIONS 


CM  C, 

Control  Moment  Gyro 

CSDL 

Charles  Stark  Draper  Laboratory 

DOF 

Degrees  of  Freedom 

IIAC 

High  Authority  Control 

LAC 

Low  Authority  Control 

LOS 

Line  of  Sight 

LQG 

Linear  Quadratic  Gaussian 

LSS 

Large  Space  Structures 

PPM 

Pivoted  Proof  Mass 

PSD 

Power  Spectral  Density 

PZT 

Piezoelectric  Type 

RAM 

Random  Access  Memory 

RMS 

Root  Mean  Square 

RRI 

Riverside  Research  Institute 

RSS 

Root  Sum  Square 

SASS 

Six-Axis  Space  Sensor 

TAARA 

Three -Axis  Angular  Rate  and  Acceleration 

VCOSS 

Vibration  Control  of  Space  Structures 

xi 


LIST  OF  TABLES 


Table 


Sensor  Types  and  Applications 

Actuator  Types  and  Applications 

Characteristic  Parameters  of  Prototype 
PPM  Actuator 

Inertial  Grade  Space  Accelerometers 

Steady-State  LQG  Controller  Form 

Eigenvalue  Perturbation  and  Stability  Condition 

Minimum  Number  of  Bits  to  Ensure  Controller 
Stability 

Deleted 

Changes  in  Lumped  Masses 

Mirror  Unrestrained  Degrees  of  Freedom 

VCOSS  Mass  Matrix  Comparison 

VCOSS  Modal  Frequency  Comparison  (Hz) 

Line-of-Sight  (LOS)  Equations 

Results  of  Line-of-Sight  (LOS) 

Modal  Rankings  Suggested  by  the  Analysis 
Closed-Loop  Control  Spectrum  of  10-Mode  Model 
Closed-Loop  Filter  Spectrum  for  10-Mode  Model 
Closed-Loop  Control  Spectrum  for  14-Mode  Model 


Closed-Loop  Filter  Spectrum  for  14-Mode  Model 


LIST  OF  TABLES  (Cont) 


Table 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


Stochastic  Performance  Analysis  of  VCOSS 
Controllers 

Individual  Actuator  RMS  Effort  for  Each 
Controller  (N) 

Normalized  Stochastic  Performance  Analysis 
of  VCOSS  Controllers  -  6  Independent 
Disturbances 

Stochastic  Performance  of  VCOSS  Controllers- 
2  Independent  Disturbances 

RSS  Actuator  Effort  -  2  Independent  Distur¬ 
bances 

Normalized  Stochastic  Performance  of  VCOSS 
Controllers  -  2  Independent  Disturbances 

Corrected  Actuator  Effort  (N) 

Actuator  Sizing 

Actuator  Sizing 

Frequencies  of  10-Mode  Design  Model 

Control  Design  Spectrum 

Filter  Design  Spectrum 

34-Mode  Performance  Evaluation 

Full  Closed-Loop  LAC  Spectrum 

34-Mode  Performance  Evaluation 

Stability  and  Performance  Robustness  Margins 

Passine  System  Response  to  Complete  PSD 

Rev  03  Closed-Loop  Performance  Without  Discrete 
Frequency  Power  Spectra 


xni 


Page 

105 

106 

108 

110 

111 

112 

114 

119 

121 

130 
132 

131 
136 
149 
164 
171 
178 
178 


SECTION  I 


CONTROL  HARDWARE  ASSESSMENT 


1.  DISTURBANCE  SOURCE  AND  PERFORMANCE  METRIC  CHARACTERIZATION 

Typical  high  performance  large  flexible  spacecraft  of  the 
type  which  have  motivated  this  study  are  subject  to  a  wide  range 
of  disturbance  loads  characterized  in  Fig.  1.  The  selection  of 
control  architectures  and,  subsequently,  sensor /actuator ,  and 
processor  hardware  depends  strongly  on  the  precise  form  of  the 
system  performance  metric,  the  vehicle  structural  characteristics, 
and  the  relative  disturbance  levels.  Without  these  descriptions, 
hardware  selection  is  not  possible.  Many  performance  metrics  may 
be  selected,  depending  on  mission  requirements,  including  line-of- 
sight,  wavefront,  structural  amplitude,  and  rigid  body  pointing 
errors.  For  the  systems  of  interest  in  the  VCOSS  program,  only 
steady-state  dynamic  line-of -sight  (LOS)  errors  will  be  con¬ 
sidered.  This  type  of  error  will  form  the  basis  for  evaluating 
benefits  of  control-configured  spacecraft  design.  Limiting  the 
performance  metric  in  this  way  widens  the  choices  in  actuator  and 
sensor  selection  by  eliminating  any  direct  need  to  measure  or 
control  rigid  body  modes  or  absolute  motion  with  respect  to  iner¬ 
tial  space.  Static  LOS  offsets  due  to  shape  and  alignment  errors 
are  generally  controlled  with  specialized  actuator/sensor  systems 
that  operate  near  DC;  actuator/sensor  placement  is  governed  by 
spatial  frequency  error  and  optical  alignment  requirements  which 
are  not  generally  compatible  with  controllability /observability 
measures  for  dynamic  system  modes.  Thus,  it  is  usually  the  case 
that  dynamic  error  suppression  systems  must  be  synthesized  sepa¬ 
rately  from  (static)  figure  control  systems  and  the  distinction 
between  static  and  dynamic  figure  control  should  be  clearly  made 
in  any  discussion. 
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With  these  ideas  in  the  background,  it  is  now  possible  to 
qualitatively  evaluate  the  effects  of  the  disturbances  in  Fig.  1. 
For  structures  in  the  30m  size  range,  significant  modes 
below  0.1  Hz  are  unlikely.  Thus,  orbit  rate  based  disturbances 
induced  by  gravity-gradient  torques,  magnetic  dipole  moments,  and 
solar  pressure  effects  all  look  like  static  input  loads  to  a 
"stiff”  structure  and  will  be  compensated  by  static  figure  and 
alignment  subsystems.  These  disturbances  are  spectrally  separated 
from  those  which  interact  with  the  vehicle  dynamics  and,  thus, 
they  do  not  influence  requirements  for  dynamic  error  suppression 
systems  (i.e.,  these  disturbances  will  be  ignored  in  this  study). 
The  only  possible  exception  to  this  scenario  will  be  thermal  and 
solar  pressure  transients  during  solar  eclipse  with  the  earth. 

For  certain  applications  these  inputs  will  influence  the  required 
modal  damping  in  the  lower  frequency  modes.  Usually,  however, 
stabilization  of  the  attitude  control  system  places  similar 
damping  requirements  on  these  modes  so  that  transient  settling  at 
low  frequencies  is  generally  not  treated  as  a  separate  design 
problem. 

Vibration  sources  which  excite  structural  modes  may  then  be 
considered  to  originate  on-board  the  spacecraft  in  the  form  of 
cryo-cooler  pumps,  CMG  rotors,  reaction  wheels,  coolant  systems, 
and  reaction  jets.  Roughness  in  digital  controllers  used  to 
implement  stability  augmentation  and  disturbance  rejection  systems 
may  also  produce  unwanted  modal  excitation.  Momentum  storage 
devices  and  pumps  generally  produce  narrow-band  excitation  at 
multiple  frequencies  corresponding  to  rotor  and  bearing  retainer 
motion.  Coolant  flow,  reaction  jet  operation,  and  combustion 
processes  usually  produce  a  wide-band  uncorrelated  noise.  The 
exact  force  levels  these  types  of  excitation  produce  depend  on 
types  and  sizes  of  devices;  in  many  cases,  supporting  data  for 
broad-band  disturbance  levels  is  non-existent.  Meaningful  evalu¬ 
ation  of  structural  response  must  then  depend  on  point  design 
solutions  along  with  assumptions  of  system  linearity.  For  optical 
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systems,  some  questions  regarding  non-linearity  of  joints  in  the 
micro-vibration  region  have  been  raised  which  have  yet  to  be 
resolved  via  ground  or  flight  experiments  (necessary  to  remove 
the  Ig  pre-load).  For  the  present  evaluations  these  questions 
will  be  put  aside. 

These  considerations  can  be  invoked  to  create  a  strawman 
disturbance  spectrum  (developed  by  Riverside  Research  Institute 
(RRI)  and  Charles  Stark  Draper  Laboratory  (CSDL)  shown  in  Fig. 

2.  Draper  Model  Number  Two  had  two  sine  disturbances — one  at 
ION  and  10  Hz  on  the  top  of  the  truss  and  one  at  5  Hz  and  20N  on 
the  aft  body.  Because  of  the  discrete  nature  of  the  distur¬ 
bances  the  control  system  synthesis  addressed  only  part  of  a 
realistic  disturbance  rejection  design  problem.  In  addition, 
the  disturbances  were  oriented  along  the  stiff  axis  of  the  struc¬ 
ture.  To  avoid  these  problems  a  new  model  for  the  disturbances 
is  proposed  in  Fig.  2.  The  flat  Power  Spectral  Density  (PSD) 
was  selected  so  that  disturbances  occurred  across  the  many  modes. 
The  level  selected  was  chosen  on  the  basis  of  supplying  as  much 
disturbance  energy  to  the  system  as  the  original  two  discrete 
systems  supplied  if  a  damping  ratio  of  0.002  was  assumed.  The 
15  Hz  rolloff  was  selected  to  include  more  modes  than  before  and 

to  extend  disturbances  out  into  the  area  called  the  Low  Authority 

2 

Control  (LAC)  region.  The  1/f  rolloff  was  specified  to  ease 
the  total  bandwidth  and  keep  the  problem  in  hand.  The  three 
discrete  frequencies  are  selected  to  represent  motors,  etc.  The 
10  and  20  Hz  frequencies  were  selected  to  present  the  problem  of 
control  of  frequencies  at  multiples  of  each  other.  The  8  Hz  fre¬ 
quency  was  selected  to  provide  a  beat  frequency  against  the  10  Hz 
frequency  that  was  in  the  structural  control  bandwidth.  The  ±  1 
Hz  on  the  8  Hz  frequency  was  to  insure  a  mode  would  be  directly 
excited  by  a  sine  source.  Since  non-linear  elements  are  not 
included  in  the  models,  frequency  mixing  will  not  be  an  issue 
here. 
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The  controls  problem  posed  in  the  VCOSS  study  is  the  active 
rejection  of  this  disturbance  (Fig.  2)  as  measured  by  a  line-of- 
sight  performance  metric.  Other  measures  could  be  used,  such  as 
wave  front  error,  which  would  alter  the  form  of  the  synthesized 
control  systems.  These  issues  are  discussed  later  in  Section 
3  on  model  integration.  For  now,  the  assumptions  under  which 
actuator/sensor  hardware  is  to  be  selected  are  summarized  as  follows 

(1)  rigid  body  motions  are  controlled  by  other  systems 

(2)  static  figure  control  is  ignored 

(3)  only  system  dynamic  line-of-sight  error  (tilt)  is  to 
be  controlled. 

These  restrictions  now  permit  some  definite  selections  to  be 
made. 

2.  HARDWARE  REQUIREMENTS  ASSESSMENT 

In  this  section  we  attempt  to  specify  the  kinds  of  hardware 
necessary  to  address  the  problem  discussed  above.  In  subsection 

3,  generic  requirements  are  examined  with  respect  to  what  is 
actually  available. 

a.  Sensor  Systems 

Sensors  for  vibration  measurement  may  be  classified  into 
inertial  measurement  and  relative  measurement  categories.  The 
first  category  includes  accelerometers  and  instrument  gyros 
while  the  second  includes  all  optical  angular  and  distance  mea¬ 
suring  devices.  Additional  sensors  to  measure  point  source  dis¬ 
turbances  may  also  be  required.  Sensor  types  must  be  matched  to 
performance  requirements,  control  architectures,  and  mission  con¬ 
straints.  Some  of  these  considerations  are  listed  in  Table  1. 


TABLE  1  SENSOR  TYPES  AND  APPLICATIONS 


Sensor 

Type 

Measurement 

Class 

Applicable 

Control 

Architecture 

Comments 

Acceler¬ 

ometer 

Inertial 

Co-located/ 
non-colocated 
state  feed¬ 
back 

Servo-type  useful  > 

0.2  Hz  in  micro¬ 
vibration  region 

Rate 

gyros 

Inertial 

Co-located/ 
non-colocated 
state  feed¬ 
back 

Very  high  cost  for 

micro-vibration 

applications 

Load 

cells 

Relative- 
strain- 
guage  type 

Vibration 
load , 
actuator 
force 

measurement 

Temperature 

sensitivity /reliability 
problems 

Laser 

angular 

systems 

Laser 

phase 

systems 

Relative 

Non-< 

state 

back* 

colocated 
i  feed- 
► 

♦Colocated  mechaniza¬ 
tions  possible  with 
base  motion  compen¬ 
sation 

Multi¬ 

color 

interfer¬ 

ometers 

Absolute 

distance/ 

non-inertial 

Shape 

control 

Bandwidth  limited  for 
real-time  control  for 
scanning  systems 

Colocated  systems,  where  actuators  and  sensors  must  have  coinci¬ 
dent  effector/sensor  axes  and  locations,  are  most  easily  mecha¬ 
nized  using  inertial  measurement  instruments  or  optical  systems 
with  inertial  base  motion  compensation.  For  some  applications,  it 
has  been  demonstrated  experimentally  that  optical  angular 
measurement  can  provide  colocated  measurement  over  a  limited  fre¬ 
quency  range  (1).  Optical  systems  inherently  provide  the  widest 
bandwidth  (down  to  DC)  sensing  and  the  best  signal-to-noise  ratio 
of  any  system.  The  measurements  are  inherently  non-colocated , 
however,  and  computation  of  the  measurement  matrix  is  usually 
complex.  System  constraints  may  also  not  permit  deployment  of 
such  sensors  on  certain  parts  of  the  spacecraft. 

b.  Actuator  Systems 

Actuators  which  produce  constant  forces  and  moments  are 
necessary  for  attitude  control  but  not  for  vibration  control. 

Thus,  additional  actuators  capable  of  producing  only  oscillating 
loads  may  be  considered  for  the  VCOSS  application.  A  represen¬ 
tative  list  is  shown  in  Table  2  organized  along  the  lines  of  the 
previous  sensor  table. 

The  range  of  controllable  structural  amplitudes  as  a  function 
of  actuator  type  is  illustrated  in  Fig.  3.  Requirements 
imposed  by  the  current  VCOSS  strawman  configuration  and  distur¬ 
bances  limit  actuator  selection  to  a  subset  of  those  illustrated 
in  the  table.  Rejection  of  sinusoidal  and  broadband  disturbance 
to  produce  a  "quiet"  structure  may  require  actuator  effectiveness 
over  a  30  Hz  bandwidth  with  linear  well-defined  response  charac¬ 
teristics.  Reaction  wheels  and  reaction  jets  are  generally  not 
capable  of  this  type  of  response.  Member  dampers,  which  impose 
forces  and  moments  on  points  in  the  structure  are  generally  band¬ 
width  limited  by  the  additional  structure  necessary  to  inter¬ 
connect  two  nodes.  This  also  makes  actuator  modeling  more  diffi¬ 
cult.  It  appears  at  present  that  such  devices  are  more  suitable 
for  low-frequency  shape  control. 
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TABLE  2  ACTUATOR  TYPES  AND  APPLICATIONS 


Actuator 

Type 

Actuator 

Loads 

Applicable 

Control 

Architecture 

Comments 

Proof  - 
mass 

Inertial , 

AC  only 

Co 1 oca  ted /n on - 
colocated  state 
feedback 

Microvibration  control 
>_  1  Hz,  possible  non¬ 
linear  behavior 

CMG 

Inertial , 
limited 

DC  torque 

Attitude 
non-colc 
state  fe 

i  control , 

>cated 

;edback* 

♦Most  common  use,  high 
bandwidth  torque 
unavailable  with  cur¬ 
rent  designs,  very  high 
torque  possible 

Reaction 

wheel 

Inertial , 
limited 

DC  torque 

Low  bandwidth  ~  1  Hz , 
low  torque,  difficult 
vibration  isolation 

Member 

damper 

Intra- 
structural 
loads  (to 

DC) 

Non-colocated 
state  feedback 

Complex  control  distri¬ 
bution  matrix,  complex 
actuator  dynamics, 
usable  for  static  shape 
control 

Reaction 

jet 

Inertial 
to  DC 

Attitude  control, 
colocated/non- 
colocated  state 
feedback 

Inherent  vibration 
sources,  principal  use 
in  maneuvering,  and 
large  amplitude  control 

Figure  3.  General  Overview  of  Damping  Devices'  Domain  of  Operation 


3.  EXISTING /PROPOSED  HARDWARE  DEVICES 


The  kinds  of  hardware  shown  in  subsection  2  represent  either 
specialized  applications  of  existing  devices  or  entirely  now 
devices,  most  of  which  have  never  been  tested;  In  spacecraft 
system  design,  measurement  of  surface  modal  deflections,  and 
control  of  dynamic  mode  shapes  and  amplitudes  have  not  histori¬ 
cally  been  of  interest.  Hardware  development  has  not  been  spe¬ 
cifically  directed  toward  such  applications.  Consequently,  few 
directly  usable  devices  exist.  Paper  studies,  prototype  devel¬ 
opment,  and  modification  of  existing  hardware  have  occurred, 
nonetheless,  and  some  of  this  work  was  reported  to  DARPA  by  CSDL 
and  LMSC  in  a  joint  briefing  in  late  1980.  Very  little  new  work 
exists  to  date,  and  an  overview  of  this  device  survey  is  now  pre¬ 
sented  as  representative  of  current  thinking  in  the  hardware 
development  area  and  as  a  basis  for  the  device  models  suggested 
in  the  next  section. 

The  state  of  the  art  in  inertial  instrument  design  is  quite 
advanced  and  devices  which  are  now  flight  qualified  or  which 
could  be  flight  qualified  are  presented  in  Appendix  B.  The  list 
of  usable  servo-type  accelerometers  is  reasonably  long  and  instru¬ 
ments  similar  to  the  Bell  XI  are  listed  also  in  Section  II. 
Piezo-electric  type  accelerometers  do  not  currently  perform  well 
enough  for  the  VCOSS  application  because  of  their  high  noise 
level  and  poor  low  frequency  performance.  Two  new  instruments, 
derived  from  existing  CSDL  hardware,  the  Three-Axis  Angular  Rate 
&  Acceleration  (TAARA)  and  Six-Axis  Space  Sensor  (SASS)  sensors, 
provide  three  axis  angular  measurement  and  integrated  angular  and 
translational  measurement  respectively  in  a  single  instrument. 

Such  devices  can  be  colocated  with  a  cluster  of  Control  Moment 
Gyro  (CMG)  actuators  and/or  proof-mass  drivers.  The  character¬ 
istics  of  these  devices  are  described  in  Appendix  B.  A  new  low- 
cost  accelerometer  using  a  flexible  substrate  chip  technology  may 
have  possible  applications  to  VCOSS  if  sensitivity  and  signal -to- 
noise  can  be  improved  from  current  levels.  Prototype  devices  do 
show  some  promise  as  indicated  in  the  tables. 
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The  use  of'  optical  systems  to  measure'  structural  deflections 
was  originally  motivated  by  use  of  multi-color  interferometers  to 
compute  figure  error  on  large  surfaces  (such  as  mirrors).  Abso¬ 
lute  figure  measurement  is  unnecessary  for  dynamic  figure  error 
suppression  but  some  form  of  high-bandwidth  microvibration  rela¬ 
tive  velocity  and  position  measurement  will  be  needed  for  some 
applications.  Stand-off  measurement  capability  offered  by  these 
systems  is,  however,  difficult  to  use  because  the  effective 
measurement  matrix  usually  reflects  base  motion  of  the  sensor  as 
well  as  motion  of  the  desired  structural  components.  In  some 
instances,  an  accurate  measurement  matrix  can  be  constructed,  in 
others,  base  motion  may  be  restrained  inertially  (active  or  pas¬ 
sive),  and  for  some  applications,  base  motion  reference  can  be 
provided.  The  inertially  stabilized  retro-reflectors  shown  in 
the  charts  accomplish  this  last  objective  by  providing  a  refer¬ 
ence  surface  against  which  other  vibration  measurements  from  the 
optics  are  compared.  This  allows  colocated  measurement  and  con¬ 
trol  if  the  actuators  are  located  at  the  surface  measurement 
sites . 

A  representative  10  channel  microphase  optical  system  is 
then  described,  built  out  to  the  prototype  test  stage,  which 
might  be  used  in  conjunction  with  such  stabilized  retro-reflectors. 

Actuators  usable  for  VCOSS  systems  are  essentially  all  in  the 
early  development  stage.  The  basic  proof-mass  driver  using  elec¬ 
trodynamic  motor  technology,  is  illustrated.  A  proposed  high- 
bandwidth  CMG  which  uses  electrodynamic  gimbal  torquers  and  sus¬ 
pension  components  to  isolate  high  frequency  rotor  vibration  from 
the  controlled  structure  is  also  sketched.  No  prototypes  of  this 
device  currently  exist.  The  last  item  illustrated  is  a  sketch  of 
a  member  damper  concept  using  either  piezoelectric  or  electro¬ 
dynamic  drivers.  It  must  be  emphasized  that  such  devices  have 
only  been  tested  in  the  form  of  external  actuators  used  to  man¬ 
euver  test  specimens  and  have  not  actually  been  used  as  intra- 
structural  "dampers".  Small  piezoelectric  actuators  have  been 
used  as  in tras tructural  dampers  by  R.  L.  Forward  at  Hughes. 


4.  REQUIRED  HARDWARE  DEVELOPMENT 


Actuators  which  appear  most  promising  for  vibration  control 
include  proof-mass  (linear,  inertial  force)  and  control  moment 
gyro  types.  (These  devices  provide  band-limited,  inertial 
structural  loads.)  Characteristics  and  sizing  of  these  devices 
will  be  discussed  in  the  next  section.  Proof-mass  devices  exist 
only  as  laboratory  test  articles  and  must  be  developed  into 
flight  qualified  hardware.  CMG  flight  hardware  has  been  limited 
to  low-bandwidth  attitude  control  actuators  which  can  control 
low  frequency  bending  but  which  are  inadequate  to  affect  broad 
band  noise.  In  addition  to  current  designs,  new  high-bandwidth 
prototype  designs  need  to  be  brought  to  the  ground  test  stage  of 
development . 

Accelerometers  are  now  highly  developed  and  the  servo  types 
can  yield  nano-g  resolution  above  0.5  Hz.  Optical  systems  and 
instruments  are  well-developed  but  use  on  Large  Space  Structures 
(LSS)  will  be  very  case  dependent  and  sensitive  to  specific 
geometry.  System  level  tests  are  needed  to  bring  such  instru¬ 
ments  into  wide-spread  use  at  the  ground  test  level  before  con¬ 
figuring  them  for  subsequent  flight  tests. 

The  hardware  illustrated  in  the  previous  section  represents 
recent  thinking  by  CSDL,  LMSC ,  Integrated  Systems,  and  other 
industry  practitioners.  Many  of  these  devices  are  in  the  pre¬ 
development  or  early  development  stage.  Successful  deployment 
of  LSS  type  systems  in  the  next  5  to  10  years  will  depend  on  the 
rigorous  development  of  these  and  other  devices. 

In  Section  II  the  system  design  aspects  of  an  LSS  vibration 
control  design  are  explained,  including  hardware  selection  and 
placement  methodology.  Basic  models  for  selected  hardware  are 
described . 


SECTION  II 


SYSTEM  DESIGN 


1.  OVERALL  CONTROL  SYSTEM  DESIGN 

a.  Introduction 

This  section  will  attempt  to  define  a  rational  procedure  to 
specify  active  vibration  control  devices  and  integrate  them  in  the 
design  of  a  control  configured  space  structure,  given  certain 
performance  requirements  in  the  presence  of  on-board  or  environ¬ 
mental  disturbances.  However,  because  of  the  present  scarcity  of 
applicable  actuators  and  sensors  and  technological  limitations,  it 
is  expected  that  some  trade-off  studies  will  always  have  to  be 
made  between  actuator/sensor  performance  and  system  performance. 

It  is  not  unusual  that  a  small  relaxation  of  mission  requirements 
could  yield  a  significant  decrease  in  the  required  performance 
for  the  hardware  and  make  the  system  much  more  realizable.  Thus 
the  procedure  described  here  should  be  envisioned  more  as  a 
guide-line  or  as  one  step  in  a  more  complex  iterative  process. 


The  second  problem  encountered  in  the  design  of  a  Vibration 
Controlled  Space  Structure  is  that  the  starting  point  is  not  well 
defined  in  practice.  Although  it  would  seem  that  the  overall 
design  of  the  spacecraft  should  be  started  as  a  whole,  i.e., 
including  the  control  system  (the  so-called  Control  Configured 
Spacecraft),  it  is  more  usual  that  the  structural  elements  are 
first  defined  and  then  control  systems  are  built  around  them. 

This  second  approach  is  considered  in  this  report,  especially  in 
view  of  the  specific  application  to  the  model  spacecraft  developed 
by  CSDL  discussed  in  Section  III. 


.  System  Design  Main  Steps 


The  overall  process  is  sketched  in  Figure  4.  Once  the 
physical  structure  of  the  spacecraft  has  been  designed  to  sat¬ 
isfy  more  or  less  the  mission  requirements,  a  structural  (fin¬ 
ite-element,  FEM)  model  is  usually  obtained.  From  this  FEM  model,  the 
mathematical  model  for  the  distrubances ,  and  the  various  performance 
metrics  definitions,  a  state-space  model  is  obtained  which  will  be  the 
basis  of  the  open- loop  analyses. 

These  analyses  constitute  a  very  important  step  because 
it  is  from  them  that  actuator  and  sensor  selection  can  proceed. 

Various  methods  are  used  for  this  selection  (discussed  later) 
and  they  result  in  the  definition  of  the  location  of  the  actu¬ 
ator,  the  type  of  actuation  (force,  torque,  bending  moments, 
etc.)  and  the  direction  of  the  effort.  Similar  results  obtain 
for  the  sensors. 

At  this  point  the  set  of  actuators  and  sensors  are  still 
conceptual,  with  ideal  characteristics.  From  this  set,  and  the 
original  spacecraft  model,  the  system's  model  is  derived  to 
which  control  synthesis  techniques  are  applied. 

The  definition  of  the  control  laws  and  their  subsequent 
evaluation  furnish  information  concerning  the  control  effort 
and  the  sensing  sensitivity  which  generate  requirements  for 
the  actuators  and  sensors  to  be  designed.  If  these  require¬ 
ments  are  easy  to  meet,  the  process  could  stop  there.  In  prac¬ 
tice,  as  was  stated  earlier,  many  iterations  may  be  required. 

A  strong  driver  for  this  iterative  process  is  that  the  physical 
actuator  (or  sensor)  has  non-ideal  characteristics  (limited 
bandwidth,  nonlinear  behavior,  noise)  which  need  to  be  modeled 
and  evaluated  for  their  impact  upon  the  total  system  perfor¬ 
mance. 

This  last  process  of  actuator  design  could  vary  from  just 
selecting  an  "off-the-shelf"  item  to  designing  an  entirely  new 
actuator  or  sensor. 
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Finally,  the  control  design  leads  also  to  requirements  for 
the  digital  (and  analog)  hardware  which  in  turn  drive  the  choice 
or  design  of  the  microprocessors  and  electronic  compensators. 

c.  Actuator  Selection 

The  actuator  selection  process  is  generally  dependent  upon 
the  type  of  performance  which  is  sought.  In  general,  the  first 
step  is  to  "look"  at  the  structure  when  it  is  being  excited  open- 
loop,  and  determine  what  degrees  of  freedom  are  moving  the  most 
and  what  are  those  which  are  the  most  responsible  for  the  degrada¬ 
tion  of  the  performance.  This  "look"  can  be  either  graphical 
(movies,  computer  graphics,  displays,  etc.)  or  merely  mathematical 
in  nature  (controllability  maps,  etc.).  In  any  case  it  must  be 
done  with  a  keen  sense  of  the  physics  of  the  problem,  commonly 
known  as  engineering  judgment.  For  instance,  in  an  optical  sys¬ 
tem,  the  line-of-sight  (LOS)  error  will  be  a  direct  function  of 
the  motion  of  a  certain  critical  degrees-of- freedom  (dof)  of  the 
structure.  This  relationship  is  termed  the  performance  metric. 

If  the  problem  is  to  stabilize  the  LOS,  the  actuators  should  be 
placed  on  the  corresponding  nodes  and  act  in  the  direction  of  the 
excited  dof's.  This  method  leads  to  far  more  superior  results 
than  stuffing  the  structure  with  actuators  at  all  the  places  which 
seem  to  be  moving.  Certainly  the  latter  will  ensure  a  reduction 
of  the  vibration,  but  will  require  a  tremendous  amount  of  power 
and,  because  it  was  not  tailored  to  the  specific  problem,  the 
actual  LOS  performance  may  not  be  as  good.  However,  if  the  full 
surface  of  a  mirror  is  to  be  kept  within  optical  tolerances,  the 
same  performance  metric  method,  used  for  the  LOS,  may  indicate  the 
need  for  many  actuators  on  that  surface. 

The  more  serious  question  concerns  the  possibility  of  further 
reducing  the  number  of  the  actuators.  Despite  the  theoretical 
result  from  optimal  control  stating  that  if  all  the  modes  of 
interest  are  controllable  with  one  actuator  (this  just  implies 
that  those  modes  do  not  have  a  vibration  node  at  the  exact 


location  of  the  actuator)  then  only  one  actuator  is  enough  to 
do  the  .job,  there  are  several  other  considerations  which  will 
mandate  the  use  of  somewhat  more  "redundant"  set.  Between  the 
extremes  of  one  actuator  for  the  whole  system  and  as  many  actu¬ 
ators  as  there  are  degrees  of  freedom,  the  optimal  number  can¬ 
not  yet  be  determined  by  a  push-button  technique  but  rather 
by  trials,  errors  and  trade-off  studies. 

In  summary,  location  will  be  determined  by  the  performance 
metric  (i.e.,  at  the  nodes  which  have  a  direct  influence  on  the 
performance),  while  the  type  will  depend  upon  the  kind  of  motions 
involved.  Displacements  will  call  for  force  actuators,  rotations 
for  torque  actuators  (corresponding  to  proof-mass  actuators  or 
CMG's  for  instance).  But  again,  because  of  the  necessary  redun¬ 
dancy  and  non-unique  solution,  there  is  no  hard-line  choice  in 
general  but  only  trade-offs  at  the  system's  level.,  and  the  initial 
decision  may  be  modified  on  the  basis  of  the  next  step  of  the 
process,  where  requirements  are  defined  in  terms  of  power,  band¬ 
width,  etc. 

d.  Sensor  Selection 

Now  that  the  control  system  has  been  given  muscle,  the 
sensing  system  must  be  selected  so  that  the  loop  could  eventually 
be  closed.  Sensor  selection  is  more  difficult  than  its  actuator 
counterpart,  because  it  is  related  to  what  the  "error  signal"  is 
going  to  be,  which  is  part  of  the  control  design. 

Thus,  although  the  chart  of  Fig.  4  shows  this  process 
appearing  after  the  open-loop  analyses,  the  dotted  line  coming 
from  the  control  design  is  very  important.  From  open-loop  con¬ 
siderations,  one  may,  as  for  the  actuators,  place  sensors  where 
the  motion  is  maximum  and  choose  the  type  which  will  measure 
either  the  local  displacements  or  local  rotations.  In  an  optimal 
control  design,  the  goal  is  to  obtain  a  good  estimate  for  the 


state  vector,  and  observability  is  the  criterion.  However,  in 
some  special  cases,  specific  sensors  are  sought.  For  instance,  in 
the  case  of  on-board  sinusoidal  disturbance,  where  frequency¬ 
shaping  techniques  play  an  essential  role  in  the  control  design, 
it  is  important  to  measure  the  LOS  error,  because  it  is  a  direct 
input  to  the  control  system.  Certainly  the  LOS  error  can  be  esti¬ 
mated — reconstructed  from  nodal  motions  measurements — but  it  is 
very  advantageous  to  obtain  it  from  a  focal  plane  image  processor, 
for  instance,  or  any  other  direct  means. 

The  most  immediately  applicable  sensors  are  the  acceler¬ 
ometers  and  the  angular  and  linear  optical  sensors.  The  first 
give  inertial  measurements  and  are  straightforward,  but  the 
optical  sensors  give  only  relative  indications  and  thus  their 
implementation  is  less  obvious  because  it  involves  differential 
motions  between  parts  of  the  structure.  With  optical  sensors, 
the  idea  of  a  decentralized  colocated  system  has  to  be  abandoned 
or  generalized  to  bi-located  systems  which  do  not  lend  themselves 
to  simple  analytical  or  practical  methods  of  selection.  The 
process  of  selection  is  much  more  case  dependent  and  has  a  greater 
element  of  arbitrariness. 

In  summary,  location  and  type  of  sensors  are  determined  from 
both  open-loop  simulations  and  preconceived  ideas  on  the  control 
system  itself,  and,  like  actuators,  will  depend  on  the  type  of 
performance  sought. 

2.  Actuator  Selection  and  Modeling 

Structural  actuators  can  be  classified  to  reflect  the  duality 
between  translation/force  and  rotation/torque,  as  well  as  the 
criterion  which  distinguishes  where  the  reaction  forces  or  torques 
are  transmitted.  If  these  reactions  are  transmitted  back  into  the 
structure,  the  actuators  are  intrastructural ,  while  if  they  are 
transmitted  to  inertial  space,  the  actuators  are  inertial  .  The 
corresponding  classification  matrix,  with  some  generic  examples, 
are  shown  in  Figure  5. 
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Figure  5.  Classification  Matrix  of  Actuators  for  Vibration 
Control  and  Generic  Examples 


A  general  property  of  inertial  actuators  is  their  inability  to 
transmit  structural  vibration  energy  from  one  part  of  the  struc¬ 
ture  to  another,  as  can  occur  when  several  intrastructural  member 
dampers  are  used  simultaneously.  Nevertheless,  non-inertial 
dampers  (other  than  member  dampers)  may  be  required  for  very  low 
frequencies  or  global  reshaping  (e.g.,  curvature  control),  rather 
than  just  for  vibration  suppression. 


Low  frequency  and  large  amplitudes  are 
intrastructural  actuators,  while  lower  ampl 
quencies  will  use  inertial  devices.  Vibrat 
have  a  very  typical  range  of  operations  as  i 


The  main  design  drivers  are  the  travel  and  force  limits.  The 
shape  of  the  travel  limit  boundary  depends  upon  the  type  of 
actuator  considered  and  leads  to  various  types  of  optimization. 
Figures  7  and  8  show  typical  curves  for  CMG  and  Pivoted  Proof- 
Mass  (PPM)  actuators  respectively.  In  the  case  of  the  CMGs,  two 
more  constraints  were  considered,  the  maximum  gimbal  rate  (j  ) 
and  power  dissipation  (Pmax>.  The  output  force  limit  is  expressed 
in  terms  of  the  maximum  gimbal  acceleration  (o  )  that  the  gimbal 
motor  can  produce. 
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Figure  7.  Typical  Region  of  Operation  of  CMGs 
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The  boundaries  of  Fig.  8  are  derived  from  the  CMG  dynamic 
equations  and  the  equations  describing  the  structural  dynamics. 

For  the  purposes  of  illustration,  it  is  convenient  to  use  a 
simple  example.  The  linearized  equations  of  motion  for  a  single 
mode  controlled  by  a  single  CMG  are  as  follows: 


Jo  =  K(o  -  o)  +  H§ 
c 


8  +  2(  m  6  +  w  e  =  *  Ho  (2) 

n  n  n  n 

where  0  is  the  local  rotation  angle  associated  with  the  mode  and 
a  is  the  CMG  gimbal  angle;  J  and  H  are,  respectively,  the 
transverse  moment  of  inertia  and  angular  momentum  of  the  CMG;  K 
is  the  servo  loop  gain  and  a  the  commanded  rate;  wn,en,  and 
are,  respectively,  the  frequency,  damping  ratio,  and  influence 
coefficient  of  the  mode. 


Within  the  bandwidth  of  the  CMG  actuator,  o  -o  ,  and  for 

c 

w  *=  o»n,  equation  (2)  leads  to  the  relation 
*  H 

5  -  2TV-  *  <3' 

n  n 

The  various  hardware  constraints  can  now  be  expressed  in 
terms  of  inequalities  involving  u>  and  9.  For  instance: 


4>  H  a 

•  ,  •  — .  .  .  n  max 

o  <  a  =>  6  <  rr- -  - 
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,  etc. 
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Similar  equations  may  be  derived  for  proof-mass  actuators, 
but  in  this  case  instead  of  a  torque  (Her)  ,  the  actuator  produces 
a  force  equal  to  mx,  where  m  is  the  mass  of  the  proof -mass  and 
x  its  displacement.  Equation  (3)  now  has  the  form: 


2c  u> 
n  n 


where  6  is  the  local  displacement  of  the  structure  (where  the 
actuator  is  attached).  Thus  the  boundaries  shown  in  Fig.  8. 
Absolute  bandwidth  is  determined  from  the  electrical  properties  of 
the  system  (e.g.,  AC  coupling,  windings  inductance,  etc.). 

Design  (scaling)  curves  may  be  obtained  from  simplified 
actuator  models  to  help  properly  size  these  devices.  Examples  are 
derived  below  for  the  CMGs  and  the  PPMs. 

a.  CMGs 

The  model  equations  are: 


where  T^  is  the  gimbal  torque  (produced  by  the  gimbal  motor)  and 
Tv  the  torque  produced  on  the  vehicle.  The  latter  is  limited  by 
the  mechanical  strength  of  the  rotor,  which  is  directly  related  to 
the  wheel  cross-section  at  the  shaft.  If  CMGs  of  different  size 
are  considered,  assuming  that  they  are  geometrically  proportional, 
and  letting  L  be  this  geometric  scaling  factor,  Ty  will  then 
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vary  as  L2 .  Similarly  the  other  properties  of  the  CMG  may  be 
written  in  terms  of  L.  Based  on  simple  physical  considerations 
it  is  possible  to  establish  such  scaling  factors  (the  "K” 
quantities  in  the  following  expressions): 

5 

J  =  KjL  (spin  inertia) 

4 

H  =  KjjL  (momentum) 

2 

Ty  =  (output  torque) 

2 

=  KnL  (gimbal  torque) 


Other  quantities  may  be  derived  from  these  scaling  laws,  such  as 

a  or  a  which  are  more  currently  used  by  gyro  mann- 

ma  x  max 

f ac  turers : 


°ma*  =  <W  L 


o  =  (K„/K  .)  L 
max  G '  J 


Assuming  that  the  maximum  angular  travel  o  is  used,  and 

m&  a 

a  sinusoidal  excitation,  (4)  and  (5)  yield: 

A  O 

T„  =  Ho  a)  =  K„o  L  u>  <  KtL 
V  max  H  max  -  T 


_  T  2  «  I  ^  ^  ^  t  ^ 

G  max  J  max  — 


Thus  the  conditions: 


L  “  <  <VKH> 


A2  <  (K^Kj)  =m„ 


A  nomogram  such  as  the  one  shown  in  Fig.  9  may  then  be  used 


for  sizing. 


operating 

frequency 


Figure  9.  CMG  Sizing 


PPMs 


1) . 

Similar  curves  may  be  obtained  for  the  PPMs.  The  dynamic 
model  shown  in  Fig.  10  leads  to  the  following  linearized 
equations : 


Varying  d  above  or  below  this  value  will  change  the  sign  of  the 
output  torque  as  well  as  its  magnitude.  The  mass  of  the  PPM 
actuator  and  its  CM  may  be  adjusted  by  changing  the  position  and 
mass  of  a  lead  piece  situated  at  the  distance  l  from  the  pivot 
point. 


Let  bQ  ,  mQ  and  I  be  respectively  the  position  of  the 
CM,  the  mass  and  the  inertia  with  respect  to  the  pivot  of  the 
"unloaded"  actuator,  and  the  mass  of  the  lead  piece  with 

which  it  is  loaded.  Then 


=  m  + 

m. 
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1 

m  b 
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(10) 


and  using  Eqs.  (6)  and  (8)  shows  that 
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(11) 


There  are  two  important  parameters  in  the  design  and  use  of 
PPM  actuators  besides  the  mass/inertia  properties  mentioned 
previously: 


(1)  f  =  maximum  electrodynamic  force 
em 


(2) 


a 

m 


maximum  angular  displacement 


Since  the  angular  displacement  a  depends  upon  the  frequency, 
optimal  choices  for  the  design  parameters  will  depend  upon  the 
frequency  range  of  application.  From  the  two  conditions 
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and  Eqs.  (6), 

(10) 

and  (11),  two  new  conditions  are  found: 
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when  u)/2u  is  the  frequency  of  the  output  force. 

These  two  limits  are  shown  on  Fig.  11  as  functions  of  i/r 
for  a  prototype  PPM  actuator  whose  characteristic  parameters  are 
shown  in  Table  3. 


TABLE  3  CHARACTERISTIC  PARAMETERS  OF  PROTOTYPE  PPM  ACTUATOR 


Mass  Geometry  Electrodynamic 


I  =1.5  10"5 
o 

.  2 
kgm 

b  =  .016 
o 

m 

m  =  .088 

o 

kg 

r  =  .021 

m 

f  =  1  N 

em 

e  =  .  067 

rad 
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Figure  11.  PPM  Actuator  Design  -  Performance  Region  (fj) 


An  optimal  value  for  £/r  is  obtained  at  the  intersection  of 
the  two  curves.  Assuming  that  this  optimal  value  is  always  used 
in  the  design,  it  is  then  possible  to  proceed  with  design  curves 
by  relating  the  various  parameters  to  the  geometric  scaling  factor 
L.  Fig.  12  shows  a  set  of  such  curves  for  various  values  of  the 
ratio  x  =  mi/m  •  The  value  L  =  1  corresponds  to  the  existing 
prototype  actuator.  For  desired  values  of  f  and  o> ,  &  point  is 

chosen  on  the  right  vertical  axis.  Intersecting  the  corresponding 
horizontal  line  with  the  first  set  of  curves  ascending  gives 
values  for  L.  Then  the  optimal  £/r  ratio  is  obtained  with  the 
other  set  of  curves  (descending). 

c.  Actuator  Models 

(1)  CMG  Model 

This  model  is  linear  and  includes  the  interaction  between  the 
structure  and  the  CMG.  Defining  the  gimbal  axis  by  the  unit 
vector  g  and  the  nominal  position  of  the  momentum  vector  by  the 
unit  vector  h  ,  the  following  equations  hold: 

Jo  =  T  +  Hc*w  (14) 

g 

Ty  -  -co  (15) 

where  c  =  i  x  h  and  u  is  the  local  angular  rate  of  the  struc¬ 
tural  element  to  which  the  gyro  is  attached. 

(2)  PPM  Model 


This  linear  model  is  as  follows: 


PPM  Actuator  Design  Curves  (x  =  /m  as  a  Paramete 


Assuming  that  the  application  point  is  chosen  according  to 
Eq.  (9)  (no  torque  produced)  the  PPM  equations  may  be  written  as: 

mx  +  k(x-6)  *  f  (16) 

fv  *  -*  (17) 

where  k  is  the  ilex  pivot  spring  constant,  x  the  inertial 
displacement  of  the  center  of  mass  of  the  PPM  actuator  moving  arm 
6  the  local  displacement  of  the  structure,  and  f  the  reaction 
force  seen  by  the  structure;  f  is  related  to  the  electrodynamic 
force  f  by  the  equation  (11). 


3.  SENSOR  SELECTION  AND  MODELING 


The  primary  sensors  selected  for  the  VCOSS  design  study  are 
accelerometers.  These  devices  provide  the  most  practical  and 
lowest  cost  sensing  alternative  for  vibration  rejection  and  are 
relatively  easy  to  model.  Optical  systems  used  outside  the 
laboratory  environment  will  have  to  deal  with  complex  geometries 
base  motion  compensation,  and  other  issues  too  complex  to  model 
in  the  present  study.  Trade  studies  conducted  in  Phase  II  with 
these  basic  sensors  will  provide  a  foundation  for  studying  more 
complex  mechanizations  later  in  the  program. 

a.  Background 

Accelerometers  of  the  piezoelectric  type  (PZT)  have  been 
extensively  used  for  structural  modal  testing  and,  recently,  for 
feedback  control  associated  with  large  space  structure  stability 
augmentation.  Although  PZT  mechanizations  offer  rugged,  low- 
cost  instrumentation,  high  noise  PSD  at  low  frequencies  may  make 
them  unsuitable  for  LSS  control  applications.  For  system  iden 
tification,  where  ensemble  averaging  and  advanced  signal  pro¬ 
cessing  techniques  can  improve  signal-to-noise  ratios,  the 
PZT  devices  may  still  be  adequate,  though  for  real-time  control 
between  0.5  and  10  Hz  these  devices  may  be  excessively  noisy. 
Inertial  grade  servo  force  rebalance  and  PIGA  (pendulus  inte¬ 
grating  gyro  accelerometer)  types  can  provide  substantially 
higher  performance  at  lower  frequencies  if  one  to  two  orders 
of  magnitude  cost  increase  can  be  justified  over  the  laboratory 
grade  PZT. 

b.  Typical  Inertial  Grade  Accelerometer  Performance 

The  noise  PSD  for  a  generic  inertial  grade  accelerometer  is 
shown  in  Figure  13.  This  behavior  is  representative  of  all  high 
performance  accelerometers  such  as  the  Bell  model  7  and  others 
listed  in  the  CSDL  data  of  Table  1.  Above  0.5  Hz  (3  rad /sec) 
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the  noise  spectrum  is  essentially  white  with  an  RMS  acceleration 
resolution  on  the  order  of  10'8  m/sec2.  Below  0.5  Hz  the  response 
is  dominated  by  1/f  noise  down  to  about  10-3  rad/sec.  Below  this 
frequency,  bias  and  drift  stability  issues  dominate  the 
performance,  which  varies  from  instrument  to  instrument,  and  such 
values  are  difficult  to  measure.  For  LSS  control,  the  region 
above  0.5  Hz  is  probably  the  most  practical  one. 


accelerometer 
output  \ 


2  4 

m  /sec  /rad/sec 


10-4  to  10  7  g  bias/stability  depending 
on  instrument 


1/f  noise 


band-limited  white 
/spectrum 


(see  text) 
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LSS  control  range 


3  rad/sec 


Figure  13.  Typical  Inertial-grade  Accelerometer  Noise  PSD 


TABLE  4  INERTIAL  GRADE  SPACE  ACCELEROMETERS 


SyttfOft  4841  S«rvo  fore#  balance  1.5  ■  1  *  1.0  200  1.5  p*«k  2  i»9  Delta,  Trtd«nt 

Oo*n»r  pandulttfl  I  h,  (0  Sub,  Lance, 

Titan  UK  0. 


For  the  white  portion  of  the  spectrum,  the  spectral 
intensity,  R,  is  typically  (BW  denotes  bandwidth) 

02  «  R(BW) 

(10-8  m/sec2)2  =  R  (0.1  rad/sec) 

R  =  10~15(m  sec2  )2 /rad/sec  . 


Thus,  over  a  100  rad/sec  typical  useful  operating  range  for  LSS 
control  (above  0.5  Hz),  the  RMS  noise  observed  is 

o2  =  [10-15  m/ sec2  see-1]  100  rad/sec 

o  =  3  x  10~7  m/sec2  =  30  ng 

Values  as  low  as  0.3  ng  may  be  achievable  with  instrument  modifi¬ 
cations  and  greatly  increased  costs.  Cost  of  "standard"  instru¬ 
ments  should  range  from  $5K  to  $15K  which  probably  does  not 
include  space  qualification  and  certainly  does  not  include  special 
modifications. 
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For  space  structure  control  including  vibration  rejection  on 
precision  optical  surfaces,  figure  error  at  near  IR  is 

X  2  7 

|Ay|  =  Displacement  =  20  =  2Q~  =  0,135  um  . 

For  sinusoidal  deflection  of  the  surface, 

Ay  =  | Ay |  e^ut 

and  |  Ay  |  =  u>2  |  Ay  | 

=  (0.135  pm)  w2 

at  b  =  0.5  Hz  (  3  rad/sec) 

| Ay |  =  lpm/sec2  =  0.1  pg  =  100  ng. 

Thus,  precision  optical  surfaces  should  be  controllable  to  X/20 
tolerances  down  to  0.5  Hz  with  a  signal  to  noise  ratio  of  better 
than  3  to  1  at  low  frequencies.  This  S /N  improves  rapidly,  as  the 
frequency  squared,  at  higher  frequencies.  The  degraded  signal  to 
noise  at  low  frequency  will  probably  not  have  a  significant  impact 
on  system  performance  due  to  the  overlap  between  figure  control 
and  vibration  control  systems  at  low  frequency.  Figure  sensing 
systems  will  control  deformation  below  0.5  Hz  using  absolute 
position  measurement  with  S /N  >  103 .  For  non-precision  structural 
deformations,  control  tolerances  can  be  relaxed  by  at  least  an 
order  of  magnitude  which  will  yield  an  S /N  30. 

d.  Preliminary  Conclusions 

Although  microvibration  control  will  be  difficult  with  PZT 
type  instruments,  high  performance  should  be  attainable  with 
currently  available  inertial  grade  accelerometers,  especially 
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those  designed  lor  grad i omot or  applications.  Signal  processing 
electronics  may  require  some  redesign  but  noise  levels  well  be  low- 
50  ng  should  be  achievable  in  the  frequency  range  0.5  -  lb  He. 
Commercial  versions,  possibly  already  flight  qualifiable,  should 
cost  between  $5K  and  $15K.  Ultra-high  performance  types  such  as 
Bell  Mesa  might  cost  $100K/unit.  Sensing  below  0.5  Hz  will  be 
handled  by  the  figure/alignment  control  system  sensor  which, 
typically,  will  be  an  optical  measurement  instrument. 

In  the  frequency  range  of  interest,  sensor  dynamics  are  not 
of  any  particular  concern  since  actuator  dynamics  will  dominate 
system  behavior  much  earlier  in  the  spectrum.  These  devices  can 
be  modeled  effectively  using  only  bias,  scale  factor,  and  white 
additive  measurement  noise  terms. 

4.  SIGNAL  PROCESSING 

Difficulties  experienced  in  VCOSS  Phase  IA  with  digital 
implementation  have  led  to  new  signal  processing  architectures  to 
achieve  high  sample  rates  for  multivariable  controllers.  New 
architectures,  discussed  near  the  end  of  this  section,  will  remove 
sampling  rate  constraints,  but  in  so  doing,  raise  the  issue  of 
filter  stability  vs.  wordlength.  These  issues  are  discussed 
first.  The  basic  signal  processing  issues  for  LSS  control  are: 

•  Wordlength  and  algorithm  selection  for 

-  Minimum  pole  shift  (robustness) 

-  Mechanization  stability 

•  Sample  rate  selection 

-  Constraints  imposed  by  stability 

-  Constraints  imposed  by  architectures 

-  Architecture  assessment/memory  requirements 

•  Requirement  for  anti-aliasing  filters 
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•  Other  investigations 

-  Quantization  effects  and  controller  roughness 

-  Throughput  capacity  and  efficiency 

-  Flexibility,  fault-tolerance 

-  Emulation  of  large  spacecraft  model 

Some  of  these  issues  are  discussed  below, 
a.  Multivariable  Flight  Control 

Multivariable  control  methods  have  been  highly  successful  for 
designing  complex  control  systems  for  multi-input/multi-output 
systems.  To  apply  multivariable  control  methods  to  flight  con¬ 
trol,  the  spacecraft  is  generally  represented  by  the  equations  of 
Table  3. 

Multivariable  control  design  is  based  on  the  quadratic  cost 
functional  shown  in  the  table  which  produces  the  closed-loop  feed¬ 
back  control  law. 

When  the  state  x  is  not  measured,  x  must  be  estimated 
from  sensor  outputs  using  the  Kalman  filter  also  shown  in  the 
table. 

In  the  steady-state,  the  control  and  filter  gains  are 
constant  and  can  be  solved  off-line  for  C  and  K.  C  and  K 
can  then  be  used  to  compute  u  and  x  in  real  time.  Therefore, 
the  real  time  calculations  required  to  implement  a  steady-state 
multivariable  controller  are  given  in  Table  3. 

The  real  time  implementation  can  be  performed  by  analog  or 
digital  circuits.  However,  the  digital  implementation  is  often 
more  flexible  since  it  can  be  easily  modified. 
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By  examining  the  steady-state  controller  equations  it  is 
clear  that  only  matrix-vector  multiplication  is  required  for 
implementation.  The  matrix  elements,  which  are  constant  in  the 
steady-state  designs,  can  be  stored  in  memory  for  each  flight 
condition.  These  observations  form  the  basis  for  the  architecture 
described  later. 

b.  Stability  Analysis 

Because  of  modeling  uncertainties,  it  is  important  to  examine 
the  margin  of  stability  associated  with  the  computations.  Figure 
14  gives  a  geometric  interpretation  of  a  stability  margin  in  the 
z-plane.  |z|  represents  the  magnitude  of  poles  of  H(z)  and 
| A z |  is  the  distance  from  |z|  to  the  edge  of  the  unit  circle. 

The  magnitude  of  | A z |  is  dependent  on  a  number  of  factors  such 
as  modeling  uncertainties  and  finite  wordlength  considerations. 
From  an  analysis  viewpoint,  it  is  convenient  to  lump  these  con¬ 
siderations  into  a  perturbation  of  Fjj  and  analyze  the  stability 
of  the  eigenvalues  of  i  +  AFu» 

If  FX1  is  perturbed  to  Fj  j  +  AFni  then  the  corresponding 
change  in  the  eigenvalues  of  Pj j ,  can  be  estimated  from  the 
formula  shown  in  Table  6. 

c.  Wordlength  Considerations 

The  formula  in  Table  7  represents  an  effective  numerical  test 

for  rapidly  determining  the  minimum  number  of  bits  required  to 

implement  a  given  digital  flight  control  law.  The  magnitude  of 

the  maximum  eigenvalue  of  Fj j  is  the  parameter  which  defines  the 

margin  of  stability.  The  equation  was  evaluated  for  different 

values  of  A  and  n.  The  results  are  shown  in  Figure  15. 

max 

Note  that  the  number  of  bits  required  for  stability  increases 
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TABLE  7  MINIMUM  NUMBER  OF  BITS  TO  ENSURE  CONTROLLER  STABILITY 
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This  formula  represents  an  effective  numerical  test  for 
rapidly  determining  the  minimum  number  of  bits  required  to  imple¬ 
ment  a  given  digital  flight  control  law.  The  magnitude  of  the 
maximum  eigenvalues  of  Fii  is  the  parameter  which  defines  the 

margin  of  stability  =  1  -  X  =  AX.  The  equation  was  evaluated 

max 

for  different  values  of  X  and  n.  The  results  are  shown  in 

max 

Figure  15.  Note  that  the  number  of  bits  required  for  stability 

increases  substantially  as  ♦  1  and  as  n  increases. 

max 

Because  X^^  is  directly  related  to  the  sampling  rate  (1/At), 
Figure  15  indicates,  therefore,  that  at  high  sampling  rates,  a 
large  number  of  bits  are  needed  to  keep  the  computations  stable. 
Since  the  improvement  in  stability  is  marginal  beyond  16  bits,  a 
16  bit  arithmetic  processor  is  used. 

d.  Architecture  of  the  Controller 

The  architecture  of  the  processor  is  centered  around  a  high¬ 
speed  matrix-vector  multiplier  (see  Figure  17).  The  control 
design  matrix  elements  are  stored  in  data  memory  and  then  down¬ 
loaded  into  a  high-speed  RAM.  The  analog  measurement  data  are 
converted  to  digital  form  by  a  high-speed,  12-bit,  A/D  converter. 
The  results  are  stored  in  a  high-speed  RAM  along  with  the  state 
estimates.  The  control  law  is  then  evaluated  by  multiplying  the 
control-design  matrix  by  the  state  estimates  and  measured  data  in 
RAM.  A  high-speed,  16-bit  hardware  multiplier /adder  performs  the 
matrix  multiplications  and  additions  required.  All  addresses  for 
the  RAMs  and  hardware  multiplier/adders  are  generated  by  the  firm¬ 
ware.  The  digital  controls  are  then  sent  to  a  12-bit  D/A  con¬ 
verter.  The  D/A  converts  the  digital  controls  to  analog  signals 
which  are  stored  in  sample-and-hold  circuitry.  Once  the  controls 
have  been  updated,  the  controls  are  simultaneously  sent  to  the 
actuators . 

In  phase  II  of  the  VCOSS  program,  digital  emulations  like  the 
illustration  above  will  be  used  to  assess  signal  processing  hard¬ 
ware  effects  on  the  CSDL  #2  system. 


SECTION  III 


MODEL  INTEGRATION 

1.  VCOSS  MODEL  DESCRIPTION 

This  section  describes  the  evolution  of  the  VCOSS  systems 
model  from  the  original  CSDL  data.  Development  of  complex 
control  design  and  evaluation  models  has  historically  been  a 
team  effort  involving  many  validation  processes  to  assure  that 
structural  data  and  performance  matrices  correctly  represent 
the  physics  of  the  problem.  This  model  development  process  has 
been  followed  for  VCOSS  and  some  of  the  modeling  revisions  are 
discussed  now.  This  validation  process  is  essential  to  the  syn¬ 
thesis  of  meaningful  and  physically  implementable  controls  and 
to  the  accurate  assessment  of  control  system/vehicle  performance. 

The  structural  model  used  for  the  VCOSS  control  design  is 
the  Charles  Stark  Draper  Laboratory  (CSDL)  model  tt 2  (revision  #3) 
(Figures  18  and  19).  Model  #2  was  described  in  CSDL  report 
C-5437  dated  September  1981  (ref.  7)  and  revision  #3  was  des¬ 
cribed  in  CSDL  memo  entitled  "VCOSS  Design  Model  and  Disturbances" 
dated  November  12,  1981  (ref.  8). 

Listed  here  (in  items  la  -  h)  are  additional  points  of 
information  that  are  needed  to  create  the  VCOSS  revision  tt 3  model 
which  are  not  covered  by  either  of  the  aforementioned  CSDL  docu¬ 
ments. 
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a.  CSDL  model  #2  (ref.  7)  had  lumped  masses  to  provide  the 
total  mass  distribution.  The  VCOSS  model  (ref.  8)  uses  both 
lumped  masses  and  structural  member  weights  to  define  the  total 
mass  distribution.  The  mass  changes  needed  for  reference  7 
model  to  become  reference  8  model  are  shown  in  Table  9.  Essen¬ 
tially,  these  changes  represent  changes  in  the  mirrors  and  their 
supports. 


TABLE  9  CHANGES  IN  LUMPED  MASSES 


Node 

Number 


Translation  (x,y,z) 
Mass  (kg.) 


(ex,ey,6z) 

Rotation 

Mass 

(kg-M2) 


27 

69.5 

28 

6.74 

29 

69.5 

30 

6.74 

32 

6.74 

33 

6.74 

34 

69.5 

35 

69.5 

9 

67.4 

10 

67.4 

11 

67.4 

12 

67.4 

44 

3500. 

20611,10500, 

28777 

48 

81.91 

270 

50 

163.82 

540 

52 

73.82 

270 

53 

73.82 

270 

55 

163.82 

540 

57 

81.91 

270 

3 

b.  Material  properties  now  include  a  value  of  P  =  1720  kg/m 
(ref.  8)  for  all  beam  members  except  the  tubular  truss  members 
#116  and  #185  which  use  a  value  of  p=  2579.701  kg/m^  (ref.  8). 

c.  The  value  of  torsional  inertia  (J)  for  the  beam  members 
is  equal  to  twice  the  bending  moment  of  inertia  (J  =  21). 

d.  The  equipment  section  was  separated  30  centimeters  from  the 
truss  structure.  This  means  that  node  numbers  43,  44,  45,  48, 

49,  50,  51,  52,  53,  54,  55,  56,  57  have  a  z  coordinate  of  -1.30  m. 

e.  The  changes  listed  in  reference  8  for  the  structural  beams 
were  listed  by  beam  member  numbers  rather  than  by  the  node  num¬ 
bers  for  the  ends  of  the  beam. 

f.  Section  properties  were  shown  in  reference  8  for  each  of  the 
beam  members.  Many  of  the  same  section  properties  were  repeated 
for  different  beam  members.  In  order  to  make  a  concise  section 
property  table,  a  unique  set  of  these  section  properties  (20  in 
all)  was  constructed  and  assigned  section  numbers  for  each  beam 
member  to  reference  more  efficiently. 

g.  The  new  mirror  models  (Figure  20)  described  in  reference  8 
assume  that  each  surface  is  a  planar  rigid  body  which  is  connec¬ 
ted  to  the  support  structure  by  kinematic  mounts  which  decouple 
elastic  motion.  The  CSDL  NASTRAN  model  uses  rigid  body  kine¬ 
matic  equations  to  relate  the  motions  of  the  mirror  center  of 
mass  to  its  support  restraints.  The  restraint  conditions  are 
shown  on  the  sketches  for  each  mirror  model  in  that  memo.  LMSC 
ASTRO  model  uses  very  stiff  massless  tubes  (outer  diameter  =  1  M, 
inner  diameter  =  .9  M)  to  connect  the  mirror  center  of  mass  to 
its  support  points.  Guyan  reduction  was  used  on  these  beam 
element  stiffness  matrices  to  eliminate  the  mirror  support 
degrees  of  freedom  which  are  not  attached  to  the  backup  struc¬ 
ture. 
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The  unattached  d.o.f.  which  provide  a  "free"  motion  for 
the  mirrors  are  listed  below. 


TABLE  10  MIRROR  UNRESTRAINED  DEGREES  OF  FREEDOM 


Mirror 


Unrestrained 
Degrees  of  Freedom 


PRIMARY 


TERTIARY 


SECONDARY 


34 

35 
2830 

27 

29 

3233 

910 

1112 

40 

11 
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x  &  0. 
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Vez 

Vez 
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y  z 
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e  ,0 
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0  ,0 
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0  ,6 
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FOCAL 

PLANE 


9 

40 


h.  The  LOS  error  model  from  ref.  7  was  updated  to  produce 
the  correct  line-of-sight  equations  for  the  VCOSS  model.  These 
equations  use  the  restrained  degrees  of  freedom  of  the  primary, 
tertiary,  and  focal  plane  surfaces  (Fig.  21)  together  with  the 
six  degrees  of  freedom  of  the  secondary  mirror  center  of  mass.  The 
CSDL  NASTRAN  NODE  100  DOF  #1  represents  the  LOSX,  DOF  #2  represents 
the  LOSY  and  DOF  #3  represents  the  DEFOCUS.  The  DEFOCUS  in  the 
NASTRAN  results  was  corrected  due  to  an  incorrect  model  coordinate 
for  node  40,  i.e.,  node  40  (x,y,z)  coordinates  of  (0,5. 0,2.)  should 
have  been  coordinates  of  (0,2. 5,2.).  This  error  was  on  the  LOS 
error  model  only  and  did  not  affect  the  NASTRAN  structural  model 
mode  shapes  (see  Table  13). 

i.  The  equipment  section  was  modeled  by  CSDL  in  NASTRAN  using 
rigid  body  equations  to  connect  node  44  to  nodes  42,  45,  47,  46,  43. 
This  section  was  modeled  by  LMSC  in  ASTRO  using  very  stiff  massless 
tubes  (outer  diameter  =  1M,  inner  diameter  =  0.9M)  to  connect  node 
44  to  nodes  42,  45,  47,  46,  43.  In  addition,  these  tubes  were  used 
around  the  perimeter  of  the  equipment  section  to  provide  more 
rigidity  by  connecting  nodes  41-45-47-46-43-42. 

j.  All  degrees  of  freedom  were  retained  and  used  in  the  modal 
analysis  by  CSDL.  Only  the  degrees  of  freedom  which  had  masses 
assigned  to  them  were  retained  and  used  in  the  modal  analysis  by 
LMSC.  The  massless  degrees  of  freedom  (DOF)  were  reduced  out  of 
the  stiffness  matrix  using  Guyan  reduction.  This  was  done  on  the 
massless  DOF  because  they  have  no  dynamic  effect  and  reducing 

them  produces  a  smaller  "STATICALLY  EQUIVALENT"  stiffness  matrix 
which  can  be  numerically  manipulated  with  much  greater  speed. 

k.  The  Mass  Matrix  for  the  VCOSS  model  at  coordinates  (0,0,0) 
is  shown  (Table  11)  for  both  LMSC  and  Draper  models  and  the  results 
are  in  excellent  agreement.  This  verifies  all  of  the  model  co¬ 
ordinates  and  mass  distribution. 
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Figure  21.  VCOSS  Optical  System 
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Table  12.  VCOSS  MODAL  FREQUENCY  COMPARISON  (Hz) 


MODE 


LMSC 

DRAPER 

-0.0534 

0.0000010 

-0.0462 

0.0000084 

-0.0290 

0.0000031 

-0.0269 

0.0000038 

-0.0215 

0.0000047 

-0.0154 

0.0000048 

0.1131 

0.1140 

0.1469 

0.1469 

0.1490 

0.1496 

0.1741 

0.1752 

0.4549 

0.4554 

0.5568 

0.5573 

0.5953 

0.5962 

0.6131 

0.6148 

0.6351 

0.6363 

0.6403 

0.6416 

0.8151 

0.8151 

0.8160 

0.8162 

0.8235 

0.8235 

0.9152 

0.9156 

0.9703 

0.9720 

1.158 

1.159 

1.551 

1.551 

1.773 

1.773 

2.250 

2.250 

2.254 

2.254 

3.434 

3.434 

3.452 

3.452 

3.957 

3.957 

3.986 

3.987 

4.051 

4.052 

4.337 

4.338 

6.550 

6.563 

8.057 

8.065 

8.434 

8.453 

8.840 

8.854 

8.991 

9.010 

MODE 

tt 

LMSC 

DRAPER 

38 

10.25 

10.34 

39 

10.51 

10.52 

40 

11.48 

11.55 

41 

11.65 

11.65 

42 

11.65 

11.65 

43 

11.65 

11.65 

44 

11.66 

11.66 

45 

11.74 

11.84 

46 

13.34 

13.50 

47 

14.10 

14.19 

48 

14.89 

15.03 

49 

16.51 

16.74 

50 

17.14 

17.32 

51 

17.51 

17.61 

52 

17.77 

17.77 

53 

17.77 

17.77 

54 

17.77 

17.77 

55 

17.78 

17.78 

56 

21.08 

21.32 

57 

21.66 

21.68 

58 

22.02 

22.16 

59 

23.47 

23.89 

60 

23.89 

23.94 

61 

24.35 

24.56 

62 

24.48 

24.66 

63 

25.46 

25.69 

64 

25.94 

26.11 

65 

26.44 

26.79 

66 

28.12 

28.41 

67 

28.72 

29.02 

68 

31.85 

32.38 

69 

32.55 

33.11 

70 

33.21 

33.56 

71 

35.03 

35.37 

72 

36.06 

36.60 

73 

38.76 

39.44 

74 

39.90 

40.17 

75 

40.04 

40.63 

57 


A  normal  modes  analysis  of  the  VCOSS  finite-element  model  was 
performed  using  the  ASTRO  finite  element  program.  The  resulting 
list  of  frequencies  up  to  40  Hz  have  been  tabulated  (Table  12). 

The  LMSC  model  was  created  with  the  ASTRO  computer  program  for  the 
following  reasons: 

a.  Was  used  as  a  verification  check  on  the  CSDL  model. 

This  procedure  enabled  CSDL  and  LMSC  to  eliminate  certain  discrep¬ 
ancies  and  inconsistencies  in  the  original  model  and  to  converge 
to  a  final  version  that  was  considered  as  a  baseline.  The  minor  re 
maining  numerical  differences  between  the  two  models  are  incon¬ 
sequential  as  far  as  the  control  studies  are  concerned. 

b.  LMSC  Model  interfaces  with  other  LMSC  control  analysis 
programs. 

c.  LMSC  Model  can  produce  16mm  movies  to  observe  modal  sine 
sweeps  and  LOS  jitter  on  ray  traces  from  mirrors. 

Comparisons  between  Draper  and  LMSC  models  show  very  good 
agreement  with  minor  differences  arising  from  different  modeling 
techniques  used  for  "rigid"  body  substructures  and  numerical 
round-off  caused  by  matrix  reduction. 

The  performance  of  this  system  is  measured  by  the  ability  to 
maintain  line-of-sight  (LOS)  and  defocus  within  specified  toler¬ 
ances.  The  results  of  the  LOS  for  sixty-nine  modes  (elastic  mode 
numbers  7-75  renumbered  as  mode  numbers  1-69)  are  tabluated  for 
each  mode  to  show  LOS-x,  LOS-y,  and  defocus  for  the  optical  axis 
(Table  14). 
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ACTUATOR /SENSOR  LOCATIONS  AND  TYPES 


Selection  of  actuator /sensor  typos  and  locations  is  a  necessary 
first  step  in  the  actual  process  of  control  design.  Before  this 
selection  can  be  made,  certain  ground  rules  must  be  established 
which  loosely  define  the  control  objectives  and,  thus,  partially 
motivate  the  hardware  selection.  For  the  VCOSS  system,  as  well  as 
for  many  others,  these  ground  rules  or  design  assumptions  play  a 
crucial  role  in  determining  system  performance,  viz.,  these  assump¬ 
tions  are  far  more  important  than  selection  of  control  design  algo¬ 
rithms  because  they  set  physical  constraints  on  achievable  perfor¬ 
mance.  For  this  study,  the  principal  guidelines  are:  (1)  vibra¬ 
tion  sources  cannot  be  isolated,  (2)  actuator /sensor  dynamics  must 
fall  in  the  bandwidth  of  the  evaluation  model,  (3)  control  design 
models  must  use  no  more  than  50#  of  the  evaluation  model  modes,  and 
(4)  more  than  one  type  of  control  system  will  be  required  (at  the 
system  level).  These  assumptions  very  substantially  impact  achiev¬ 
able  performance.  The  first  assumption  means  that  active  vibration 
isolation,  an  important  subject,  will  not  be  studied  now.  The  study 
objective  here  is  structural  control;  vibration  isolation  does  not 
raise  the  LSS  control  design  issues  (and  hardware  limitations)  like 
robustness,  spillover  insensitivity,  etc. ,  and  is  therefore  con¬ 
sidered  to  be  a  separate  design  exercise. 

The  second  assumption  means  that,  for  example,  very  high  band¬ 
width  actuators  may  not  be  either  feasible  or  desirable  and  the 
effects  of  such  limitations  will  be  investigated  and  not  ignored. 

The  third  assumption  is  critical  with  respect  to  rejection  of 
broadband  noise  in  a  structure.  Since  so  many  modes  are  excited, 
large  evaluation  models  are  necessary  to  evaluate  system  perfor¬ 
mance.  It  is  not  practical  (and  perhaps  not  even  possible)  to 
identify  large  portions  of  such  models.  Thus,  practical  control 
design  must  be  based  on  realistic  subsets;  this  sharply  limits 
closed-loop  performance  from  many  orders-of-magnitude  reduction 


factors  to  factors  on  the  order  of  100.  Very  high  performance 
against  broadband  noise  will  require  active  isolation  of  the 
disturbance  directly,  when  allowed  by  system  design  constraints. 

The  last  assumption  strongly  influences  interpretation  of 
results.  In  particular,  the  VCOSS  structure  exhibits  low-frequency 
bending  near  0.1  Hz.  Distrotion  of  LOS  performance  this  close  to 
DC  (zero  frequency)  will  be  controlled  by  mirror/optical  train 
alignment  systems  which  may  use  member /actuator  devices.  In  the 
analytical  work  below  only  one  type  of  actuator  dynamics  was 
studied  to  more  closely  conform  to  program  resources.  Since  con¬ 
trols  power  is  necessary  over  a  broad  range  of  frequencies,  mul¬ 
tiple  actuator  types  will  be  necessary  for  real  systems,  and 
several  different  actuator  models  will  habe  to  be  included  in  the 
desing.  These  complexities  are  left  to  future  work. 

These  guidelines  were  introduced  and  agreed  upon  during  the  first 
TD  meeting  review  at  AFWAL  on  June  16,  1982  and  reaffirmed  in  the 
second  TD  meeting  on  March  8,  1983.  It  is  important  to  understand 
the  reasons  for  these  guidelines.  For  the  CSDL  #2  problem  considered 
in  this  study,  isolation  of  the  vibration  sources  obviates  any  need 
to  control  structural  modes  because  no  other  disturbance  loads  are 
considered.  (This  is  not  the  case  in  a  real  system).  Structural 
control  solutions  which  use  all,  or  nearly  all,  of  the  evaluation 
model  as  a  basis  for  control  design,  can  produce  any  performance 
level  desired  (assuming  sufficiently  good  sensors  and  activators). 

It  is  unfortunate,  in  this  case,  that  the  broadband  noise  spectrum 
effectively  excites  the  entire  evaluation  model.  This  means  that 
high  performance,  near  the  specification  limit,  cannot  be  achieved 
without  isolation. 

Without  these  assumptions  on  isolation  and  synthesis  models, 
structural  control  performance  evaluations  are  meaningless  for  the 
VCOSS/CSDL  problem.  It  was  also  agreed  that  the  rigid  body  modes 
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would  be  eliminated  from  the  control  design  exercises.  This  means 
that  all  very  low  frequency  control  mechanizations  (rigid  body 
sensing,  surface  figure  sensing,  and  alignment)  are  eliminated 
from  consideration  (since  they  all  interact),  i.e.,  low  frequency 
control  cannot  be  studied  with  the  rigid  body  modes  deleted. 

For  the  baseline  VCOSS  control  system  implementation,  analysis 
of  the  line-of-sight  performance  indicates  that  linear  proof-mass 
actuators  installed  to  control  mirror  tilt  will  directly  influence 


the  desired  performance  metric  and  provide  adequate  modal  control. 


The  actuator  locations  are  shown  in  Fig.  20.  Each  driver  acts  m 
the  z-direction  which  gives  x  and  y  tilt  plus  defocus  control. 
Colocated  accelerometers  are  placed  at  each  proof-mass  actuator 
location.  Thus,  the  baseline  system  provides  9  colocated  actuators 
and  sensors  plus  a  focal  plane  processor/alignment  system  output 
which  provides  a  direct  measurement  of  LOS  error.  All  LAC  mechani¬ 
zations  require  the  colocated  sensor  outputs  while  high  authority 
control  (HAC)  and  adaptive  systems  require  the  focal  plane  sensor 
to  evaluate  the  performance  metric.  Note  that  the  actuators  are 
mounted  to  structural  nodes  and  not  directly  to  the  mirrors. 


In  phase  II,  some  of  these  actuators  are  removed  to  evaluate 
fault  tolerance  and  redundancy  requirements. 


3.  OPEN-LOOP  SYSTEM  EVALUATION  FOR  CONTROL  DESIGN 


a.  Linear  First-Order  Evaluation  Model 

A  linear  VCOSS  evaluation  model  is  formed  by  taking  the  first 
54  structural  modes  and  transforming  to  a  108th  order  linear  sys¬ 
tem.  Six  disturbance  inputs  are  used:  one  each  in  the  x,  y, 
and  z  directions  at  each  disturbance  location.  The  disturbance 
is  modeled  as  six  independent  colored  noise  sources  having  the  PSD 
given  in  Figure  2.  The  resulting  linear  evaluation  model  is: 
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This  gives  a  114th  order  linear  evaluation  model. 

Open-loop  frequencies  for  this  model  are  plotted  in  Figure  22 
(distregarding  the  noise  states).  Significant  frequency  "clustering" 
of  the  modes  is  evident.  The  effect  of  the  colored  noise,  as  seen 
through  the  line-of-sight ,  is  shown  in  Figures  23  through  27  which 
give  the  PSD  and  integrated  PSD  for  each  LOS  component.  The  LOSy 
level  is  an  order  of  magnitude  higher  than  that  of  LOSx.  The  domi¬ 
nant  modal  contributions  are  below  2  Hz  for  LOSy  and  below  1  Hz 
for  LOSx.  The  modal  properties  of  the  disturbance  are  evaluated  in 
detail  later. 

The  computation  of  these  curves  was  done  on  the  digital  com¬ 
puter  and  based  on  the  following  procedure.  The  closed-loop  system 
equations  are  represented  by: 

x  =  Fc  x  +  Dlu>l  +  D2u2  (22) 

LOSx  =  Hx  x 
LOSy  =  Hy  x 

where  x  is  the  state  vector,  u>l  and  u: 2  the  disturbance  inputs, 
and  Fc,  Dl,  D2,  Hx  and  Hy  constant  coefficient  matrices.  These 
equations  can  be  used  to  define  the  transfer  functions  between 
disturbances  and  LOS  errors.  For  instance 

,  LOSx/ul  =  Hx  (ju>I  -  Fc)  Dl  =  Txl  (u.')  (23) 
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the  PSD  of  LOSx  is  then  obtained  by: 

2  2  2  2 
PSDx  (u>)  =  |Txl(u>)  |  |w(u)l  +  |Tx2(u)  |  |w(u>)| 

2 

where  |w(u)l  represents  the  disturbance  PSD  as  defined  in 
Section  I,  Figure  2.  (wl  and  w2  have  some  PSD  and  are  uncorrelated) 
The  LMSC  analysis  software  computed  the  transfer  functions  Txl, 

Tx2,  Tyl,  Ty2  and  the  expression  which  was  plotted  as  a 
function  of  frequency  the  total  LOS  error  is  given  by 

L0S2  =  L0Sx2  +L0Sy2 


thus  the  total  PSD  was  obtained  by  simply  summing  PSDx  and  PSDy. 
Finally,  the  integrated  value  was  computed  by  the  integral 


PSD(u)  du 


(25) 


(Note  that  the  total  RMS  error  is  simply  yS(°°)) 


r 


,o.i 


0.  10.  20.  30.  40.  50. 

MODE  NUMBER 

Figure  22.  VCOSS  Evaluation  Model  Frequency  vs.  Mode  Number 
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Figure  25.  Integrated  PSD  of  Disturbance  to  Total  LOS 
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b.  Modal-Analysis  for  Mode  Selection-Theory 


The  size  of  the  model  and  uncertainty  in  the  high-frequency 
modes  makes  model-reduction  necessary  for  control  design. 

Although  a  114th  order  controller  poses  no  numerical  difficulties, 
high-frequency  feedback  based  upon  poorly  modeled  modes  is 
undesirable.  Several  criteria  for  mode  selection  are  used.  To 
account  for  model  uncertainty  at  high-frequency,  it  is  assumed 
that  the  model  is  unreliable  above  6  Hz.  In  other  words,  the 
first  25%  of  the  24  Hz  bandwidth  is  considered  accurate  enough  for 
control  design.  This  is  an  ad  hoc  limit. 

(1)  Internal  Balancing 

To  determine  the  most  important  modes  for  control  design, 
many  criteria  must  be  considered  including  controllability,  dis- 
turbability,  observability  in  performance,  and  observability  in 
the  measurements.  Any  mode  which  is  highly  controllable, 
observable,  and  disturbable  must  clearly  be  included  in  the  design 
model;  however  highly  controllable-but-unobservable  modes,  for 
example,  are  difficult  to  judge.  Moore  (Ref.  1)  has  devel¬ 
oped  an  "internal  balancing"  approach  whereby  asymptotically 
stable  linear  models  are  transformed  to  an  essentially  unique 
coordinate  representation  for  which  controllability  and  obser¬ 
vability  rankings  are  identical.  The  definition  of  internally 
balanced  coordinates  follows: 

•  De f ,  An  asymptotically  stable  model  j x  =  Ax  +  Bu  \ 

(y  =  cx  )  is 
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Notice  that  the  balanced  representation  is  such  that  the 
controllability  Gramian  and  observability  Gramian  are  equal  and 
diagonal.  The  ck's  are  termed  "second-order  inodes." 

In  general,  the  required  transformation  "scrambles"  the 
original  coordinate  system  such  that  the  physical  meaning  of  the 
states  is  lost. 
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However,  for  lightly  damped  structural  models  with  decoupled 
dynamics,  the  internally  balanced  coordinate  representation  is 
approximately  equal  to  a  scaled  representation  of  the  modal 
states.  Thus  it  is  possible  to  write  approximate  formulae  for  the 
{o2}  in  terms  of  the  original  model.  Three  modal  rankings  are 
considered : 
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•  disturbance  inputs  to  LOS 

•  actuator  inputs  to  LOS 

•  actuator  inputs  to  sensor  outputs 
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For  a  linear  first  order  differential  model  of  the  form: 
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The  approximate  second  order  modes  are  given  by  (see  Appendix  A) 
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(32) 


A  detailed  discussion  of  the  application  of  internal  balancing 
theory  to  lightly-damped  structural  models  is  provided  in  Appendix  A 


(2)  "Modal  Cost"  [  2  ] 


The  "modal  costs"  for  a  linear-quadratic-gaussian  system  are 
simply  the  diagonal  elements  of  the  cost  matrix  C  given  by 


Fx  +  Dw 

(33) 
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(34) 
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These  numbers  give  an  indication  of  the  relative  contributions  of 
each  mode  to  the  total  cost.  Although  mathematically  less  precise 
than  the  second-order  modes  (e.g.,  modal  costs  can  be  negative), 
modal  cost  analysis  gives  an  inexpensive  method  of  ranking  which 
can  also  be  used  for  closed-loop  analysis  where  the  low-damping 
property  no  longer  holds. 
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c.  Modal  Analysis  of  VCOSS  Mode 


Using  the  approximation  equations  (30)  through  (32)  the 
approximate  second-order  modes  are  computed  for  the  VCOSS  evalu¬ 
ation  model.  Open-loop  modal  costs  are  computed  using  the 
colored-noise  disturbance.  Here  the  absolute  values  of  the  modal 
costs  are  used.  The  RMS  second-order  modes  and  modal  costs  are 
plotted  versus  mode  number  in  Fig.  28.  Immediately  evident  is 
the  clustering  of  these  modal  phenomena.  The  disturbance  effect 
as  seen  through  the  line-of -sight  is  constrained  to  clusters  of 
modes  as  is  the  ability  to  measure  and  control  the  model.  The 
coincidence  of  the  controllable  clusters  and  disturbable  clusters 
indicates  a  favorable  actuator /sensor  configuration  for  the  prob¬ 
lem.  Table  15  gives  the  quantitative  modal  ranking  prescribed  by 
each  method. 
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Table  15  MODAL  RANKINGS  SUGGESTED  BY  THE  ANALYSIS 


DIST/LOS 

ACT/LOS 

ACT/ SEN 

MODAL  COST 

Mode 

°1 

Mode 

°1 

17 

0.0242 

17 

0.0366 

1 

0.0202 

1 

0.0130 

16 

0.0141 

7 

0.0119 

6 

0.0116 

34 

0.0081 

4 

0.0110 

27 

0.0080 

7 

0.0099 

3 

0.0078 

8 

0.0085 

16 

0.0074 

3 

0.0084 

39 

0.0072 

9 

0.0059 

41 

0.0066 

15 

0.0056 

40 

0.0051 

32 

0.0038 

6 

0.0042 

34 

0.0034 

28 

0.0042 

10 

0.0033 

50 

0.0041 

27 

0.0032 

18 

0.0039 

18 

0.0028 

8 

0.003  8 

31 

0.0025 

31 

0.0038 

39 

0.0024 

32 

0.0034 

28 

0.0024 

42 

0.0034 

Mode  o  j  Mode  ji^ 


5 

0.1836 

17 

0.0010 

17 

0.1813 

16 

0.0003 

7 

0.1037 

1 

0.0002 

2 

0.0941 

6 

0.0001 

41 

0.0831 

7 

0.0001 

34 

0.07  86 

8 

0.0001 

27 

0.0702 

4 

0.0001 

39 

0.0700 

32 

0.0001 

18 

0.0582 

15 

0.0000 

50 

0.0570 

34 

0.0000 

40 

0.0564 

9 

0.0000 

29 

0.0510 

3 

0.0000 

1 

0.0466 

27 

0.0000 

28 

0.0426 

39 

0.0000 

10 

0.0417 

31 

0.0000 

53 

0.0401 

28 

0.0000 

31 

0.0379 

42 

0.0000 

42 

0.0367 

41 

0.0000 

(1)  Actuator  Influence 

To  lend  further  insight  into  the  actuator/LOS  relationship, 
the  approximate  second-order  modes  were  computed  for  each  actuator 
with  respect  to  the  line-of-sight  and  are  plotted  in  Figs.  29,  30, 
and  31.  Actuators  2,  3,  5,  6,  8,  and  9,  have  very  similar  profiles 
to  each  other  and  to  the  total;  while  actuators  1,  4,  and  7  appear 
to  have  little  LOS  authority.  Thus  in  a  magnitude  sense  it  appears 
that  the  number  of  actuators  required  to  control  the  LOS  can  be 
reduced  to  as  few  as  3  (one  for  each  mirror),  although  the  con¬ 
sequent  loss  in  low  authority  control  (LAC)  authority  and  the  in¬ 
creased  force  requirement  might  make  this  reduction  desirable. 


Figure  29.  Controllability/LOS  Levels  for  Actuators  1,  2,  and 


Controllability/LOS  Levels  for  Actuators  4,  5,  and 


Figure  31.  Controllability/LOS  Levels  for  Actuators  7,  8,  and 


4.  CONTROL  DKGIGNS 

Three  controllers  are  designed  and  analyzed  for  the  VCOSS 
model.  The  first  uses  10  low-frequency  modes  from  the  model  and 
relatively  low  gains,  i.e.,  a  cautious  design.  Based  upon  closed 
loop  analysis  of  this  controller,  four  additional  modes  are 
selected  and  a  second  controller  is  designed.  The  bandwidth  of 
this  controller  is  allowed  to  expand  into  an  adjacent  "dead  zone" 
of  the  model  under  the  assumption  that  the  insensitive  modes 
cannot  cause  spillover.  A  low  authority  control  is  designed  and 
added  to  the  hi-gain  controller  to  give  the  third  design.  Full 
closed  loop  modal  and  stochastic  analyses  are  presented. 

a.  "Cautious"  Controller 

(1)  Mode  selection 

Using  the  modal  analyses  of  Fig.  28,  a  10-mode  design 
model  is  selected.  These  modes  are  shown  in  Table  15,  which 
gives  the  quantitative  modal  rankings  implied  by  each  method.  The 
first  10  "most  disturbable"  modes  with  respect  to  L£)S  were  used 
for  the  design.  Figure  32  shows  the  modal  selections. 

Inspection  of  the  second  column  shows  that  several  controllable 
modes  have  been  omitted  in  the  design  model  bandwidth.  It  is 
reasonable  to  expect  that  any  undesirable  effects  of  the  con¬ 
troller  will  be  related  to  these  modes. 

(2)  Control  design 
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The  10-mode  structural  design  model  transforms  into  a 
26th  order  linear  design  model  (including  the  disturbance). 

Using  standard  linear  quadratic  gaussian  (LQG)  techniques, 
a  state- feedback  control  is  found  which  minimizes 
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where  is  the  LOS  output  of  the  reduced  model.  Several  values 

of  b  were  tried  until  a  reasonable  closed-loop  bandwidth  was 
obtained  using  b  =  10-11.  Table  16  shows  the  resulting  closed- 
loop  design  spectrum.  The  significant  effect  of  this  control  is 
on  modes  17  and  8.  The  damping  on  mode  K  is  increased  by  nearly  a 
factor  of  100  hut  the  frequency  increases  by  less  than  one  per¬ 
cent.  Mod<>  17  experiences  a  damping  increase  of  a  factor  of  20 
and  the  frequency  increase  is  21%.  This  is  between  the  fre¬ 
quencies  of  modes  18  and  19  which  are  not  modeled.  Referring 
again  to  Fig.  28,  modes  19,  20,  and  21  are  apparently  uncon¬ 
trollable  with  respect  to  line  of  sight  as  well  as  sensor  place¬ 
ment.  Mode  18,  however,  has  some  controllability  (it  ranks  14th 
out  of  54).  The  control  bandwidth  is  therefore  likely  to 
influence  mode  18. 

Similarly,  a  standard  Kalman  filter  is  found  for  the  design 
model  using  the  actual  disturbance  PSD  and  LOS  "measurement  noise” 
of  intensity  5  x  10-1°.  This  choice  gives  a  bandwidth  which  is 
roughly  equivalent  to  the  control  bandwidth.  The  filter  spectrum 
is  shown  in  Table  17  as  is  the  control.  The  frequency  increase 
in  mode  17  is  sufficient  to  encomoass  unmodeled  mode  18. 

(3)  Evaluation 

The  10-mode  controller  was  implemented  with  the  full 
114th  order  VCOSS  model.  The  closed-loop  system  is  stable. 

Figures  33  through  36  show  the  closed-loop  and  open-loop  spectra 
with  varying  degrees  of  magnification  about  the  origin.  Four 
poles,  corresponding  to  modes  2,  5,  10,  and  18  move  to  the  right. 
The  spillover  in  mode  18  is  expected  since  its  frequency  lies 
within  the  controller  bandwidth  and  it  is  somewhat  controllable 
(see  Fig.  28).  Modes  2,  5,  and  10  are  unmodeled  modes  within 
the  open-loop  design  model  spectrum.  Although  their  modal  con¬ 
trollability  and  disturbability  is  small,  the  control  and  filter 
gain  are  sufficient  to  cause  an  "in-band"  spillover  at  these  fre¬ 
quencies  . 


K4 


*  **• 


Table  16  CLOSED-LOOP  CONTROL  SPECTRUM  OF  10-MODE  MODEL 


Seal 

part 

Iaag Inary 

part 

Frequency 

(Hz) 

Damping 

ratio 

-4.8001 

10.8240 

1 . 8845 

.4054 

-0.0270 

7.3328 

1.1671 

.0037 

-0.0127 

6.0956 

0.9701 

.0021 

-0.0115 

3.9841 

0.6341 

.0029 

-0.0114 

3.8621 

0.6147 

.0030 

-0.2178 

3.7528 

0.5983 

.0579 

-0.0114 

3.4864 

0.5549 

.0033 

-0.0033 

1.0962 

0.1745 

.0030 

-0.0222 

0.9388 

0.1495 

.0237 

-0.0203 

0.7341 

0.1169 

.0276 
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Table  17  CLOSED-LOOP  FILTER  SPECTRUM  OF  10-NODE  MODEL 


Real 

part 

Iaag inary 

part 

Fr eqnency 

(Ha) 

Daaping 

ratio 

-5.5453 

11.2239 

1.9925 

.4430 

-0.5170 

7.3093 

1.1662 

-0.0829 

6.0898 

0.9693 

.0136 

-0.0654 

3.9831 

0.6340 

.0164 

-0.2152 

3.9445 

0.6287 

.0545 

-0.4231 

3.6162 

0.5795 

.1162 

-0.0944 

3.6258 

0.5773 

.026  0 

-0.0877 

1.0982 

0.1753 

.07  96 

-0.0699 

0.9299 

0.1484 

.0750 

-0.1559 

0.7796 

0.1265 

.1961 

Figure  33.  //VCOSS  (  +  )  Open-Loop  Spectrum,  (X)  Closed-Loop 
Spectrum,  26th  Order 
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Spectrum,  26th  Order 


(X)  Closed-Loop 


The  stochastic  performance  with  this  controller  is  included 
In  T;il»  I  <•  '/JO.  Using  a  total  RMS  control  effort  of  12. 2N  the  RMS 
LOS  error  Is  reduced  to  16%  of  the  open-loop  value.  The  closed- 
loop  (absolute)  RMS  "modal  costs"  with  this  controller  are  shown 
in  Fig.  37.  Notice  that  the  effect  is  to  reduce  the  distur¬ 
bance  levels  significantly  in  modes  17,  16,  and  1;  which  are  the 
largest  open-loop  contributors.  Also  note  that  the  disturbance 
contribution  of  the  "spillover  modes"  2,  5,  10,  and  18  does  not 
increase. 

b.  High-Gain  Controller 

(1)  Mode  Selection 

The  10-mode  controller  caused  the  damping  in  the 
unmodeled  controllable  modes  2,  5,  and  10,  and  "high"  frequency 
mode  18  to  decrease.  To  account  for  this  effect  and  improve 
closed-loop  performance,  these  four  modes  were  added  to  the  design 
model;  resulting  in  a  34th  order  linear  control  design  model.  The 
14  modes  used  are  shown  in  Fig.  38.  The  in-band  unmodeled 
modes  11  through  14  have  very  low  controllability  and  are  not 
required.  From  Fig.  28  it  is  also  evident  that  modes  19 
through  26  are  highly  insensitive  to  control  and  disturbance 
inputs.  This  fact  is  used  to  justify  a  closed-loop  controller 
bandwidth  extending  into  this  "dead  zone." 

(2)  Control  Design 

As  before,  a  standard  linear  LOS  regulator  and  Kalman 
filter  is  designed  using  the  34th  order  design  model.  For  a  con¬ 
trol  penalty  of  5  x  10-12I  and  "measurement  noise"  of  5  x  10_11I, 
the  resulting  controller  has  bandwidth  of  3.5  Hz  which  is  well 
into  the  "dead  zone."  The  control  and  filter  design  spectra  are 
given  in  Table  18  and  19,  respectively. 
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Figure  37.  Closed-Loop  Modal  Costs  for  10-Mode  Controller 
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Table  18 

CLOSED -LOOP  CONTROL  SPECTRUM  FOR  14  MODE-MODEL,  b  =  S  x 


Real 

Imugi nary 

Frequency 

Dampi  n 

Part 

Part 

(Hz) 

Rat  io 

-13.7611 

16.8251 

3.4594 

0.6331 

-13.7611 

-16.8251 

3.4594 

0.6331 

-0.0310 

11.1483 

1.7743 

0.0028 

-0.0310 

-11.1483 

1.7743 

0.0028 

-0.0172 

7.3410 

1.1684 

0.0023 

-0.0172 

-7.3410 

1.1684 

0.0023 

-0.0129 

€.0954 

0.9701 

0.0021 

-0.0129 

-6.0954 

0.9701 

0.0021 

-0.9724 

3.8742 

0.6357 

0.2434 

-0.9724 

-3.8742 

0.6357 

0.2434 

-0.0178 

2.8514 

0.4538 

0.0062 

-0.0178 

-2.8514 

0.4538 

0.0062 

-0.0077 

3.4803 

0.5539 

0.0022 

-0.0077 

-3.4803 

0.5539 

0.0022 

-0.0095 

3.8664 

0.6154 

0.0025 

-0.0095 

-3.8664 

0.6154 

0.0025 

-0.0083 

4.0291 

0.6413 

0.0020 

-0.0083 

-4.0291 

0.6413 

0.0020 

-0.0097 

3.9788 

0.6329 

0.0024 

-0.0097 

-3.9768 

0.6329 

0.0024 

-0.0078 

0.7501 

0.1194 

0.0104 

-0.0078 

-0.7501 

0.1194 

0.0104 

-0.0027 

1.0983 

0.1748 

0.0025 

-0.0027 

-1.0983 

0.1748 

0.0025 

-0.0747 

0.9771 

0.1560 

0.0762 

-0.0747 

-0.9771 

0.1560 

0.0762 

-0.0023 

0.9227 

0.1469 

0.0025 

-0.0023 

-0.9227 

0.1469 

0.0025 
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Table  19 


CLOSED-LOOP  FILTER  SPECTRUM  FOR  14  MODE-MODEL.  v  =  5  x  10 


Real 

Part 

Imaginary 

Part 

Frequency 

(Hz) 

Damping 

Ratio 

-12.3240 

15.9457 

3.2075 

0.6115 

-0.0354 

11.1465 

1.7740 

0.0032 

-0.0354 

-11.1465 

1.7740 

0.0032 

-0.5400 

7.3256 

1.1691 

0.0735 

-0.5400 

-7.3256 

1.1691 

0.0735 

-0.0900 

6.0858 

0.9687 

0.0149 

-0.0906 

-6.0858 

0.9687 

0.0149 

-1.0419 

3.7793 

0.6239 

0.2658 

-1.0419 

-3.7793 

0.6239 

0.2658 

-0.3461 

3.9007 

0.6359 

0.0066 

-0.3461 

-3.9807 

0.6359 

0.0866 

-0.0222 

4.0329 

0.6419 

0.0055 

-0.0222 

-4.0329 

0.6419 

0.0055 

-0.0720 

3.9663 

0.6314 

0.0182 

-0.0720 

-3.9663 

0.6314 

0.0182 

-0.0771 

3.6160 

0.5756 

0.0213 

-0.0771 

-3.6160 

0.5756 

0.0213 

-0.0143 

2.8490 

0.4534 

0.0050 

-0.0143 

-2.8490 

0.4534 

0.0050 

-0.1039 

1.0994 

0.1758 

0.0941 

-0.1039 

-1.0994 

0.1758 

0.0941 

-0.1944 

0.7966 

0.1305 

0.2370 

-0.1944 

-0.7966 

0.1305 

0.2370 

-0.1713 

0.8836 

0.1432 

0.1903 

-0.1713 

-0.8836 

0.1432 

0.1903 

-0.0019 

0.9227 

0.1469 

0.0020 

-0.0019 

-0.9227 

0.1469 

0.0020 
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Using  the  34th  order  controller,  a  closed-loop  evalu¬ 
ation  with  the  full  VCOSS  model  was  performed.  Figure's  39 
through  41  show  the  resulting  closed-loop  spectrum.  No  "in- 
band"  spillover  occurs,  since  the  candidate  modes  were  included  in 
the  model.  Mode  27  experiences  a  significant  right-shift.  It  is 
the  closest  unmodeled  mode  to  the  controller  bandwidth  having 
significant  controllability  (see  Fig.  28).  Notice  that  as 
expected  no  problems  occur  in  the  "dead  zone"  modes.  Higher  fre¬ 
quency  spillover  effects  are  evident  in  modes  53,  52,  and  31. 

These  are  well  above  the  bandwidth  of  model  certainty  and  cannot 
be  explicitly  handled.  In  the  next  section  a  low-authority  con¬ 
trol  is  used  for  the  prevention  of  spillover  at  unknown  high  fre¬ 
quencies. 

The  performance  with  this  controller  is  included  in  Table  20. 
With  a  total  RMS  control  effort  of  29N  the  RMS  LOS  error  is 
improved  by  an  order  of  magnitude.  The  closed-loop  modal  costs 
with  this  controller  are  shown  in  Fig.  42.  The  modal  LOS 
disturbance  contributions  are  reduced  significantly  in  the  con¬ 
troller  bandwidth.  Note  that  the  modal  cost  for  mode  27  increses 
slightly  over  the  open-loop. 

c.  High-Gain  Controller  with  Low-Authority  Control 

The  purpose  of  low-authority  control  (LAC)  is  to  add 
damping  to  structural  modes  in  the  high-frequency  uncertainty 
region.  This  method  is  described  fully  in  reference  [3]  where  it 
was  used  successfully  to  eliminate  high-frequency  spillover.  For 
the  nine  colocated  actuators  and  sensors  of  the  VCOSS  model,  a  LAC 
is  designed  for  4  percent  damping  in  the  structural  modes. 


96 


Figure  40.  (+)  Open-Loop  Spectrum,  (X)  Closed-Loop  Spectrum 

34th  Order 


Figure  41.  (+)  Open-Loop  Spectrum,  (+)  Closed-Loop  Spectrum 

34th  Order 


0.A010 


Figure  42-  VCOSS  Open-Loop  (Open)  and  Closed-Loop  (Shaded) 
Modal  Costs  for  34th  Order  Controller 


(1)  Evaluation 


The  high-gain  controller  with  low-authority  rate  feed¬ 
back  was  evaluated  with  the  full  VCOSS  model.  The  closed-loop 
spectrum  is  shown  in  Figs.  43  through  45.  All  high-fre¬ 
quency  spillover  effects  appear  to  be  suppressed.  The  stochastic 
performance  is  shown  in  Table  20,  and  a  breakdown  of  actuator 
effort  is  given  in  Table  21.  There  is  a  20  percent  improvement 
in  the  RMS  LOS  error  over  the  high-gain  controller  used  alone, 
while  the  required  RMS  control  effort  is  reduced  to  25N.  Appar¬ 
ently,  the  LAC  causes  a  redistribution  of  required  control  effort 
among  the  actuators. 

One  possible  reason  for  this  may  be  found  in  the  sensor- 
actuator  placement.  Referring  again  to  Fig.  28,  modes  5  and 
17  are  very  significant  with  respect  to  sensor-actuator  effect. 

The  LQG  controls  may  be  working  very  hard  to  reduce  the  stochastic 
levels  in  these  modes  even  though  the  LAC  can  do  so  much  more 
cheaply. 
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Figure  43.  (+)  Open-Loop  Spectrum,  (X)  Closed-Loop  Spectrum 

HAC/LAC 


Figure  45.  (+)  Open-Loop  Spectrum,  (X)  Closed-Loop  Spect 
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Table  20 


STOCHASTIC  PERFORMANCE  ANALYSIS  QF  VCOSS  CONTROLLERS 


OPEN 

LOOP 


"CAUTIOUS" 

10-MODE 

CONTROL 


HIGH  GAIN 

14-MODE 

CONTROL 


RSS  LOS  x 
(y  RAD) 


131.7 


62.4 


69.0 


RSS  LOS  y 
( y  RAD) 


1109.5 


175.9 


87.5 


TOTAL  RSS 
LOS  (y  RAD) 


1117.3 


186.6 


111.4 


HIGH-GAIN 
CONTROL 
WITH  LAC 


26.4 


75.6 


80.1 


RSS  CONTROL 
EFFORT  (N) 


0 


12.2 


28.6 


25.3 


Table  21 

INDIVIDUAL  ACTUATOR  RMS  EFFORT  FOR  EACH  CONTROLLER  (N) 


Actuator 

10-Mode 

14-Mode 

14-Mode 

Number 

Control 1 er 

Controller 

Controller  with  LAC 

1 

6.1714*41 

1.32869+00 

2.0412 

2 

5.23329+00 

1.18129+01 

10.1947 

3 

5.08190+00 

1.18099+01 

10.1013 

4 

1.07119+00 

2.34089+00 

2.8637 

5 

4.92689+00 

1.18789+01 

10.4065 

6 

4.92749+00 

1.18239+01 

10.3886 

7 

1.98049-01 

5.64099-01 

1.7616 

8 

4.80799+00 

1.12139+01 

10.2338 

9 

4.72169+00 

1.11089+01 

10.0071 
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5.  IMPACT  OF  DISTURBANCE  MODEL  CHANGES 


The  preceding  analysis  presumes  six  equal  and  independent 
colored-noise  disturbances  —  as  modeled  by  equations  (18)  to 
(21).  It  should  be  noted  that  the  results  obtained  may  be 
scaled  to  give  a  solution  for  any  6th  order  white  noise  of  the 
form  (normal  distribution)  w  <—  N(0,wl).  This  scaling  does  not 
affect  the  relative  modal  rankings,  design  model  selection,  or 
controller  designs.  The  normalized  performance  analysis  for 
this  form  of  disturbance  is  shown  in  Table  22. 

If  for  computational  convenience,  the  number  of  disturbances 
is  reduced  to  two  (one  along  the  trisector  of  the  axes  at  each 
disturbance  location),  then  the  entire  solution  is  affected. 

The  modified  disturbance  distribution  will  change  the  filter 
gains  (as  well  as  the  number  of  noise  states  in  the  design 
model);  moreover,  the  altered  open-loop  modal  rankings  of  dis- 
turbance-to-LOS  and  disturbance-to-sensors  may  dictate  a  differ¬ 
ent  design  model  selection.  Thus,  in  general,  a  modified  dis¬ 
turbance  distribution  requires  a  modified  control  design.  These 
modifications  produce  minor  changes  in  controller  performance, 
i.e.,  this  is  not  a  robustness  issue. 


To  demonstrate  the  effect  of  disturbance  order  reduction 
upon  the  14-mode  controller  of  the  previous  section,  the  stoch¬ 
astic  performance  analysis  is  repeated  using  the  two-disturbance 
PSD  described  in  [5].  The  results  are  shown  in  Tables  23  and 
24.  Notice  that  the  total  RMS  LOS  reduction  is  slightly  better 
than  that  of  the  six  independent  disturbance  case:  a  factor  of 
17  for  the  former,  13  for  the  latter.  Since  closed-loop  perfor¬ 
mance  is  dominated  by  truncated-mode  response,  improvement  due 
to  controller  redesign  will  be  negligible  for  this  example  pro¬ 
vided  the  modal-fidelity  ground  rules  are  followed.  Therefore, 
no  new  design  is  generated.  Table  25  gives  the  normalized  per¬ 
formance  analysis  for  the  two-disturbance  case  where  w~N(0,wI). 
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;  Table  23 

1  STOCHASTIC  PERFORMANCE  OF  VCOSS  CONTROLLERS  -  2 

INDEPENDENT  DISTURBANCES 


■ 

HIGH  GAIN 

14-MODE 

CONTROL 

HIGH-GAIN 
CONTROL 
WITH  LAC 

RMS  LOS  x 
(p  RAD) 

77.72 

45.82 

16.19 

RMS  LOS  y 
(p  RAD) 

968.04 

55.33 

53.64 

TOTAL  RMS 

LOS  (p  RAD) 

971.16 

71.84 

56.03 

CONTROL 
EFFORT  (N) 


0 


17.12 


16.59 


ACTUATOR  NUMBER 

HAC 

HAC/LAC 

1 

1.0307 

1.2339 

2 

7.0623 

6.7832 

3 

7.3079 

6.9289 

4 

1.7503 

2.2331 

5 

7.0328 

6.7616 

6 

7 

6.9560 

0.4196 

6.6699 

1.3320 

8 

6.6080 

6.4624 

9 

6.6393 

6.4029 

TOTAL  RSS 

17.1223 

16.5913 
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Table  25 


NORMALIZED  STOCHASTIC  PERFORMANCE  OF  VCOSS  CONTROLLERS 
2  INDEPENDENT  DISTURBANCES 


H  =  -30  I  sec  1 

HIGH  GAIN 

HIGH  GAIN 

o 

OPEN 

14-MODE 

CONTROL 

W  =  wl  N  /sec 

LOOP 

CONTROL 

WITH  LAC 

RSS  LOS  x 
(y  RAD) 

0. 184  /m 

0. 109  /$r 

0.038  /f 

RSS  LOS  y 
( y  RAD) 

2.297  /iT 

0.132  /if 

0.127  /if 

TOTAL  RSS  LOS 
(Vi  RAD) 

2.304  /f~ 

0.170  /if 

0.133  /if 

RSS  CONTROL 

EFFORT  (N) 

_ 2 _ 

0.041  /¥ 

0.039  /if 

FcwSI 


’ -V-V 


6. 


LOW  FREQUENCY  ACTUATOR  LIMITATIONS 


a.  Force  Requirements 

Since  the  major  contributor  to  the  LOS  error  was  shown  to  be 
the  broadband  component  of  the  disturbance  PSD,  the  corresponding 
control  effort  is  also  mainly  due  to  this  component.  In  order  to 
evaluate  this  control  effort,  the  system  was  driven  with  the 
broadband  PSD  only  and  the  actuator  force  PSD  was  computed  for 
each  actuator.  The  total  RMS  forces  for  each  of  the  nine  actua¬ 
tors  are  shown  in  Table  26. 

One  may  already  notice  that  actuators  1,  4  and  7  require 
much  less  force  output  than  the  others.  This  results  from  the 
fact  that  the  controllability  of  the  LOS  with  respect  to  these 
three  actuators  was  much  less  than  the  other  six,  and  that  the 
optimal  controller  will  tend  to  use  more  the  most  efficient 
actuators.  The  evaluation  of  force  requirements  is  thus  narrowed 
by  only  considering  the  actuators  with  the  largest  RMS  in  each 
group  (i.e.,  primary,  secondary  and  tertiary).  Thus,  actuators 
3,  5  and  8  are  selected.  The  PSDs  corresponding  to  these  three 
actuators  are  respectively  given  in  Figs.  46,  47  and  48 
for  frequencies  ranging  from  0  to  10  rad/sec.  As  could  be  expected, 
the  maximum  control  effort  is  found  at  the  frequencies  of  the  con¬ 
trolled  vibration  modes.  The  force  requirements  may  now  be  broken 
down  for  each  frequency,  by  locally  integrating  the  area  under  the 
curve.  In  Fig.  49  for  instance,  the  scale  has  been  expanded 
to  give  a  better  display  of  actuator  tt 3  PSD  curve  around  0.75  rad/ 
sec.  Since  the  energy  is  very  localized,  one  may  separate  the 
requirements  at  different  frequencies.  Moreover,  since,  for  sizing 
purposes,  an  upper  bound  is  sought,  a  worst  case  level  of  effort 
table  can  be  built  by  determining,  at  each  peak  frequency,  what  was 
the  largest  effort  required  among  the  actuators  3,  5  and  8.  Table 
27  recapitulates  these  results.  The  bandwidth  of  the  curve  for 


Table  26 


CORRECTED  ACTUATOR  EFFORT  (N) 


LOCATION 

ACTUATOR  NUMBER 

HAC 

HAC/LAC 

1 

1.0307 

1.2339 

Primary 

2 

7.0623 

6.7832 

3 

7.3079 

6.9289 

, 

A 

1.7503 

2.2331 

Secondary 

5 

7.0328 

6.7616 

6 

6.9560 

6.6699 

7 

0.4196 

1.3320 

Tertiary 

8 

6.6080 

6.4624 

9 

6.6393 

6.4029 

TOTAL  RSS 

17.1223 

16.5913 
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ACTUATOR 


Figure  47.  VCOSS/PPM  Actuator  Power  Requirements 


Figure  48.  VCOSS/PPM  Actuator  Power  Requirements 
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Table  27 


which  the  KMS  value  was  computed  is  displayed,  with  the  corre¬ 
sponding  actuator  number.  For  instance  at  3.5  rad/sec,  actuator 
5  requires  the  most  force,  an  average  of  3.1  N  in  a  bandwidth  of 
0.033  rad/sec.  As  can  be  seen,  force  levels  are  relatively  small, 
around  or  below  5  N. 

b.  Proof-Mass  Actuator  Sizing 

In  proof-mass  actuators  in  general,  where  m  is  the  mass  of 
the  proof-mass  and  x  its  displacement,  the  force  delivered  is 
in  first  approximation  given  by: 

f  =  mx 

Thus,  if  f  occurs  about  the  frequency  to,  the  corresponding 
amplitude  of  the  displacement  is  given  by: 


For  a  given  force  level,  the  displacement  will  thus  depend 
upon  the  choice  of  m,  and  upon  the  frequency.  Since,  in  proof- 
mass  actuators,  the  displacement  is  limited,  the  mass  m  must 
become  large  when  the  frequency  is  low.  Table  28  shows  the 
displacements  required  at  the  various  frequencies  studied  earlier, 
for  three  values  of  the  mass.  As  expected,  the  strong  dependence 
upon  the  frequency  makes  the  0.75  rad/sec  case  the  worst  one. 

Even  with  a  30  kg  mass,  the  displacement  needed  is  18  cm. 
Although  such  mass  is  still  very  small  compared  to  the  total  mass 
of  the  spacecraft,  it  drives  the  design  of  the  actuator  toward 
larger  dimensions.  In  Fig.  50  are  shown  two  pivoted  proof-mass 
actuator  designs,  one  based  on  the  0.75  rad/sec  frequency,  one 
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Figure  50.  Actuator  Sizing 


based  on  the  next  frequency  (3.5  rad/sec),  the  so  called  "high" 
frequency  design.  Both  designs  utilize  a  50  N  electrodynamic 
motor,  but  different  masses,  length  and  angular  travel.  As  can 
be  seen,  the  "high"  frequency  design  is  more  compact  and  lighter 
and  will  be  easier  to  make  than  the  "low"  one.  However,  the 
latter  is  still  feasible,  with  a  lever  arm  below  1  meter.  A 
conceptual  arrangement  for  this  actuator  is  shown  in  Fig.  51. 

c.  Conclusions 

The  preceding  results  show  that  a  pivoted  proof-mass  (PPM) 
actuator  could  be  designed  to  cover  the  whole  range  of  frequencies 
required  to  control  the  LOS.  The  design  of  such  actuators  is  en¬ 
tirely  driven  by  the  lowest  frequency  to  be  controlled,  because  of 
o 

the  1/w  effect.  Considering  that  the  lowest  frequency  in  this 

case  was  0.75  rad/sec,  i.e.,  about  a  l/10th  of  a  Hertz,  it  appears 

that  it  falls  within  the  range  of  a  static  alighnment  system. 

Such  a  system  must  be  implemented  anyway  to  correct  for  structural 
deformations  due  to  launch  or  thermal  effects,  and  will  be  more 
appropriate  to  use  than  a  PPM.  The  alignment  system  will  generally 
use  position  actuators  to  move  the  mirror  surfaces  with  respect  to 

the  structure  and  will  have  enough  bandwidth  to  include  the  0.1  Hz 

mode.  The  higher  modes,  however,  can  be  taken  care  of  with  small 
size  PPM  actuators.  The  reader  should  review  the  initial  design 
guidelines  stated  earlier  to  place  these  conclusions  in  perspective. 
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7.  EFFECTS  OF  ACTUATOR  DYNAMICS 

Having  demonstrated  a  procedure  for  high-performance  control 
design  (HAC)  and  stability  augmentation  (LAC)  assuming  infinite 
actuator  bandwidth,  we  turn  now  to  the  problems  imposed  by  real¬ 
istic  actuators  upon  HAC/LAC  design,  stability,  and  performance. 
It  has  been  shown  [6]  that  the  transfer  functions  of  LAC  with  a 
pivoted-proof-mass  (PPM)  actuator  are  described  by 


f(s)/x(s)  = 


(1  +  t1s)(1  +  t2s) 


(c  >  o) 


(37) 


where  f(s)  is  the  force  on  the  structure,  c  is  the  active 

control  gain,  and  x(s)  is  the  structural  displacement  at  the 
actuator  location.  For  frequencies  smaller  than  (1/2tt)t“1  ,  the 
system  behaves  like  a  damper,  then  more  like  a  spring,  and  a 
final  roll-off  occurs  after  (\l 2^)?  ^  with  the  response  going  to 
zero  with  a  -1  slope,  as  shown  in  Fig.  52.  Since  the  "idealized" 

LAC  transfer  function  is 

f(s)/x(s)  =  cs, 

the  dynamic  effect  of  a  PPM  upon  a  commanded  force  f  is  given 

c 


f(s) 

fc(s) 


(1  +  T1S)(1  +  TgS) 


(38) 


Thus,  the  dynamics  impose  a  first-order  rolloff  on  the  con¬ 
trol  effort  at  frequency  (2ir  t^)-1,  and  above  (2rr  t2)_1 
exhibit  second-order  behavior. 

To  investigate  the  impact  of  these  actuators  upon  HAC/LAC 
controls,  the  following  assumptions  are  made: 
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Figure  52.  Frequency  Response  of  Active 
and  Passive  Dampers 

(1)  t1  and  x2  are  known  quantities.  They  are  readily 

determined  experimentally. 

(2)  (2tt  t^)  ^  is  much  larger  than  the  highest  known 

(accurately  modeled)  frequency  in  the  structural  model. 
Otherwise,  the  actuator  dynamics  can  (and  should)  be 
included  in  the  control  design  model  explicitly. 

The  degree  of  impact  upon  controlled  performance  depends  upon  the 
presence  of  unknown,  unmodeled  structural  modes  with  which  the 
actuator  might  interact.  An  analysis  using  a  class  of  "evaluation 
models"  containing  different  truncated  modes  would  show  the 
effects,  but  the  lack  of  a  common  reference  would  make  comparisons 
difficult.  Therefore,  a  different  approach  is  taken  here:  A 
single  "ideal"  controller  is  designed  and  evaluated  using  a 
representative  structural  model,  with  HAC  and  LAC  gains  chosen 
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to  reflect  a  "worst  case"  pathology;  namely,  destabilizing  HAC 
spillover  in  the  highest  mode  compensated  by  the  LAC.  Parameter 
studies  are  conducted  by  evaluating  stability  and  performance 
degradation  as  a  function  of  actuator  rolloff  characteristics 
using  this  baseline  design.  The  results  give  a  clear  picture 
of  actuator  effects;  although  in  realit  the  unknowns  are  high 
frequency  structural  modes  rather  than  actuator  dynamics. 

a.  Idealized  Baseline  Design  Example 

A  baseline  structural  HAC/LAC  control  design  is  produced 
using  a  34-mode  evaluation  model  and  two  colocated  actuator /sen¬ 
sor  pairs.  Referring  to  Fig.  28,  34  modes  give  a  credible 
representation  of  "known"  and  "unknown"  structural  modes  pro¬ 
vided  the  known  bandwidth  is  defined  as  0-2  Hz.  Actuator/ 
sensor  locations  2  and  3  (see  Fig.  20)  are  selected  based 
upon  the  performance  controllability  rankings  given  in  Figs. 

29  through  31.  These  plots  show  that  the  actuator/sensor 
complements  on  each  mirror  have  almost  identical  authority  to 
the  tota1  complement  (compare  with  ACT/LOS  in  Fig.  28),  and 
further  that  this  authority  is  dominated  by  the  actuators  on  the 
mirror  corners.  The  disturbance  model  used  is  the  two-disturbance 
PSD  described  in  [5]. 

Figure  53  shows  the  open-loop  modal  rankings  obtained 
using  this  model.  As  in  the  full  design,  mode  selections  for  a 
control  design  model  are  made  using  these  rankings;  however,  the 
range  of  model  certainty  is  presumed  to  be  0-2  Hz.  The  selected 
modes  are  shown  in  Fig.  54,  and  the  design  model  spectrum  is 
given  in  Table  29. 

Using  this  10-mode  design  model,  an  LQG  controller  is  com¬ 
puted  with  control  and  measurement  weights  selected  such  that 
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Figure  53.  Open-Loop  Modal  Analysis  of  34-Mode  Model 


Figure  54.  Open-Loop  Frequencies  of  34-Mode 
Evaluation  Model  and  Design  Model 
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Table  29 

FREQUENCIES  OF  10-M0DE  DESIGN  MODEL 


REAL 

I  MAG. 

FREQ. 

DAMPING 

PART 

PART 

(Hz) 

RATIO 

-0.0014 

0.7105 

0.1131 

0.0020 

-4.0019 

0.9359 

0.1490 

0.0020 

-0.0022 

1.0939 

0.1741 

0.0020 

-0.0070 

3.4984 

0.5568 

0.0020 

-0.0075 

3.7404 

0.5953 

0.0020 

-0.0077 

3.8519 

0.6131 

0.0020 

-0.0060 

3.9904 

0.6351 

0.0020 

-0.0146 

7.2774 

1.1582 

0.0020 

-0.0195 

9.7442 

1.5508 

0.0020 

-0.0223 

11.1387 

1.7728 

0.0020 
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the  closed-loop  design  spectrum  of  the  controller  goes  to  3  Hz. 
The  control  design  and  filter  design  spectra  are  given  in  Tables 
30  and  31,  respectively.  This  choice  of  gains  represents  an 
"optimistic"  high-gain  HAC  controller  since  the  closed-loop  band¬ 
width  exceeds  that  of  the  design  model;  and  indeed  the  resulting 
spillover  destabilizes  the  full  closed-loop  system.  A  LAC  of  the 
form 


u  = 


-1000 

0 


0 

-1000 


y 


(39) 


is  used  to  compensate  for  this  spillover. 

Thus,  the  desired  pathological  baseline  is  achieved.  The 
HAC  alone  causes  a  right-shift  in  the  unmodeled  modes,  actually 
destabilizing  mode  34,  and  the  LAC  is  required  for  system  sta¬ 
bility.  The  HAC  and  HAC/LAC  spectra  are  shown  in  Figs.  55 
and  56.  A  closed-loop  stochastic  analysis  is  given  in  Table 
32.  Note  that  the  total  vibration  reduction  is  only  a  factor 
of  7,  owing  to  the  limited  control  authority  of  the  two  actuator/ 
sensor  pairs.  This  is  further  reflected  in  the  component-cost 
comparison  of  Figs.  57  and  58,  which  show  a  higher-than-usual 
contribution  from  the  in-band  modes. 

b.  First-Order  Actuator  Effects 

The  first  parameter  study  is  performed  assuming  t2  =  0  in 
equation  (38),  representing  a  case  where  actuator/controller 
interaction  is  dominated  by  the  first-order  rolloff.  After 
appending  the  first-order  actuator  model  to  the  closed-loop 
evaluation  model,  the  closed-loop  HAC/LAC  spectrum  and  stochastic 
performance  are  evaluated  for  a  range  of  given  by: 

(2tt  Tj)'1  =  20,  18,  ...,  6,  4,  2,  1  (Hz)  (40) 
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Table  30 

CONTROL  DESIGN  SPECTRUM 


REAL 

I  MAG. 

PART 

PART 

-4.9711 

10.8875 

-4.9711 

-10.8675 

-0.0237 

11.1870 

-0.0237 

-11.1670 

-0.0282 

7.3414 

-0.0282 

-7.3414 

-0.3348 

3.7835 

-0.3348 

-3.7835 

-0.0117 

3.4855 

-0.0117 

-3.4855 

-0.0108 

3.9813 

-0.0108 

-3.9813 

-0.0087 

3.8578 

-0.0087 

-3.8578 

-0.0201 

0.7348 

-0.0201 

-0.7348 

-0.0297 

0.9422 

-0.0297 

-0.9422 

-0,0027 

1.0954 

-0.0027 

-1.0954 

-94.2478 

0.0000 

-94.2478 

0.0000 

FREQ. 

DAMPING 

(Hz) 

RATIO 

1.9020 

0.4160 

1.9020 

0.4160 

1*7773 

0.0021 

1.7773 

0.0021 

1.1884 

0.0038 

1.1684 

0,0038 

0.6013 

0.0886 

0.6013 

0.0886 

0.5547 

0.0034 

0.5547 

0.0034 

0.6336 

0.0027 

0.6336 

0.0027 

0.6140 

0.0023 

0.6140 

0.0G23 

0.1170 

0.0274 

0.1170 

0.0274 

0.1500 

0.0315 

0.1500 

0.0315 

0.1743 

0.0024 

0.1743 

0.0024 

15.0000 

1.0000 

15.0000 

1.0000 
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Table  31 


FILTER  DESIGN  SPECTRUM 


REAL 

PART 

-<U  Kfifi 

vTiLUvv 

-94.2478 

-5.4517 

-5.4517 

-0.0228 

-0.0228 

-0.5181 

-0.5181 

-0.2895 

-0.2895 

-0.1778 

-0.1776 

-0.0401 

-0.0401 

-0.0528 

-0.0528 

-0.1486 

-0.1488 

-0.0308 

-0.0308 

-0.0053 

-0.0053 


I  MAG. 
PART 

0.0000 

0.0000 

11.1664 

-11.1684 

11.1495 

-11.1495 

7.3096 

-7.3096 

3.5397 

-3.5397 

4.0071 

-4.0071 

3.9039 

-3.9039 

3.6795 

-3.6795 

0.7955 

-0.7955 

0.9365 

-0.9365 

1.0973 

-1.0973 


DAMPING 

RATIO 

1.0000 

1.0000 

0.4387 

0.4387 

0.0020 

0.0020 

0.0704 

0.0704 

0.0815 

0.0815 

0.0443 

0.0443 

0.0103 

0.0103 

0.0144 

0.0144 

0.1836 

0.1836 

0.0329 

0.0329 

0.0049 

0.0G49 


FREQ. 

(H.Z) 

15.0014 

15.0000 

1.9777 

1.9777 

1.7745 

1.7745 

1.1663 

1.1663 

0.5652 

0.5652 

0.6384 

0.6384 

0.6214 

0.6214 

0.5857 

0.5857 

0.1288 

0.1288 

0.1491 

0.1491 

0.1746 

0.1746 


233 


pop  Modal  Coordinate  "Component  Costs 
Corrected  Disturbance 


FLEXIBLE-MODE  NUMBER 


Open-Loop  Modal  Coordinate  "Component-Costs 
with  Corrected  Disturbance  —  Expanded  Scale 


Figure  59  shows  that  the  high-frequency  closed-loop  poles  move 
to  their  open-loop  values  as  rolloff  frequency  decreases,  which 
is  intuitive  since  the  actuator  dynamics  are  attenuating  the 
high-frequency  control  signals.  The  next  four  figures  zoom  in 
on  the  high-frequency  locus  to  show  how  these  poles  move. 

Figure  60  shows  the  first-order  actuator  effects  on  mode 
31.  Since  the  HAC  action  upon  mode  31  is  negligible,  this  figure 
shows  the  effect  of  the  first-order  rolloff  on  the  LAC  alone. 
There  is  no  destabilization  at  all;  the  locus  moves  uniformly  to 
the  open-loop  value  without  "overshoot."  It  will  be  shown 
analytically  in  the  next  section  that  the  first-order  rolloff 
cannot  destabilize  a  LAC. 

The  effects  on  modes  32  and  33  are  shown  in  Figs.  61  and 
62.  Here  LAC  is  compensating  for  pole-shifts  to  the  right  due 
to  HAC  spillover.  Notice  that  the  locus  crosses  0.2%  open-loop 
damping  line  before  reaching  the  open-loop  pole.  This  effect  is 
even  more  pronounced  for  mode  34,  shown  in  Fig.  63.  Since 
mode  34  is  unstable  by  design  under  HAC  alone,  and  stabilized  by 
LAC  to  less  than  the  open-loop  damping  value,  the  entire  locus 
lies  to  the  right  of  the  0.2%  damping  line.  Therefore,  insta¬ 
bilities  can  result  whenever  the  first-order  rolloff  frequency 
is  near  to  the  potential  spillover  region. 

It  is  clear  from  these  results  that  the  first  order  term 
"turns-off"  any  controller  action  beyond  its  bandwidth.  This 
applies  to  LAC  action  as  well  as  HAC  spillover  effects.  However, 
in  doing  so  the  actuator  can  permit  spillover  to  become  prob¬ 
lematic  where  LAC  action  is  critical.  There  are  two  ways  to 
alleviate  this  problem: 

(1)  Select  actuators  such  that  (2ir  is  sufficiently 

above  the  controller  bandwidth,  so  that  HAC  spillover 
cannot  be  present. 
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Figure  59.  Full  Closed-Loop  Spectrum:  X=HAC,  0=HAC/LAC, 
*=Locus  as  First-Order  Rolloff  Decreases 
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Figure  61.  Effect  of  First-Order  Actuator  Rolloff  Upon  Mode  32 


142 


(2)  Impose  a  separate  rolloff  on  the  HAC  control  signal 

to  suppress  spillover  effects  at  the  actuator  rolloff 
frequency. 

Either  way,  the  actuator  effect  will  not  be  destabilizing. 

The  effect  of  the  first-order  actuator  rolloff  on  closed- 
loop  stochastic  performance  is  shown  in  Fig.  64.  Clearly  the 
degradation  is  negligible  (less  than  7 1  at  6  Hz)  for  any  real¬ 
istically  high  actuator  bandwidth. 

c.  Second-Order  Actuator  Effects 


This  is  the  more  realistic  case,  where  Ti  ^  T2  >  ^ 

has  been  shown  [6]  that  the  LAC  root-shift  using  a  second-order 
actuator  is 


.rfe(dX  ) 
n 


1  "  TlT2U)n 


2  1  +  t 


k2)( 


i  ,  2  2 

1  +  T2“n 


2 

an 


(41) 


and  thus  that  in  the  absence  of  passive  damping,  a  LAC  will  always 
destabilize  structural  modes  whose  frequencies  are  beyond 
(2it/t~t^)"’  .  (Note  that  when  Tg  *  0  the  sign  of  d?e(dXn)  can¬ 
not  change;  a  condition  suggested  in  the  previous  section.)  For 
LAC  alone,  stability  can  be  guaranteed  if  sufficient  passive 
damping  is  present.  If  the  natural  structural  damping  is  too 
small,  additional  passive  damping  can  be  achieved  by  means  of  a 
passive  device  mounted  in  parallel  with  the  actuator.  Equation 
(41)  then  takes  the  following  form: 


f (s)/x(s) 


(42) 


where  p  is  the  passive  gain  of  the  damping  actuator,  and  the 
LAC  root  shift  equation  becomes: 


Figure  64.  Performance  Degradation  Due  to  First-Order  Actuator  Rolloff 
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>’  ^  2  Pa^L  »  P  I 


»?_  =  2  £  w 
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and  £on  is  the  natural  structural  damping.  A  necessary  and 


sufficient  condition  is  given  by  the  following  theorem: 


LAC  Stability  Theorem  (Aubrun):  Unconditional  stability  of  a  LAC 


system  is  guaranteed  if  and  only  if.  for  each  mode  n.  the  active 


damping  ratio  £cn  is  less  than  a  certain  maximum  C*n»  This 


maximum  active  damping  ratio  in  any  of  the  modes  within  the  band- 
width  of  an  active  LAC  controller  is  proportional  to  the  sum  of 


the  natural  structural  damping  Con  and  the  damping  £pn  intro¬ 


duced  by  a  passive  damper  mounted  in  parallel  with  the  actuator. 


Specifically 


**cn  ^  ^cn  ^max^on  +  ^pn^ 


(44) 


where  the  value  of  the  proportionality  constant  Rmax  is  given 


Rmax  =  ”’in  K  +  2  /K*  TlVTo 


(45) 


where 


K  =  ti/t2  +  t2/t1  +  2 


(46) 


and  where  1/t^,  l/Tg  are  the  P°les  of  the  active  damper,  and 


1/t  is  the  pole  of  the  passive  damper. 


The  proof,  along  with  a  detailed  exposition  on  LAC  design  and 
implementation,  may  be  found  in  [6], 


It  is  assumed  for  the  present  study  that  no  passive  damper 
is  present  ( Cpn  =  0),  and  that  =  ig.  This  is  the  "worst 


case"  of  second-order  actuator  since 


is  minimized. 


(1)  Second-Order  Effects  on  LAC  Alone 

To  fully  understand  the  impact  of  second-order  actuator 
dynamics  upon  HAC/LAC  controllers,  it  is  illuminating  to  look  at 
the  effects  upon  LAC  alone.  The  closed-loop  LAC-controlled  34- 
mode  spectrum  is  shown  in  Table  33.  For  the  "worst  case"  of 
t1  =  Tg  =  t,  closed-loop  spectra  are  computed  using  the  second- 
order  actuator  model  and  the  following  range  of  rolloff  fre¬ 
quencies: 


=20,  18,  . . . ,  6,  4,  2,  1 


(Hz) 


The  resulting  locus  is  shown  in  Fig.  65.  As  the  rolloff  freq- 
quency  decreases,  the  high-frequency  closed-loop  poles  move  to 
the  right  of  the  open-loop  0.2%  damping  line  before  reaching 
their  open-loop  values.  Figs.  66  through  68  show  this  motion 
more  clearly. 


There  are  two  points  to  be  made  from  these  figures.  Since 
the  active  term  in  the  root-shift  equation  (Table  33)  is  zero 
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Figure  65.  Effect  of  2nd-0rder  Actuator  Dynamics  Upon  LAC 
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Figure  67.  Effect  of  2nd-0rder  Actuator  Dynamics  Upon  Mode  33 
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when  ru>n  =  1,  the  closed-loop  damping  and  open-loop  damping  in  mode 
n  should  be  equal  when  the  rolloff  frequency  equals  the  modal 
frequency.  By  inspection  of  the  figures,  this  is  indeed  the 
case.  Secondly,  it  is  straightforward  to  show  that  equation 
(35)  achieves  its  minimum  when  Twn  =  /3*,  and  this  minimum  is 
evident  in  the  figures.  For  example,  the  11  Hz  mode  33  has  a 
predicted  minimum  when  the  rolloff  frequency  is  6.4  Hz;  and  this 
agrees  with  Fig.  67. 

(2)  Second-Order  Effects  Upon  HAC/LAC  Controller 

Using  the  "worst  case"  second-order  actuator  model,  the 
closed-loop  spectrum  and  stochastic  performance  of  the  baseline 
HAC/LAC  controller  are  computed  for  each  of  the  following  rolloff 
frequencies : 

(  2ttt  )-1  =  30,  28,  26,  ...,  6,  4,  2,  1  (Hz)  (43) 

The  locus  of  closed-loop  poles  is  shown  in  Fig.  69,  and  the 
high-frequency  loci  are  shown  in  expanded  scale  in  Figs.  70 
through  73. 

The  movement  of  mode  31  (Fig.  70)  is  essentially  the 
actuator/LAC  interaction  effect,  since  the  HAC  does  nothing  to 
this  mode.  In  fact,  the  0.2%  damping  crossing  and  the  minimum 
damping  occur  at  xun  =  1  and  xun  =  /3~  ,  respectively.  How¬ 
ever,  the  degree  of  right-shift  in  the  higher-frequency  modes 
is  much  more  pronounced  than  for  the  LAC  alone.  The  most  dramatic 
effect  is  upon  mode  34,  shown  in  Fig.  73.  Recall  that  the 
first-order  rolloff  caused  a  right-shift  in  this  mode  because 
the  destabilizing  HAC  overpowered  the  LAC  compensation  (see 
Fig.  63).  The  second-order  actuator/LAC  interaction  com¬ 
pounds  this  effect,  producing  a  more  severe  reduction  in  damping. 


Figure  69.  Effect  of  2nd-0rder  Rolloff  Upon  HAC/LAC: 

X=HAC,  0=HAC/LAC,  *=Locus  as  Rolloff 
Frequency  Decreases 
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Figure  70.  Effect  of  2nd-0rder  Rolloff  Upon  HAC/LAC:  MODE  31 
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Figure  72.  Effect  of  2nd-0rder  Rolloff  Upon  HAC/LAC:  Mode  33 


Thus,  the  LAC  stability  issue  exacerbates  the  loss  of  stability 
robustness  when  spillover  compensation  occurs  near  the  actuator 
rolloff  frequency. 

As  with  the  first-order  case,  this  problem  is  solved  if 
(2ht)-*  is  sufficiently  above  the  controller  bandwidth,  or  if 
the  HAC  is  rolled-off  artifically  well  below  the  actuator  fre¬ 
quency.  In  either  case,  however,  LAC  destabilization  is  possible. 
Sufficient  passive  damping  should  be  added  to  the  actuator  such 
that  if  an  unmodeled  mode  were  present  at  the  frequency  where 
the  root-shift  (equation  41)  is  minimum,  the  closed-loop  damping 
would  be  acceptable. 

The  impact  of  second-order  actuator  dynamics  upon  HAC /LAC 
performance  for  this  example  is  shown  in  Fig.  74.  Although  the 
degradation  is  slightly  higher  than  for  the  first  order  case,  it 
is  still  minimal  for  any  reasonable  actuator  rolloff. 


SECOND— ORDER  ROLLOFF  FREQUENCY  (HZ) 


Figure  74.  Performance  Degradation  Due  to  Second-Orde 
Actuator  Rolloff 
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SENSOR  MODELING  EFFECTS 


The  principal  terms  of  sensor  modeling  errors  will  be  bias, 
scale  factor,  and  misalignment  errors  and  these  errors  will  create 
LSS  control  mechanization  difficulties  at  low  frequency.  Since 
rigid-body  controls  are  being  put  aside  for  the  present  study,  only 
the  high-frequency  effects  for  accelerometer  mechanizations  will  be 
examined  here.  Optical  sensing,  used  extensively  in  large  space 
structures  (LSS)  control,  generally  exhibits  high  signal-to-noise 
performance.  These  systems  are  prone,  however,  to  measurement 
matrix  errors.  This  is  a  more  general  type  of  error  which  is 
discussed  (in  broad  context)  in  subsection  9. 

For  high-performance  accelerometers  used  above  0.5  Hz  or  so, 
broadband  noise  is  the  principal  hardware  characteristic  which  may 
be  important  for  control  performance  evaluation.  Typical  noise 
power  spectra  for  such  instruments  are  shown  in  Fig.  75.  The 
equivalent  rate  noise  is  also  estimated  in  the  figure.  Integrated 
noise  figures  of  30  ng  (la)  are  currently  attainable  at  modest 
cost  for  units  like  the  Bell  model  7.  For  distributed  sensing, 
the  accelerometer  noise  figure  (at  high  frequency)  does  not  cause 
significant  performance  degradation  as  shown  in  Table  34.  Here 
the  equivalent  white  rate  noise  is  propagated  through  the  system 
using  the  34  mode  evaluation  and  the  two  actuator/four  sensor 
"reduced  set"  described  in  subsection  7.  To  see  an  effect,  however 
the  white  rate  intensity  was  increased  by  about  9  orders-of-magni- 
tude  as  shown.  Thus,  it  can  be  concluded  that  accelerometer  usage 
for  LSS  control  will  not  be  limited  by  broadband  accelerometer 


noise. 


TOTAL  RSS 
LOS  (v  RAD) 


970 . 72 


132 . 93 


13A.58 
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EFFECTS  OF  GENERAL  ACTUATOR /SENSOR  MODEL  INC,  ERRORS 


a.  Robustness  of  Control  Design 

Since  the  spacecraft  model  is  a  baseline  representation  of 
the  actual  system,  the  control  design  must  be  able  to  accommodate 
errors  in  the  model.  Typical  sources  of  model  error  in  spacecraft 
systems  include: 

•numerical  errors  due  to  approximate  modeling  tech¬ 
niques,  e.g.,  high  order  NASTRAN  models. 

•actual  parameter  changes  in  the  LSS,  e.g.,  thermal 
effects,  gravity,  spacecraft  and  antenna  dimensions, 
mass  distributions,  etc. 

•unmodeled  dynamics,  e.g.,  effect  of  reduced  order 
modeling,  neglected  residual  modes  ("spillover"), 
neglected  actuator /sensor  dynamics,  non-linearities,, 
etc . 

•incomplete  data  obtained  from  on-earth  testing,  e.g., 
partially  assembled  structures  in  simulated  zero-g. 

Therefore,  uncertainty  in  the  baseline  model  arises  from  both 
actual  causes  and  intentional  approximations  of  complicated  phe¬ 
nomena.  In  many  cases,  these  are  indistinguishable. 

b.  Modeling  and  Model  Uses 

In  a  broad  sense,  model  error  refers  to  measurable  differ¬ 
ences  between  two  systems.  Normally,  one  of  these  is  the  actual 
system,  and  the  other  is  a  "model"  of  the  actual  system.  However, 
since  models  are  developed  for  different  uses,  there  can  be  sev¬ 
eral  models  of  an  actual  system;  consequently,  model  error  can 
also  refer  to  measurable  differences  between  two  models  of  the 
same  system.  For  example,  the  model  used  in  the  evaluation  phase 
of  a  control  design  is  a  high  fidelity  representation  of  the 
spacecraft;  this  model  is  usually  referred  to  as  the  baseline 


model.  The  model  used  for  control  design,  however,  is  consider¬ 
ably  simpler;  it  is  usually  obtained  from  the  baseline  model  by 
any  number  of  methods,  e.g.,  modal  truncation,  linearization,  etc. 


Thus,  there  are  several  models  for  any  system--dif f erent 
models  are  used  for  control  design,  analysis,  and  evaluation.  In 
fact,  there  is  no  unique  model  of  a  system,  because  models  are 
determined  by  their  use.  But  model  use,  in  turn,  is  determined 
by  performance  requirements.  In  other  words,  the  limitations  or 
errors  associated  with  a  particular  model  are  intimately  con¬ 
nected  with  the  performance  demands  of  control  system  design  based 
on  the  model. 

c.  Model  Error 

To  obtain  a  quantitative  bound  on  model  error,  the  source  of 
the  error  must  be  determined,  and  then  a  suitable  test  procedure 
devised.  Consider  a  controlled  spacecraft,  as  depicted  in  Fig.  76, 
with  the  following  model: 

Sensor  Model 

z  =  y  +  n 

Dynamic  Model 

y  =  M(s)u  +  d 

Controller 

u  =  C(s)e,  e  =  r  -  z  .  r  =  reference 

where  M(s)  is  a  finite  dimensional  transfer  function  matrix 
model  representing  the  dynamics  of  the  actuators  and  spacecraft; 
C(s)  is  the  transfer  function  matrix  of  the  controller;  d  is 


gure  76.  Block  Diagram  of  Controlled  Spacecraft 


the  effect  of  external  disturbance  forces  on  the  sensed  outputs; 
and  n  the  sensor  noise.  Let  the  actual  dynamics  be 

y  =  P(s)u  +  d  (44) 

where  P(s)  is  not  necessarily  finite  dimensional.  For  example, 
P(s)  can  be  either  a  high  order  NASTRAN  model  or  represent  data 
from  the  actual  system,  whereas  M(s)  is  the  reduced  order  con¬ 
trol  design  model.  Thus,  the  controlled  output  is: 

y  =  (I  +  PC)-1  d  +  (I  +  PC)-1  PC  (r-n) 


“  Hyd  d 


+  Hyr  (r-n) 


(45) 


(The  complex  variable  's'  has  been  suppressed  for  brevity  of 
notation,  unless  needed  for  clarification.)  Note  that 


+  H 


yr 


I 


Since  M  is  a  reduced  order  model  of  P,  let 


(46) 


P  =  M  +  A„ 
r 


(47) 


where  Ar  represents  the  effect  of  neglected  residual  modes. 


Thus: 

H  .  *  H  .  -  (I  +  MC)  X[I  +  A  C(I  +  MC)  ]  A  C( I  +  MC) 
yd  yd  1  r  r 

H  ,  -  ( I  +  MC)-1  (48) 

yd 

where  H  , 


is  the  nominal  transfer  function  with  no  model  error 


can  be  obtained  from  the  simple 


A  bound  on  the  size  of  A 

r 

input /output  test: 


6r<w) 


RMS  (y-Mu) 
RMS  (u) 


(49) 


u„  sin  ait 
o 


where  RMS(x)  is  the  root-mean-square  (RMS)  value  of  x.  Thus, 
injecting  sinusoidal  signals  into  the  actual  system  and  the  model 
and  comparing  the  outputs  gives  a  bound  on  Ar«  In  fact,  it  can 
be  shown  that 


a  [Ar(jw)]  <  6r(W)  (50) 

where  a(A)  is  the  maximum  singular  value  of  the  complex  matrix 
A,  i.e.,  the  square  root  of  the  maximum  eigenvalue  of  A*A. 


Similarly,  we  can  examine  the  way  in  which  other  kinds  of 
model  error  enter  into  the  closed-loop  dynamics.  The  spacecraft 
model,  for  example,  may  not  include  actuator  dynamics.  This 
omission  can  be  represented  by  the  model  error  form, 


P  =  M(I  +  A  )  (51) 

a 

where  A  represents  the  deviation  from  the  dynamics  of  an  actu- 

£L 

ator  with  infinite  bandwidth.  In  this  case,  the  model  error  test 
is: 

(52) 

u  =  u  sin  cot 
o 

The  effect  on  closed-loop  response  is: 

Hyd  =  Ryd  “  <I+MC)_1  M[ I+A&( I+CM)_1CM] -1AaC( I+MC )_1  (53) 


q[Aa(  ju>)] 


6a(u))  = 


RMS  (y-Mu) 
RMS  ( Mu ) 


d.  Stability  Robustness 


The  model  error  bounds  can  be  used  to  quantitatively  deter¬ 
mine  the  robustness  of  stability,  i.e.,  stability  margins.  If  the 
nominal  spacecraft  model  M  is  stabilized  by  the  controller  C, 
then  the  actual  spacecraft  is  stable  if: 

(2)  For  reduced  order  model  errors  &  , 

r 

6r(w)a{C( jw)[I+M( ju))C( ju>)]-1}  <1,  w  >  0  (54) 

or  (1)  For  actuator  model  errors  A. 

a 

6a(a))a{[I+C(ja))M(ja)>]"1C(ja))M(ja))}  <1,  w  >  0  (55) 

Table  35  summarizes  these  robustness  tests  for  generic 
model  errors  bounded  by 

o[A(.jaj)]  <  6  (  u> ) ,  (56) 

where  <5(w)  is  determined  from  input/output  tests,  e.g.,  RMS 
tests.  Table  35  shows  the  stability  margins,  denoted  by  6sm. 
defined  as  the  maximum  bound  on  model  error,  which  (at  the  speci¬ 
fied  location,  e.g.,  actuator,  sensor,  etc.)  insures  stability. 

Thus , 


6(w)  <  <$sm(w>  ,  u  >  0.  (57) 

e.  Performance  Robustness 

Once  having  determined  frequency  dependent  bounds  on  the  model 
a  similar  procedure  can  be  used  to  bound  the  effect  of  model  error 
on  performance,  i.e.,  performance  robustness.  For  example,  define 
the  frequency  dependent  performance  measure  p(to)  as, 


(58) 


p ( u> ) :  =  o[Hyd(.jco)  -  Hyd(.ja>)]/afHyd(.ja))] 

This  gives  a  relative  measure  of  the  deviations  of  Hyd  about 

the  nominal  H  In  the  face  of  reduced  order  model  error  A  a 

yd  r 

specified  level  of  performance  robustness  p(w)  is  achievable  if 
Sr(w)a{C(  ju))}o{  [  I+M(  jto)C(  ju>)  ]-1} 

- =—  <  p(w),  u  >  0  (59; 

1  -  6r(a))c{C(ju))tI+M(jai)C(jw)]“x} 

Similar  expressions  for  performance  robustness  can  be  obtained 

as  a  result  of  other  sources  of  model  error,  e.g. ,  sensor  model 

error.  Table  35  summarizes  these  performance  robustness  tests 

for  generic  model  error.  The  table  shows  the  performance  margins, 

denoted  by  6  ,  defined  as  the  maximum  bound  on  model  error  (at 

pm 

the  specified  location)  which  insures  the  specified  performance 
tolerance.  Thus, 

6(w)  <  Spm(w)  ,  w  >  0  (60) 

guarantees  performance  robustness.  This  also  guarantees  stability, 
since, 


6__(<d)  <  6  (o> )  ,  0)  >  0.  (61) 

pm  sm 

f.  Example:  Actuator  Uncertainty 

In  this  section  we  present  an  example  which  illustrates  the 
use  of  the  robustness  analysis  techniques  just  described.  The 
problem  considered  is  to  determine  the  effect  of  actuator  band¬ 
width  on  performance  and  stability  of  the  baseline  control  design 
for  the  CSDL-II  (rev.  3)  model.  The  control  system  uses  2  actu¬ 
ators  and  4  sensors  located  as  shown  in  Fig.  20.  The  control 
design  (see  subsection  7)  is  based  on  a  model  with  no  actuator 


dynamics,  i.e.,  the  actuators  are  assumed  to  have  infinite  band¬ 
width.  For  the  purposes  of  analysis,  both  of  the  actuators  are 
modeled  identically  by, 


where  u  is  the  force  out  of  the  actuator;  u  is  the  force  com- 

3. 

mand  (control  signal);  and  H  (s)  is  the  transfer  function  of 

ci 

each  actuator  given  by 


Ha<s>  *  ' T~~  - —5 

s  +  2«a“as  +  “a 


u  =  2-nf 
a  a 


(63) 


where  f  is  the  bandwidth  and  £  is  the  damping  ratio,  which 

3.  3 

is  assumed  to  be  fixed  at  £  =1.  The  bandwidth  is  then  varied 

a 

until  instability  ensues.  From  subsection  7c  it  is  determined 
that  performance  is  degraded  by  no  more  than  10%  for  f  >12  hz, 
and  stability  is  guaranteed  for  f  >  4  hz.  We  will  now  compare 
these  exact  results  with  the  approximate  results  obtained  using 
the  formulae  in  Table  35.  The  relevant  stability  and  performance 
margins  (from  the  Table)  are, 


=  l/a[(I+CM)-1CM] 


(64) 


^pm  (Ssmp{p+l5sm^C^I+Mc^  1]a[(I+MC)  1M]/a[(I+MC)  1]| 

The  model  error  is  taken  to  be 

Aa(s):=  I  -  Ha(s)I  (65) 

Thus,  A&(s)  represents  the  error  between  the  ideal  infinite  band¬ 
width  actuators  and  the  actual  ones.  As  a  result. 


Figure  77  shows  the  effect  of  actuator  bandwidth  on  the  stability 
margins,  i.e.,  plots  of  5gm  vs.  w  and  6a  vs.  w  with  actuator 
bandwidth  f  =  3.7  Hz.  This  corresponds  to  the  lowest  actuator 

SL 

bandwidth  that  ensures  stability,  i.e.,  6  (f  )  =  6  (f  ).  Com- 

pare  this  to  3  Hz  which  was  obtained  from  an  eigenvalue  analysis 
with  the  actuator  dynamics  included.  The  small  discrepancy  is  due 
to  the  inherent  conservatism  in  the  robustness  analysis. 

Similarly,  Fig.  78  shows  the  effect  of  actuator  bandwidth 
on  the  performance  margin,  i.e.,  plots  of  6pm  vs.  u>  and  Sa  vs.  w 
with  actuator  bandwidth  f  =  160  Hz.  This  corresponds  to  the 

SL 

lowest  actuator  bandwidth  which  ensures  a  10%  performance  degrada¬ 
tion  (p  ■  0.1)  where  ( f  Compare  this  to  20  Hz 

obtained  from  the  full  system  eigen-analysis.  The  discrepancy 
here  is  large,  principally  due  to  the  difference  in  performance 
measures.  Further  basic  research  is  required  to  alleviate  this 
conservatism. 


10. 


RESULTS  COMPARISON  FOR  STIFFNESS  DESIGNS 


One  of  the  VCOSS  program  objectives  is  to  compare  performance 
results  for  a  "soft"  actively  controlled  structure  with  those  of  a 
passive  "stiff"  structure.  Such  comparisons  are  difficult  to  inter¬ 
pret,  however,  if  each  structure  is  not  optimized  to  enhance  either 
a  passive  or  active  controls  approach.  Notwithstanding  these  diffi¬ 
culties,  the  comparisons  can  now  be  examined.  For  the  passive  struc¬ 
tural  performance,  the  results  are  given  in  Table  36,  due  to  CSDL, 
for  the  combined  broadband  and  discrete  power  spectra.  Three  models 
are  shown  for  both  the  six  independent  disturbance  sources  and  the 
two  source  "equivalent"  input  (cf.  subsection  5).  Some  disagreement 
still  remains  over  the  open-loop  response  for  two  sources,  so  dis¬ 
cussion  here  will  be  restricted  to  the  six  source  results.  Rev  00 
is  the  original  ACOSS  model  while  Rev  03,  the  VCOSS  model,  is  a 
softened  version  of  Rev  00.  Open-loop  response  degrades  by  about 
a  factor  of  3  for  the  "soft"  model.  A  stiffer  version,  Rev  04, 
produces  about  a  4  to  1  performance  spread  between  Rev  04  and 
Rev  03.  Thus,  it  is  reasonable  to  expect  these  kinds  of  improve¬ 
ments  through  structural  stiffening  exercises. 

The  closed-loop  comparison  is  recapitulated  in  Table  37. 

The  open-loop  values  have  shifted  slightly  since  the  discrete  power 
spectra  are  not  included.  The  reduction  factor  in  closed-loop  is 
now  about  14,  using  6  effective  actuators.  More  actuators  can  boost 
this  value  to  about  80.  The  principal  observation  is  that  much 
larger  reduction  is  achievable  through  simple  active  controls  than 
via  simple  structural  stiffening.  Performance  sensitivity  to  model 
uncertainty  will  also  be  better  for  the  closed-loop  systems.  In 
either  case  here,  sufficient  reduction  factors  to  meet  the  speci¬ 
fication  are  not  obtainable  without  resorting  to  integrated  controls 
including  isolation  and  low-frequency  alignment  systems--all  con¬ 
sidered  out-of-scope  for  this  study.  Nevertheless,  the  simple 
actuators  used  in  this  study,  coupled  with  low-frequency  alignment 


to  handle  0.1  Hz  modes,  can  produce  significant  performance  improve¬ 
ments  without  major  impacts  on  structural  redesign  and  without  com¬ 
plex  modeling  changes  to  synthesize  the  controls. 


Table  36.  PASSIVE  SYSTEM  RESPONSE  TO  COMPLETE  PSD 


2  INPUTS 


LOSX 

LOSY 

LOS 

REV00 

161.38 

602.71 

623.94 

REVO  3 

161.61 

1949.25 

1955.94 

REVO  4 

73.77 

485.71 

491. 28 

6  INPUTS 


LOS  DEF 


LOSX  LOSY 


91.88  364.11  375.52 

135.88  1116.37  1124.61 

83.82  264.90  277.84  6.69 


Table  37.  REV  03  CLOSED-LOOP  PERFORMANCE  WITHOUT  DISCRETE 
FREQUENCY  POWER  SPECTRA 

LOSX  LOSY  LOS 

Rev  03  131.7  1109.5  1117.3  Open-loop 

26.4  75.6  80.1  Closed-loop 
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11. 


CONCLUSIONS 


This  first  VCOSS  study  focused  on  the  system  performance 
limitations  imposed  by  actuator,  sensor,  and  signal  processing 
hardware.  To  provide  some  depth  to  the  study,  restrictive 
assumptions  were  made  at  the  outset  to  bound  the  investigation. 

The  most  significant  of  these  assumptions  is  the  removal  of  the 
rigid-body  control  problem.  It  is  in  this  region  that  sensor 
modeling  is  most  important;  attitude,  alignment,  and  modal  control 
systems'  authority  also  must  be  resolved  at  very  low  frequency. 

These  complex  issues  have  not  as  yet  been  adequately  addressed 
and,  thus,  results  shown  here  must  be  properly  interpreted.  A 
system  design  methodology  which  includes  hardware  constraints  has 
been  shown  to  be  effective  in  producing  implementable  control 
system  designs.  The  process  may  be  highly  iterative  to  meet  sys¬ 
tem  design  constraints.  A  preliminary  hardware  selection  for  the 
VCOSS  studies  includes:  (1)  proof-mass  dampers,  (2)  accelerometers, 
and  (3)  a  candidate  digital  processor  architecture.  The  main  con¬ 
clusions  of  VCOSS  to  date  deal  with  high-frequency  actuator  dynam¬ 
ics,  and  this  report  makes  a  significant  contribution  in  this 
area. 


Realistic  actuation  adds  a  significant  twist  to  the  HAC/LAC 
stability  robustness  issue  due  to  the  combined  effects  of  spill¬ 
over/compensation  rolloff  and  LAC/actuator  interaction;  however, 
this  can  be  handled  in  a  straightforward  manner.  It  is  shown, 
in  fact,  that  any  desired  high-frequency  stability  margin  can  be 
guaranteed  via  HAC  frequency  shaping  and  passive  damping  actuator 
augmentation  for  any  realistic  (PPM)  actuator.  The  impact  of 
realistic  actuators  upon  performance  robustness  is  shown  to  be 
minimal,  requiring  no  design  modification  at  all. 

The  effects  of  actuator/sensor  modeling  errors  and  the  juxta¬ 
position  of  these  errors  with  spacecraft  modeling  uncertainty 


appears  to  bo  a  solvable  problem  using  frequency-domain  error 
criteria  1 i ke  singular  value  deeomposit ion  of  the  system  return 
difference.  Combinations  of  appropriate  frequency  shaping,  hard¬ 
ware  selection,  and  realistic  assessment  of  model  error  can  allow 
LSS  control  hardware  realizations  to  produce  predictable  perfor¬ 
mance.  It  remains  to  test  the  validity  of  these  statements  through 
appropriate  ground  and  flight  experiments. 
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A  new  and  coapnt a t i ona  1  ly  simple  method  for  flexible 
spacecraft  model  redaction  is  presented.  Using  an  internal 
balancing  approximation  shown  to  be  valid  when  natural  damping  is 
small,  the  method  provides  quantitative  modal  rankings  with 
respect  to  disturbance  environment,  actuator  authority,  sensor 
observability,  and  performance  objective.  These  rankings  are  used 
to  select  a  reduced  set  of  structural  modes  for  controller  design, 
and  also  to  anticipate  potential  closed-loop  performance  and 
stability  problems  resulting  from  modal  truncation.  The  method  is 
demonstrated  using  a  54-mode  spacecraft  example. 
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Tbe  deployment  of  large  apace  structure*  will  require  feed¬ 


back  control  to  meet  stringent  vibration,  1 i ne  -  o  f  -  s  i  g b t ,  jitter, 
and  surface  quality  performance  criteria.  Flexible  spacecraft 
control  design  requites  mathematical  models  of  plant  dynamics, 
a c t u a t o r / a e n s o r  locations,  and  disturbance  environment,  all 
defined  with  respect  to  performance  objectives.  These  "evaluation 
models",  obtained  by  f ini t e -e 1 ement  analysis,  can  contain  hundreds 
of  structural  modes  known  with  varying  degrees  of  accuracy.  This 
paper  addresses  open-loop  reduction  of  these  models  for  the  pur¬ 
pose  of  feedback  control  design. 

Model  reduction  for  control  design  is  motivated  by  practical 
issues:  hardware  limitations,  reduced  computation  power  on-line, 
and  robustness.  In  principle,  an  "optimal"  control  can  be  gener¬ 
ated  for  a  large  model  using  modern  control  techniques,  but  this 
control  is  at  least  as  complex  as  the  model  itself.  When  the 
"reduced  controller-order"  constraint  is  formally  imposed,  optimal 
control  methods  lose  their  attractive  closed-form  solutions. 
Therefore,  a  reduced-order  "control  design"  model  is  selected 
using  the  same  criteria  as  for  the  evaluation  model:  namely,  to 
represent  the  system  accurately  with  respect  to  the  performance 
obj  actives . 

The  typical  approach  to  structural  model  reduction  is  mode 
selection.  The  simplest,  domi nan t- f r e quency  selection1,  ignores 
the  fact  that  actuator  and  sensor  placement,  disturbances,  and 
performance  requirements  also  affect  modal  dominance.  Modal-cost 


analysis  '  addresses  the  latter  two  issues  by  prescribing  nodal 

rankings  based  upon  relative  contributions  to  a  stochastic  cost 

2 

functional.  It  is  shown  that  the  nodal-costs  decouple  asymptoti- 
cally  as  danping  approaches  zero.  The  influence  of 
actuator/ sensor  placenent  upon  nodal  doninance  is  not  treated, 
however,  and  the  problen  renains  that  feedback  control  can  couple 
the  disturbance  into  open-loop  "undisturbed"  nodes. 

Recently  a  gene r a 1/kode 1-re  doc t ion  approach  based  upon  state 
selection  fron  a  grannian-bal  anced  ("internally-balanced")  coor¬ 
dinate  representation  has  been  propoied^'*’*.  For  a  given 
input/output  configuration,  this  approach  produces  a  "balanced 
app  rox  ina  t  i  on"  to  the  large  nodel  by  defining  and  retaining  the 
"most  controllable/obaervable"  states.  Open-loop  application  of 
this  reduction  method  to  pro blens  involving  two  sets  of  inputs 
(disturbances  and  actuators)  and  outputs  (regulated  variables  and 
sensors)  is  in  general  not  possible,  because  a  different  coor¬ 
dinate  basis  is  required  for  each  input/output  pair.  However, 

results  for  1 i ght ly-dampe d  single-input  single-output  flexible 
7 

systems  have  shown  that  modal  truncation  and  balanced  approxima¬ 
tion  are  equivalent,  asymptotically,  as  damping  approaches  zero. 
The  implied  asymptotic  equivalence  of  modal  and  balanced  coor¬ 
dinates  extends  to  the  multi-input  multi-output  case,  and  is 
fundamental  to  the  present  work. 

This  paper  gives  a  systematic  m ode- se 1 e c t i on  procedure  for 
large  lightly-damped  structural  models,  based  upon  the  theory  of 

it  4 

internal  balancing"  .  The  mode  1 - r e due t i on / con t r o 1-d e s ign  problem 
is  described  in  Section  II.  Section  III  briefly  reviews  internal 


A- 4 


balancing  theory  ,  develops  an  a pp r o x in  a t i on  for  multi-input 
multi-output  1 i gh t 1 y - d am pe d  structures,  and  investigates  the 
approximation  accuracy  for  practical  problems.  Similarities  with 

g 

recent,  independently  obtained  single-input  s in g 1 e- ou tpu t  results 
are  noted.  Section  IV  attempts  a  physical  interpretation  of  the 
results  with  respect  to  the  control  design  issues,  and  proposes  a 
new  mo d e - s e 1 e c t i on  methodology  for  lightly-damped  structures.  A 
complete  design  example  it  given  in  Section  V. 


XI.  Problea  Statement 


Consider  an  evalnation  nod  el  containing  structural  dynanics 
in  aodal  fora: 


ij  +  2ZQi)  +  On  =  ^u  +  W* 


x  *.  y/jti  +.y/jH  (i 


snt  s ) 


(1) 


y  *  +  Y\ 2 n  (regulated  variables) 


J  *  lia  E(y  y)  (perforaance  aeasure) 
t-*» 


where  ntRn,  Z  *  diagtt^},  0  a  diag{u>^),  u  is  the  control  vector, 
and  w  is  the  disturbance  vector.  Defining 


bj  *  i  *  th  row  of 


d  *  i '  th  row  of  C/ 


a  j  j  =  i  *  th  col  nan  of 


Cj4  -  i'th  coluan  of 


an  equivalent  state-space  representation  of  (1)  is: 


Mode-selection  for  control  design  requires  evaluation  of  the 
relative  iaiportance  of  each  node  to  the  control  problem.  The 
specific  issue's  which  influence  this  evaluation  are: 

•  modal  fidelity 
•controllability 
• observabil ity 
•disturbance  environment 

•  performance  objective 


None  of  these  criteria  stands  alone.  For  example,  it  makes  little 
sense  to  retain  the  most  accurately  modeled  mode  if  it  is  neither 
controllable  nor  observable  in  the  mea snrement s .  Similarly,  a 
highly  controllable  mode  adds  little  to  the  controller  performance 
if  it  is  (strictly)  unobservable  in  the  sensors  and  in  the 
performance.  The  problem  is  to  find  a  reduced-order  model  which 
addresses  each  of  these  design  issues  such  that  controls  designed 
for  the  reduced  model  perform  well  on  the  evaluation  model. 
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III.  Internal  Balancing  of  L If h t ly-Dampe d  Strnctirti 


In  general,  balanced  model  reduction  is  state  selection  in  a 
special  coordinate  system.  Before  proceeding  with  mod e - s  e  1  e  c t i on , 
this  theory  is  reviewed. 


Balanced  Modal  Reduction 


Consider  a  linear,  t im e - i nv a r i a nt  ,  asymptotically  stable 
state-space  system: 


=  Ax  +  Bu 


having  controllability  grammian  W  and  observability  grammian 

c 

W2  given  by 


r  Atn„T  At. 

W  =  J  e  BB  e  dt 

c  a 


W2  ■=  /  e A  tCTCeAtdt 

0  . 


2  2  T  T 

AW  +  W  A  +  BB  =  0  (4) 

c  c 


T  2  2  T 

AW  +  W  A  +  C  C  =  0  (5) 

o  o 


The  model  (3)  is  internally  balsnced  if 


2  2  „  2 

W  =  W  ■=  l 
c  o 


where  I 


•  •  • 


diet  Iffj.  o2* 


'  #a 


and  i < j  =>  o  ^  >  o  ^  . 


2. 


The  o . *i  are  teracd  the  second-order  aodet  of  (3).  It  has  been 


shown  that  any  aiodel  of  this  fora  (3)  can  be  taken  to  internally 
balanced  form  using  similarity  transformations. 

Next  consider  the  balanced  model  partitioned  as: 


2  J 


y  =  [  C2  C2  ]  x 


with  «=  diag  {  }  and  L  *  *=  {  o*  +  1 


,o  } .  The  essence 
& 


4  6  2  2 

of  balanced  model  reduction  '  is  that  if  0^  >>  then  the 


input  affects  x2  much  less  than  it  affects  x^  ,  and  the  output  is 


affected  by  x2  much  less  than  by  x^. 


The  internally  balanced  coordinate  representation  has  a 

number  of  desirable  properties  with  respect  to  model  reduction. 

It  is  unique  (to  within  a  sign  change  on  the  basis  vectors) 

2  2 

provided  the  0 ^ * s  are  distinct.  The  e^'s  *re  similarity  in¬ 
variants  of  (A.B.C).  Most  remarkably,  any  arbitrary  subsystem  is 


guaranteed  to  be  asynptotically  stable  (subject  to  an  additional 
restriction  on  tbe  basis  vectors  in  the  nondistinct  case*’). 

Application  of  balanced  reduction  to  very  large  models  is  no 
easy  task.  Indeed,  "transformation  methods"  are  generally  un¬ 
desirable  for  large  structural  models  because  of  computational 
problems  and  loss  of  a  physically  meaningful  state  vector.  In  the 
following  it  is  shown  that  the  balanced  coordinates  are  a  special 
case  of  modal  coordinates,  and  thus  that  no  explicit  change  of 
basis  is  necessary,  provided  the  damping  is  very  small  and  the 
frequencies  are  sufficiently  distinct. 

Application  to  Light ly-Danped  Strnctnral  Models 

Consider  the  state-space  representation  of  a  structural 
dynamics  model  in  modal  form: 


Theorem :  For  any  modal  subsystem  of  (6),  the  transformation 
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A  T  ,,  T  ,  r 

Ti  "  “iC2iCli/(c2iC2i)-Ci 


produces  the  following  internally  balanced  aodal  ivbiystea: 


qi  =  “i 


— 2  C  v 

S  li 


-2t<v2; 


1/2 


fb  bT 
'  i  i 


ay  b 

li  l  i  i 


L'°2iV2ibi  J 


A- 12 


£roon  The  proof  follows  from  direct  substitution  of  into  (4) 
end  ( 5 )  . 


2-2 

Noting  tbit  Oj  .  ~  %  i  w^en  <  1  »  define  the  bl  ock-d  is  gone  1 

transformation: 


0  °i/“i 


Vi*i 


where 


■l  4  [  bibi  1  ] 


and  V1#  Vj,  and  y ^  are  as  defined  above.  Applying  (7)  to  (6),  the 
following  model  is  obtained: 


$i  “  Ai<li  *  Bj®  • 


1  V* 
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L  <i+2^i 


V  ‘  V  ) 

li  lx 
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Main  Remit 

(i)  If  <<  1.  i  =  l....,n,  then  each  model  subsystem  of  (9) 
is  approximately  internally  balanced  with 


-  s, 


A  1 


2 

a  . 
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(4Vl) 


-1 


[  Vi 


,  T 
(c.  ,c 


li  li 


*  -2  T  ,  1 


1/2 


(10) 


(ii)  If  ma  x  (  ^  i  ,  Cj  )  ■*  z  (  ••  ^  >  m  j  )  /  I  u  j  I  <<  1  »  i/j,  then  the 
entire  model  (9)  is  approximately  internally  balanced  with 


A1 


a2 1 
n 


(11) 


EXfifif 


.T, -S  ^  ■ 


(i)  The  error  in  the  balancing  approximation  for  each  modal 


tabaystem  ia: 


I2  -  Z 
Li  Li 


2  2 
"li'0! 


a  2  2 

0  °2i"°i 


1  h-[ 


,1/2 

2  T.  ^  T  J 


2  T  .  T 

(I)  c  C  +  C  c 

i  li  li  2i  2i 


i  1  1  '  [ 


2  . [Iy4 l+( 1+Y2) 1/2]  .1/ 2 

- 1 - 1 - 1 -  1  I 

1  +  dYj  +  C.)2  J 


<  1/2 

-  1-C1-3C.)  '  provided  X,  £  <  1/3 


(ii)  Consider  any  two-mode  subsystem  of  the  transformed 
model.  The  controllability  grammian  for  this  subsystem  is  ap¬ 
proximated  by: 


I  P 
*  j  i 

pji 


■  here  satisfies 


AjPji  +  PjiAi  +  BjBi  =  0 
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•  >  .  »  .  •  .  « 


Solving  for  and  taking  the  L^-nora  leads  to 


II 


Pji  11  <  2oj0i(1/Pi+1/P2>[^j“J^i“i/C<l-^j>(l-l;i))  ] 


1/2 


P1'P2  =  [«j+“>i  +  2Cj«.>j  Ci»i±2»js»i((l-cJ)  (1-e*))172] 


1/2 


Reaoving  the  subscripts  to  denote  the  saiinos ,  ve  obtain 


II  Pj.  Il/o2  i  2C(a»/Aw  +  l)/(l-C)  -  2(^w/Aw+C)  if  £<<  1 


from  which  it  follows  that 


II  'j.  II 

-  <<  1  provided  Um/Ao.)  <<  1 


Similar  results  can  be  obtained  for  the  observability  grammian. 


Rate  Output  Case 

Vhen  c2.=0,  i=l,...,n,  the  approximation  (10)  is  exact  and 
each  modal  subsystem  of  (6)  is  balanced  regardless  of  its  damping 
ratio.  To  illustrate  the  approximation  error  in  (11),  consider 
the  two-mode  subsystem  obtained  from  (6)  assuming  ^ e  £  and 
normalizing  the  first  frequency  to  unity: 
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I 

I 

I 


0 


M  L  •  •  s  * 


“2  C  (  1  +  Aw/ a) )  -(  HAw/w) 


JM  in 


1  *11  *21  I  C 1 2  *22  1 


(12) 


•here  q  =  dq/dx  and  x  =  tu.  Choosing  b ^  -  b  ^  =  c^  ~  1 
the  true  (I)  and  ipproxiaate  (I)  balanced  grammians  are  computed 
for  several  values  of  £  and  a  range  of  Au.  The  plot  of  ap¬ 
proximation  error  in  Fig.  1  shows  the  improvement  in  modal 
decoupling  as  damping  ratio  decreases.  At  t,  -  0.001,  for  ex¬ 
ample.  a  1%  frequency  difference  results  in  only  5%  approximation 
error  . 


lisnlacement  Output  Case 


When  c^  =  0,  i=l,...,n,  the  modal  balancing  error  bound  is 


I2  -  I 
Li 


<  |  l-[l-2Ci(l  +  C2)1/2-2!;2]1/2  |  <  ^(2  +  2^  +  C2 


When  X,  ^  <  .1  this  bound  is  comparable  to  than  that  given  in  Ref.  8 
for  the  SISO  case: 


.r 


Figure  1.  App r o z  is ■ t i on  Error  for  Two-Mode  Subsystem 
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4^(  1  +  2 t j)1/2(1+2^1)'1 (2  +  2;^)_1/2 


Choosing  b 


C22  in  the  displacement  output 


balancing  error  is  compnted  and  plotted  in  Fig.  1. 
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IV.  Mod*  Selection  Methodology 


The  ipproxiBite  nodal  balancing  result  enables  quantitative 
nodal  analysis  of  the  evaluation  nodel  with  respect  to  physically 
neaningful  design  issues.  Four  nodal  rankings  are  discussed, 
followed  by  a  suggested  node- se 1 e c t i on  procedure. 


Distnrbene*  to  Kegnlated  Variables: 

2 

The  t°nr  )  give  a  nodal  ranking  in  terns  of  open-loop 
L  i 

perfornance.  Large  values  indicate  high  disturbance  propagation 
to  the  output,  while  snail  values  inply  low  open-loop  perfornance 
contribution.  If  the  goal  were  sinply  to  match  open-loop  perfor¬ 
mance,  node  selection  would  be  based  on  these  rankings.  Observe 
2 

that  the  lo..  }  are  equal  to  the  balanced  steady-state 
Di  i 

covariance  when  the  nodel  is  forced  with  un i t - i n t e n s i t y  white 

2 

noise.  An  interesting  relationship  exists  between  the  {e^g  ^ 

i 

3 

and  the  modal  costs  of  the  nodel  (2).  Assuming  a  zero-mean 
un i t - int e n s i t y  white  noise  disturbance  and  light  damping,  the 
open-loop  modal  costs  are  approximated  by 


(4tl^~1["i"T<*'2icli,li*c2lc2i)] 


Using  (10),  the  following  relationship  is  obtained: 


i 


4? 


,  w 


4 

iffi 


Thus,  modal-cost  and  1 n t e r na 1 -b a  1 anc 1 ng  criteria  can  differ. 


For 


example,  if  two  modes  hive 
favors  tbe  higher  frequency  while  balancing  ranks  them  eqoally. 


cr°j 


and  then  modal-cost 


Aotuators  to  Regulated  Variables: 


I 


2 

BC 


These  rankings  show  modal  controllability  of  the  performance 

2 


A  small  valne  of 


BC 


indicates  that  the  given  actuator  con¬ 


figuration  has  little  direct  effect  upon  the  contribution  of  mode 

i  to  the  performance,  regardless  of  its  open-loop  performance 

2  2 

contribution.  In  particular*  if  crfiC  /II^BC^  8ttal*  and 

i 

2  2 

/  I  I IpC  I  I  is  large,  a  redesign  of  the  actuator  configuration 
i 

is  suggested. 


2 

Disturbances  to  Sensors: 


2 

The  ^  show  modal  observability  of  the  disturbance  in 

i 

2 

the  sensors.  A  mode  having  a  small  relative  may  be  impos- 

i 

2 

sible  to  estimate  on-line.  If  the  corresponding  Op^,  is  large, 

i 

the  selection  of  sensor  locations  or  types  is  inappropriate. 


2 

Aotuatora  to  Sensors: 


The  provide  a  modal  analysis  of  potential  controller 

authority  for  the  given  the  actuator  and  sensor  configuration. 


Ideally  the  nodes  with  large  ogM  should  align  with  those  having 

i 

2 

large  0^^  so  that  controller  authority  matches  the  performance 
i 

2 

objective.  A  Mode  with  large  onu  should  be  included  in  the 

B  M  . 

1 

2 

design  model  even  if  0^  is  low.  particularly  if  it  is  in  the 

i 

controller  bandwidth,  to  prevent  spillover  problems. 


Node-Seleotion  Procedure 


Using  the  four  modal  rankings,  the  node-selection  process  is 
performed  as  follows: 

2 

1.  Select  the  modes  having  the  largest  0^  .  These  modes 

i 

contribute  most  to  the  performance  objectives,  and  assum¬ 
ing  reasonable  actuator  and  sensor  placements,  it  should 

be  possible  to  control  each  of  the  modes  to  some  extent. 

2 

2.  Examine  the  {0  }.  Include  in  the  design  model  any 

D  M  . 


1 

highly  con t r o 1 1 ab 1 e /me  a su r abl e  modes  not  selected  in  (1), 

especially  if  they  are  close  in  frequency  to  selected 

modes.  Omission  of  these  modes  can  cause  spillover, 

which  can  destabilize  the  system. 

2  2 

3.  Examine  the  { 0^  }  and  the  ^CTgc  ^  '  Unselected  modes 

1  i 

having  large  values  in  either  of  these  rankings  indicate 

2 

actuator/ sensor  configuration  pathologies.  A  large  0 

DM  . 


1 

indicates  an  unmodeled  mode  in  the  measurements,  which 
will  inhibit  state  estimation.  An  unmodeled  mode  with 
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V.  A  Dtil|a  Example 


This  application  is  performed  on  the  Charles  Stark  Draper 

9 

Laboratory  model  til  (revision  til)  which  has  54  flexible  modes, 

each  with  0.2%  open-loop  damping.  Six  disturbance  inputs  are 

modeled  as  independent  colored-noise  sources,  each  having  a  15  Hz. 

2 

rolloff  and  a  mean-squared  value  of  600n  N  .  There  are  nine  colo¬ 
cated  f o r c e - a c t u a t o r / r a t e - s e n s o r  pairs  and  an  additional 
measurement  of  1 i n e - o f - s i g h t  (LOS)  pointing  error.  Open-loop 
modal  frequencies  are  given  in  Table  1.  The  linear  evaluation 
model  has  114  states  (108  structural  and  6  disturbance). 

This  example  follows  the  two-level  "  h i g h- a u t h o r i t y  /  1  o  w  - 
authority"  control  design  approach  described  in  Ref.  10.  The 
balanced  mode-selection  methodology  is  applied  for  control  of  LOS 
pointing  performance  (HAC),  and  " 1 ow- an thor i  ty  control  (LAC)** 
using  colocated  actuators  and  rate-sensors  is  appended  to  the 
control  law  to  prevent  high-frequency  "spillover"*^. 

Modal  Analysis 

Using  the  approximate  internal  balancing  relation  (10),  the 
following  modal  rankings  are  obtained: 


•  Disturbances 


LOS 


Table  1  Opea-loop  evaluation  aodel  frifitaoiit 


Mode  Frequency  Mode  Frequency  Mode  Frequency 

No.  (Hz)  No.  (Hz)  No.  (Hz) 


1 

0.1131 

2 

0.1469 

3 

0.1490 

4 

0.17  41 

5 

0.4549 

6 

0.5568 

7 

0.5953 

8 

0.6131 

9 

0.6351 

10 

0.6403 

11 

0.8151 

12 

0.8160 

13 

0.8235 

14 

0.9152 

15 

0.9703 

16 

1.1582 

17 

1.5508 

18 

1 .7728 

19 

2.2499 

20 

2.2541 

21 

3.4335 

22 

3.4522 

23 

3.9571 

24 

3.9861 

25 

4.0513 

26 

4.3367 

27 

6.5503 

28 

8.0570 

29 

8.4335 

30 

8.8400 

31 

8.9907 

32 

10.2474 

33 

10.5133 

34 

11 .4797 

35 

11  .6510 

36 

11.6522 

37 

11.6523 

38 

11  .6606 

39 

11  .7353 

40 

13  .3386 

41 

14.0969 

42 

14.8900 

43 

16 .5092 

44 

17.1423 

45 

17.5072 

46 

17.7673 

47 

17 .7682 

48 

17  .7686 

49 

17.7767 

50 

21 .0746 

51 

21  .6642 

52 

22.0228 

53 

23 .4744 

54 

23 . 8904 

•  Actuators  -  Rate  Sensors 


Absolute  values  of  the  open-loop  nodal  costs  are  computed 
for  comparison  using  the  c o 1  or e d - no i s e  disturbance.  The  RMS 
second-order  modes  and  modal  costs  are  plotted  versus  mode  number 
in  Fig.  2.  Immediately  evident  is  the  clustering  of  these  modal 
phenomena.  The  disturbance  effect  as  seen  through  the  line-of- 
sight  is  constrained  to  clusters  of  modes  as  is  the  ability  to 
measure  and  control  the  model.  Alignment  of  the  "controllable 


actuator/ sensor  configuration  for  the  problem.  Table  2  gives  the 
quantitative  modal  ranking  prescribed  by  each  method. 


Control  Designs 


Three  controllers  are  designed  and  analyzed  for  the  model. 
The  first  uses  10  low-frequency  modes  from  the  model  and  rela¬ 
tively  low  gaina,  i.e.,  a  cautious  design.  Based  upon  closed- 
loop  analysis  of  this  controller  and  re-evaluation  of  Fig.  2,  a 
second  controller  is  designed.  The  bandwidth  of  this  controller 
is  allowed  to  expand  into  an  adjacent  "dead  zone"  of  the  model 


spillover.  A  low  authority  control  is  designed  and  added  to  the 
high-gain  controller  to  give  the  third  design.  Full  closed-loop 
modal  and  stochastic  analyses  are  presented. 
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Caatioii  Controller 


Mod «  Selection 

Using  the  nodal  analyses  of  Fig.  2,  nodes  1,  3,  4,  6,  7,  8. 
9,  15.  16  and  17  are  selected.  These  are  the  first  ten  nost 
distnrbable  nodes  with  respect  to  LOS.  Inspection  of  the  second 
colnnn  <  Table  2  shows  that  several  controllable  nodes  in  the 
design  lodel  bandwidth  have  been  onitted,  and  it  is  reasonable  to 
expect  that  any  nndesirable  in-band  controller  effects  will  be 
related  to  these  nodes. 

Control  Design 

The  ten  selected  nodes  and  six  disturbance  states  transforn 
to  a  26th  order  linear  design  nodal.  Using  standard  LQG  tech¬ 
niques,  a  state-feedback  control  is  found  which  nininixes 

J  «  lin  E{xRxR  +  bn^iiR)  (13) 

t-t-* 

where  z „  is  the  LOS  output  of  the  reduced  nodel  and  b  is  a 
R 

tuning  paraneter.  A  reasonable  closed-loop  design  spectrun,  shown 
in  Table  3,  is  obtained  using  b  =  10  **.  The  significant  effect 
of  this  control  is  on  nodes  7  and  17.  Mode  7  danping  increases  by 
nearly  a  factor  of  30,  but  its  frequency  increases  by  less  than 
one  percent.  Mode  17  experiences  a  danping  increase  of  a  factor 
of  200  and  a  frequency  Increase  of  21%,  placing  it  between  the 
frequencies  of  unnode  1  ed  nodes  18  and  19.  Referring  again  to 


Table  2  Modal  rankings  suggested  by  the  analysis 


DIST/LOS  ACT/LOS  ACT/SEN  MODAL  COST 


Mode 

°1 

Mode 

Mode 

c  . 
i 

Mode 

1/  L 

"i 

17 

0.0242 

17 

0.0366 

5 

0.1836 

17 

0.0010 

1 

0.0202 

1 

0.0130 

17 

0.1813 

16 

0.0003 

16 

0.0141 

7 

0.0119 

7 

0.1037 

1 

0.0002 

6 

0.0116 

34 

0.0081 

2 

0.0941 

6 

0.0001 

4 

0.0110 

27 

0.0080 

41 

0.0831 

7 

0.0001 

7 

0.0099 

3 

0.0078 

34 

0.0786 

8 

0.0001 

8 

0.0085 

16 

0.0074 

27 

0.0702 

4 

0.0001 

3 

0.0084 

39 

0.0072 

39 

0.0700 

32 

0.0001 

9 

0.0059 

41 

0.0066 

18 

0.0582 

15 

0.0000 

15 

0.0056 

40 

0.0051 

50 

0.0570 

34 

0.0000 

32 

0.0038 

6 

0.0042 

40 

0.0564 

9 

0.0000 

34 

0.0034 

28 

0.0042 

29 

0.0510 

3 

0.0000 

10 

0.0033 

50 

0.0041 

1 

0.0466 

27 

0.0000 

27 

0.0032 

18 

0.0039 

28 

0.0426 

39 

0.0000 

18 

0.0028 

8 

0.0038 

10 

0.0417 

31 

0.0000 

31 

0.0025 

31 

0.003  8 

53 

0.0401 

28 

0.0000 

39 

0.0024 

32 

0.0034 

31 

0.0379 

42 

0.0000 

28 

0.0024 

42 

0.0034 

42 

0.0367 

41 

0.0000 

Real  Imaginary 
part  part 

-4.8001  10.8240 
-0.0270  7.3328 
-0.0127  6.0956 
-0.0115  3.9841 
-0.0114  3.8621 
-0.2178  3.7528 
-0.0114  3.4864 
-0.0033  1.0962 
-0.0222  0.9388 
-0.0203  0.7341 


Frequency 

(Hz) 

Damping 

ratio 

1 . 8845 

.4054 

1.1671 

.0037 

0.9701 

.0021 

0.6341 

.0029 

0.6147 

.0030 

0.5983 

.0579 

0.5549 

.0033 

0.1745 

.0030 

0.1495 

.0237 

0.116  9 

.0276 

Real 


Iaiaginary  Frequency  Damping 


part 

part 

-5.5453 

11.223  9 

-0.5170 

7.3093 

0.0829 

6.0898 

0.0654 

3.9831 

0.2152 

3.9445 

0.4231 

3.6162 

0.0944 

3.6258 

0.0877 

1.0982 

0.0699 

0.9299 

0.1559 

0.7796 

(Hz) 

ratio 

1.9925 

.4430 

1.1662 

.0706 

0.9693 

.0136 

0.6340 

.0164 

0.6287 

.0545 

0.5795 

.1162 

0.5773 

.0260 

0.1753 

.07  96 

0.1484 

.07  50 

0.1265 

.1961 

/ 


Fig.  2,  modes  19,  20  and  21  are  apparently  uncontrollable  with 
respect  to  both  1 ine-of-sight  and  sensor  placement.  Mode  18, 
however,  has  some  controllability  (it  ranks  ISth  out  of  54).  The 
control  bandwidth  is  therefore  likely  to  influence  mode  18. 

A  standard  Kalman  filter  is  found  for  the  design  model  using 
the  actual  disturbance  PSD  and  LOS  "measurement  noise"  of 
intensity  5  x  10  which  gives  a  filter  bandwidth  roughly  equiv¬ 
alent  to  the  control  bandwidth.  The  filter  spectrum  is  shown  in 
Table  4.  Here  also,  the  frequency  increase  in  mode  17  is  suffi¬ 
cient  to  encompass  unmodeled  mode  18. 

Evaluation 

The  10-mode  controller  is  implemented  with  the  full  114th 
order  model.  Fig.  3  shows  that  the  closed-loop  spectrum  is 
stable,  and  also  that  the  dominant  control  effects  are  limited  to 
a  few  poles.  The  expanded  scale  in  Fig.  4  shows  the  anticipated 
spillover  in  mode  18.  Notice  also  the  decreased  damping  in  un¬ 
modeled  modes  2,  5,  and  10,  which  are  well  within  the  controller 
bandwidth.  This  "in-band  spillover"  occurs  despite  the  relatively 
low  controller  authority  over  these  three  modes  (see  Fig.  2). 

The  stochastic  performance  with  this  controller  is  shown  in 
Table  5.  Using  a  total  RSS  control  effort  of  12. 2N,  the  RSS  LOS 
error  is  reduced  to  16%  of  its  open-loop  value.  Fig.  5  shows  the 
absolute  RMS  modal  costs  with  and  without  this  controller. 
Significant  cost-level  reductions  occur  in  modes  17,  16  and  1, 
which  are  the  largest  open-loop  contributors.  Notice  that  the 
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OPEN-LOOP  FLEXIBLE  MODE  NUMBER 


Figure  5  . 


Closed-Loop  Model  Costs  for  10-Mode  Controller 


cost  contributions  of  "spillover  nodes"  2,  5,  10,  and  18  do  not 


|  increase. 

High-Gain  Controller 

i 

Node  Selection 

To  account  for  the  obaerved  spillover  effects  and  improve 
closed-loop  performance,  modes  2,  5  and  10,  and  high  frequency 
mode  18  are  added  to  the  control  design  model.  The  fact  that 
modes  19  through  26  are  highly  insensitive  to  control  and  distur¬ 
bance  inputs  (see  Fig.  2)  is  used  to  justify  a  closed-loop 
controller  bandwidth  extending  into  thia  dead  zone. 

Control  Design 

As  before,  a  standard  linear  LOS  regulator  and  Kalman  filter 
ia  designed.  For  a  control  penalty  of  5  x  10  I  and  'measurement 
noise"  of  5  x  10  **I,  the  controller  has  bandwidth  to  3.5  Hz  which 
ia  well  into  the  dead  zone. 

Eva luat i on 

The  closed-loop  spectrum  resulting  from  the  high-gain  control 
ia  ahown  in  Fig.  6.  Aa  before,  the  dominant  control  effects  are 
limited  to  a  few  modes.  The  expanded  scale  in  Fig.  7  shows  that 
the  inclusion  of  modes  2,  5,  and  10  in  the  design  model  eliminates 
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their  in-band  spillover.  The  expected  high-frequency  spillover  is 
present;  particolarly  in  node  27,  which  is  the  closest  unmodeled 
mode  to  the  controller  bandwidth  having  significant  control¬ 
lability  (see  Fig.  2).  Notice  that  as  expected  no  problems  occur 
in  the  dead  zone  modes.  Higher  frequency  spillover  effects  are 
evident  in  modes  53,  52  and  31,  which  are  well  above  the  bandwidth 
of  model  certainty  and  cannot  be  handled  explicitly.  Fig.  8  shows 
the  modal-costs  of  the  14-mode  control  compared  with  the  10-mode 
control.  Notice  the  overall  improvement  in  performance,  with  the 
exception  of  a  sp i 1 1 o v e r - i ndn c e d  cost  increase  in  mode  27.  The 
total  RSS  LOS  error  with  this  controller  is  1%  of  the  open-loop 
value  (see  Table  5). 

High-Gain  Controller  with  Low-Amthority  Control 

The  purpose  of  low-authority  control  (LAC)  is  to  add  damping 
to  structural  modes  in  the  high-frequency  uncertainty  region. 
This  method  is  described  fully  in  Ref.  11.  The  actuator  to  rate 
sensor  rankings  in  Fig.  2  indicate  potential  LAC  modal  authority. 
For  this  example,  an  LAC  is  designed  to  add  4  percent  damping  to 
the  structural  modes. 

Evaluati on 

Fig.  9  shows  the  characteristic  damping  increase  in  the 
closed-loop  spectrum.  The  expanded  scale  in  Fig.  10  and  the 
modal-costs  in  Fig.  8  show  that  all  spillover  effects,  including 


H AC /LAC  14-MODE  HAC  10-MODE  HAC 


.00010 

.00008 

.00006 

.00004 

.00002 

.00000 
,00010 

.00008 

.00006 

.00004 

.00002 

.00000 
,00010 

.OOOOB 

.00006 

,  .00004 

^  .00002 

-  ,00000 


I 

I 

I 


Table  3  Stookaatio  pirf oraaici  analyais  of  eoatrollara 


Open 

Loop 

10-Mode 

Control 

14-Mode 

Control 

14-Mode 

Control 

with  LAC 

RSS  LOS  z 

( |i  Rad) 

131.7 

62.4 

69.0 

26 .4 

RSS  LOS  y 

( n  Rad) 

1109.5 

175.9 

87.5 

75.6 

Total  RSS 

LOS(m  Rad) 

1117.3 

186.6 

111.4 

80.1 

RSS  Control 

Effort  (N) 

0 

12  .2 

28.6 

25  .3 
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